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EXECUTIVE SUMMARY 


For many years, water supply planning and restoration projects have used regional 
or high-resolution sub regional, numerical models for analysis in the LEC, The principal 
tool used by the SFWMD for the regional analysis of its water control system is the South 
Florida Water Management Model (SFWMM). The SFWMM simulates the groundwater 
system within its boundary using a vertically aggregated, single layer to emulate the 
composite effects of the nonhomogeneous, surficial aquifer. However, the grid cell size 
is quite large (2 miles by 2 miles) and is not well suited for evaluating sub regional scale 
projects. This model was originally developed by MacVicar and others (1984) and has 
undergone numerous revisions as more information about the surface water and 
groundwater system becomes available. 


Several high-resolution, sub regional groundwater flow models were also 
developed to evaluate potential benefits and impacts of proposed changes to the water 
management system. Models developed most recently include the Martin County 
Groundwater Model, North Palm Beach Groundwater Model, South Palm Beach 
Groundwater Model, Broward Groundwater Model, North Miami-Dade Groundwater 
Model, and South Miami-Dade Groundwater Model. These models use the United States 
Geological Survey (USGS) modular three-dimensional finite difference groundwater 
flow model, commonly known as MODFLOW. These models have previously been 
used to support numerous projects in the Restudy, Florida State Senate Bills, the Water 
Preserve Feasibility Study, the Lower East Coast Water Supply Plan, permitting of the 
FEMA C-4 Impoundment, permitting and analysis of the G-160 structure, among others. 


In order to simulate the majority of the LEC planning area, the sub regional 
models were modified, updated, and combined into one model, the Lower East Coast sub 
regional Model (LECsR). Several factors influenced the decision to develop the model, 
LECsR. First, the knowledge base relating to hydrology and hydrogeology of south 
Florida is growing continually, owing to extensive monitoring networks and cooperative 
efforts, such as the Restudy and CERP. Increased knowledge of the system has tied in 
with the development of several add-on packages to MODFLOW, which simulate 
complex system interactions involving water restrictions, wetland hydroperiods, and flow 
diversions. In conjunction with the system knowledge, advancements in computer 
technology have influenced the development, maintenance, and application of one, high- 
resolution model spanning a large area (approximately 7,500 square miles). Maintaining 
and applying one, large model, rather than six, county-specific models, promotes efficient 
use of human resources, as well as consistent methodology during model development. 


The sediments composing the Surficial aquifer system and the Biscayne aquifer 
are primary marine carbonate and clastic sediments of Pleistocene to Miocene age. The 
model was discretized into three layers utilizing Chronostratigraphic correlation of the 


DRAFT LEC subRegional MODFLOW Model Documentation 


sub-aerial exposure surfaces within these sediments resulting from the periodic 
submergence and emergence of the Floridan peninsula from fluctuations in sea level 
stands. 


The LECsR encompasses approximately 7500 square miles. It ranges from C- 
44/Stuart to the north, the Atlantic Ocean on the east, Lake Okeechobee and the western 
edge of the Water Conservation Areas and the West, and Biscayne Bay/Florida Bay on 
the south. A grid spacing of 704 foot by 704 foot was chosen to match the SFWMD 
regional models which can be used to provide internal boundary conditions for future 
simulations. The model was calibrated over an extend period of 14 years from January 
1986 through September of 1999 with daily time steps and stress periods. The model 
was primarily calibrated to observed heads with an overall mean error of 0.0 feet, a mean 
absolute error of 0.54 feet and root mean square error of 0.72 feet. The calibration results 
indicate a reasonable match between observed and measured water levels in most areas of 
the model domain. The model was verified from September, 1999 through December, 
2000 and produced similar results to the calibration period. The budget from the model 
indicates that recharge to the aquifer is the primary inflow to the system with 
evapotranspiration and canal drainage as the primary outflows. 


The model can provide an understanding of the movement of water in the study 
area. The model proved quite robust considering the long term daily calibration and 
verification periods in its overall performance during the simulation period. The 
calibration period included both a 1 in 100 year drought and a 1 in 100 year wet events 
and the model did not react adversely to these extreme events. The LECsR model 
conceptualization and discretization was designed at a subregional or basin level scale. 
The special variability of the input parameters is also best described at a similar scale. 
Therefore, the model should be used for regional to subregional or basin level projects 
and interpretation of the results should also be at that scale. In addition, the model does 
provide a reasonable estimate of drawdowns associated with wellfield withdrawals and 
the ground/surface water interactions within wetland systems. 
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CHAPTER 1 


Introduction 


BACKGROUND 


Development of the Lower East Coast Subregional (LECsR) Model was initiated 
by the South Florida Water Management District (SFWMD), Model Application Section 
of the Water Supply Department to support ongoing water supply management and 
ecosystem restoration efforts with computer modeling. This modeling project supports 
the development of several authorized activities including: 


* Lower East Coast (LEC) Regional Water Supply Plan (Section 
373.0361, F.S.), 


e Minimum Flows and Levels (MFLs) for the Loxahatchee River, 
Biscayne aquifer and other watercourses (Section 373.042, F.S.), 


° Core SFWMD projects for evaluating engineering designs, permitting 
and operational rules, 


e Water Reservations (Section 373.223(4), F.S.), 


° Several components in the Comprehensive Everglades Restoration 
Plan (CERP) (USACE and SFWMD 1999), and 


° Acceler8 (Bush and Gutierrez 2004). 


Water supply plans develop strategies (based on a 20-year future planning 
horizon) to meet future water demands of urban and agricultural uses, while meeting the 
needs of the environment by identifying historically used sources that will not be 
adequate to meet future demands and evaluating other source options to meet the deficit. 
Minimum flows and levels identify the groundwater levels in an aquifer and surface 
water levels and flows at which further withdrawals would significantly harm the 
resource. MFLs must be established for surface waters and aquifers within the boundaries 
of all water management districts, which are directed to use the best available data. 
Water reservations protect fish and wildlife or public health and safety while protecting 
existing legal uses of water by reserving water from consumptive uses in designated areas 
of concern. The Water Resources Development Act of 2000 (Public Law No. 106-541, of 
the 106th Congress) provides approval to restore the south Florida ecosystem in a plan 
called CERP, which consists of 60 components that include above- and under-ground 
water storage, preserve and treatment areas, decompartmentalization of natural areas, 
operational changes to improve ecology and water conservation benefits, and reuse of 
water. Acceler8 is an initiative by the State of Florida, SFWMD, and federal partners to 
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accelerate the completion of eight, identified projects from the CERP. By accelerating the 
completion of these projects, benefits from the Everglades restoration will be experienced 
sooner than originally scheduled. These projects require numerical tools that can quantify 
the availability of current water resources and in most cases, estimate or predict the 
demands on future water resources. 


The study area focuses on the greater part of the Lower East Coast Planning 
Region, which covers roughly 10,500 square miles, and consists of a large freshwater 
lake, wetlands and estuaries, uplands, agricultural and urban areas, and coastal 
ecosystems lying within a highly managed system of canals, operational structures, 
levees, and retention ponds (Figure 1). This area continues to experience population 
growth and a subsequent increase in overall water demands, which must be balanced with 
the water needs of the environment. 


PURPOSE AND SCOPE 


The purpose of this project is to develop a numerical model capable of simulating 
the groundwater flow of the Surficial Aquifer System (SAS), wetland hydroperiods, 
water deliveries, canal-aquifer interaction, and general management of the water 
resources for the Lower East Coast of Florida. The model will be used as an interpretive 
and predictive tool and must be flexible enough to answer a wide range of questions. 


The model will be used primarily to perform predictive simulations of proposed 
water resource projects and/or operational schemes. It will also be used as an interpretive 
tool for the SFWMD by identifying data gaps in aquifer characteristics, hydrogeologic, 
stratigraphic and hydrologic parameters, and producing water budgets and groundwater 
flow maps to better understand the surface/groundwater system. In addition, due to the 
extensive data collection and quality assurance and control efforts associated with model 
development, data sets can be used for a wide variety of purposes including populating 
various databases and providing information to the public according to Florida’s Public 
Records Laws (Section 119.07, F.S.). 


The model was developed to support several current and future projects in the 
LEC. The ability to evaluate water levels for water availability, cumulative wellfield 
impacts to natural areas or existing legal uses of water, water table increases that could 
potentially affect flood protection, and the frequency and severity of water shortages was 
considered when designing this numerical model. Water supply managers evaluate urban 
and agricultural water uses and must ensure current and future reasonable beneficial uses 
while protecting and restoring the environment and water resources. Due to these 
considerations, and to account for the complexities of the aquifer and drainage systems, 
the model must adequately represent horizontal and vertical aquifer heterogeneities and 
canal hydraulic properties. 


Draft - 2 


LEC subRegional MODFLOW Model Documentation DRAFT 















$ AQ Sumteq n | _seminoteV¥usia 
Hernando "A | Lake ! 
Sca wt A 
Pasco ke … 


a¢ =-=- 


l 
deles = 1 ARilIsborougH Polk 
d $ | 
M d SN 
/ | L l Kl y Indian River 
F — — — f -RRR | 1 
Lee ` KISSIMMEE | 
pS Manatee | Hardee l | BASIN | r 
Q | M ^L d L € 
EN vz em l L 'RESG TOG | > Luci 
| | —----W.,52:9 UPPER 
lm! Gel A BASTAS 
Sarasota De Soto ia ser "cl = 


E 


— Palm Beach - x e = 

i X LOWER 7 

A BASTI S]: 
= = GOASID E ATLANTE LC. 
TP UD RS ETE OCEAN 


OF MEXICO 


— 


| Kissimmee Basin e» South Florida Water 
Management District 
Boundary s 


ALI 
4 


LI 


— armm—mEm Miles | Lower East Coast 
Upper East Coast 0 5101520 subRegional Model 


c j q 8 $5223 


= = County Boundaries Prepared by : MAPU 
Date: 02/28/06 
Map Doc.: planreg_06.mxd 


Figure 1. Boundaries for the South Florida Water Management District Planning Regions. 
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Other considerations were made to accommodate ecosystem restoration goals and 
objectives. The ability to simulate hydropatterns in natural systems, water storage 
scenarios (e.g., siting a reservoir, delivering water to natural or urban areas as necessary), 
and impacts to adjacent water users was recognized when designing this model. In order 
to address these requirements, the model needed to take into account the vertical and 
horizontal flux components of the wetland-aquifer interaction, include a mechanism to 
manage and store water deliveries, quantify canal seepage from the aquifer, and represent 
partially or fully penetrating flow barriers. 


Moreover, this report documents the development of a three-dimensional, 
numerical model, which simulates transient groundwater flow and the full hydrologic 
cycle in wetlands (Ge, integrated surface water-groundwater) in the LEC Planning 
Region. This is a technical document to help familiarize engineers, hydrogeologists and 
project stakeholders and the public with the model and its potential applications. It is not 
intended to serve as a user’s manual for model applications. 


With these objectives in mind, the scope of this document covers the development 
of the model in its entirety. This chapter (Chapter 1) has introduced the purpose and 
scope of this study, while referencing previous modeling studies that supported SFWMD 
goals. Chapter 2 presents the development of the conceptual model. Available data is 
collected and assembled accordingly to define the hydrogeologic system. Simulating the 
flow system involves two aspects — code selection and model design - which are 
discussed in Chapter 3. Chapter 4 details the processes of model calibration, sensitivity 
analysis and verification. Conclusions and recommendations with respect to model 
capabilities, limitations, and future improvements of this modeling study are presented in 
Chapter 5. 


A standard protocol for model development requires completing the steps in 
Chapters 1, 2, 3 and 4 (Anderson and Woessner 1992). The modeling protocol begins by 
establishing the purpose of the model. Having a clear idea of what questions need to be 
answered will help determine the magnitude and complexity of the model. Building a 
conceptual model is the next step. This step involves collecting data and assembling data 
into a meaningful representation of the system. At this stage, the modeling team should 
know the location of the system boundaries, the water budget components, and the flow 
system. 


Selecting the mathematical and numerical models require that the governing 
equations are valid (for the processes described in the conceptual model) and that the 
computer code accurately solves those mathematical equations. Model design begins with 
the spatial discretization of the grid and also involves choosing the appropriate temporal 
discretization. This stage also includes setting initial parameter values for properties and 
system stresses, as well as defining initial and boundary conditions. 


After the initial design of the model, the model will be executed numerous times 


with the intention of reproducing field (or measured) water levels and flow rates. This 
step is called model calibration and the model is considered to be calibrated when the 
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simulated and measured (or historical) water levels and flows meet the specified 
calibration criteria. A sensitivity analysis is conducted to quantify and show the effects of 
uncertainty in the calibrated model. 


The last required step in the modeling protocol is model verification. A second set 
of field data is developed and the model is executed. If the model results meet the 
established calibration criteria, there will be greater confidence in the calibrated model. 


PREVIOUS SUBREGIONAL MODEL DEVELOPMENT 


For many years, water supply planning and restoration projects have used regional 
and/or high-resolution subregional numerical models for analysis in the LEC Planning 
Region. A regional model is defined herein as a numerical model with a low spatial 
resolution and large spatial extent. Typically, a regional model evaluates the system-wide 
water resources for planning and management of those resources. A subregional model is 
a numerical model with a high spatial resolution and moderate to large spatial extent. 
Subregional models can also be used for planning and management of regional water 
resources, though these models generally address county-wide or local concerns that a 
regional model is not designed to address. 


The principal tool currently used by the SFWMD for the regional analysis of its 
water control system is the South Florida Water Management Model (SFWMM), which 
is considered a regional model (i.e., grid resolution of 2 miles by 2 miles with daily 
stresses). This model simulates structure operations and couples surface water and 
groundwater systems within its boundary, but has a coarse spatial discretization. 
Groundwater is represented using a vertically aggregated, single layer to emulate the 
composite effects of the nonhomogeneous, SAS. This model was originally developed 
and documented by MacVicar and others (1984) and has undergone numerous revisions 
(SFWMD 2005) as more information about the regional system has become available. 


Several high-resolution (i.e., grid resolutions between 500 ft and 1320 ft with 
daily stresses), subregional groundwater flow models were developed to evaluate 
potential benefits and impacts of proposed changes to the water management system. The 
most recently developed groundwater models include those for Northern Palm Beach 
County (SFWMD 2001), South Palm Beach County (Nair et al. 2001), Broward County 
(Restrepo et al. 2001), North Miami-Dade County (Wilsnack ef al. 2000), and South 
Miami-Dade County (Restrepo el al. 2001). These models use the U.S. Geological 
Survey (USGS) modular three-dimensional finite difference groundwater flow model, 
commonly known as MODFLOW (Harbaugh and McDonald 1996; McDonald and 
Harbaugh 1988). Figure 2 depicts the boundaries of the SFWMM and subregional, 
county-specific models. Originally, the subregional models addressed county-level 
problems and the knowledge base was specific to each county and its water resources. 


These previous subregional models utilized internal, surface water boundary 
conditions from the South Florida Water Management Model for predictive simulations 
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in the Lower East Coast Water Supply Plan (SFWMD 2000) and for the Comprehensive 
Everglades Restoration Plan (USACE and SFWMD 1999). Boundary conditions 
obtained from the regional model include canal stages and structure flows. The coupling 
of the regional and subregional models was necessary to evaluate water conditions in the 
urban areas, including urban wetland systems, at a finer resolution. Due to the large 
number of wellfields close to natural areas and the saltwater interface near the coast, this 
finer resolution was required due to the inability of the regional model to adequately 
predict the cone of influence from the wellfield withdrawals. Conversely, the subregional 
models did not simulate the operation of the Central and Southern Florida Project 
(USACE and SFWMD 1999) including the operations of Lake Okeechobee (see Figure 
1). For the urban wetland systems, which are generally much smaller than the regional 
model cells, the subregional models helped address potential impacts as well as evaluate 
improved hydropatterns due to the changes in water control structure operations. 


In order to simulate the majority of the LEC Planning Region, the subregional 
models were modified, updated, and combined into one model, the Lower East Coast 
subRegional Model (LECsR). Several factors influenced the decision to develop the 
LECsR Model. First, the knowledge base relating to hydrology and hydrogeology of 
south Florida is growing continually, owing to extensive monitoring networks, 
cooperative efforts (e.g., Comprehensive Everglades Restoration Plan [USACE and 
SFWMD 1999]), and research partnerships with other agencies (e.g., FAU and USGS). 
With respect to the groundwater models, increased knowledge of the system has tied in 
with the development of several add-on packages to MODFLOW, which simulate 
complex system interactions involving water restrictions, wetland hydroperiods, overland 
flow and surface water-groundwater interaction in wetland systems and canals, and water 
deliveries (or diversions). In conjunction with the system knowledge, advancements in 
computer technology (see Chapter 3) have influenced the development, maintenance, 
and application of one, high-resolution model spanning a large area (approximately 
7,500 square miles). Maintaining and applying one large model (LECsR), rather than five 
county-specific models, promotes efficient use of human and computer resources as well 
as consistent methodology during model development and application. This model will 
be applied for both regional and basin-scale projects. 
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Figure 2. Locations of the Boundaries for the LEC County-Specific Subregional Models and the 
South Florida Water Management Model. 
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CHAPTER 2 


Model Conceptualization 


DESCRIPTION OF THE STUDY AREA 


The study area for this report includes the Lower East Coast of Florida as shown 
in Figure 3. This area encompasses southern Martin County, the majority of Palm Beach, 
Broward, and Miami-Dade Counties, and eastern portions of Collier and mainland 
Monroe Counties; however, the focus of the study area lies within Palm Beach, Broward, 
and Miami-Dade Counties. This area includes the urbanized coastal communities, 
Everglades Agricultural Area (EAA), Stormwater Treatment Areas (STAs), Water 
Conservation Areas (WCAs) and Everglades National Park (ENP). STAs cover 47,000 
acres of large constructed wetlands which are designed to reduce phosphorus levels 
entering the Everglades to 50 parts per billion (ppb). WCAs are made up of saw grass and 
island hammocks, which are divided into five, enclosed compartments and are regulated 
through a series of pumps, weirs and canals. Everglades National Park consists of 
preserved, subtropical wetlands with temperate and tropical plant communities and 
marine and estuarine environments. The Florida Keys and barrier islands are not included 
in the study area due to the limited and thin freshwater lens that exist on several 
Keys and islands, which are not connected to the primary Surficial aquifer resources of 
the area. 


The study area was further refined and delineated using surface water and 
groundwater boundaries in south Florida (Figures 3 and 4). The eastern boundary 
coincides with the Atlantic Ocean and follows the brackish Intracoastal Waterway which 
serves as a stable physical boundary. The St. Lucie Canal (C-44) from the St. Lucie 
Estuary west to Lake Okeechobee is the northern boundary and is a large, managed canal 
that is part of the regional, surface water management system. Lake Okeechobee, in the 
northwestern corner of the study area, is a managed lake and the primary surface water 
supplier for southeast Florida. Western surface water boundaries include primary water 
management district canals Oe, from north to south, L-24 and L-28) until reaching 
Everglades National Park. The southern model boundary is defined by water levels in 
ENP near the groundwater divide and along Florida Bay. 


The Surficial Aquifer System, which includes the Biscayne aquifer, forms a 
physical boundary. The extent of the Biscayne aquifer forms the primary physical, 
aquifer boundary in the majority of the study area (Figure 3). The SAS is composed of 
both the Water Table and Biscayne aquifers, which lies within, but also extend beyond 
the study area boundaries. 
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Surface Water and Groundwater Boundaries that Define the Study Area. 
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PHYSIOGRAPHIC REGIONS AND TOPOGRAPHY 


White (1970) delineated ten distinct physiographic regions in Florida. Five of 
these regions are intersected by the study area and shown in Figure 5. South Florida 18 
unique for the east coast of the United States in that the land extends to the edge of the 
Continental Shelf as a result of rapid deposition of marine limestones during high sea 
level stands. Due to the low relief in Florida, White (1970) used a combination of natural 
features to develop the regions. The primary regions include the Everglades in the west, 
the Atlantic Coastal Ridge along the Atlantic coastline, and the Eastern Valley to the 
north. The Everglades is a broad low area in the center of South Florida, flanked by 
relatively higher lands to the east, west and north. It stretches from Lake Okeechobee 
southward to Florida Bay. The Everglades soils are generally peat overlying limestone 
with elevations ranging from sea level to + 12 feet, National Geodetic Vertical Datum of 
1929 (NGVD). The Atlantic Coastal Ridge is located on the east adjacent to the Atlantic 
Ocean and is characterized by a series of beach ridge and dunes deposits laid down 
during previous high sea level stands. Elevations along the Atlantic Coastal Ridge range 
from sea level to over + 30 feet NGVD but generally are closer to +15 to + 20 feet 
NGVD. The northern part of the study area is situated in the Eastern Valley. This area 
has slightly higher elevations but it contains poorly drained isolated wetlands unlike the 
Everglades with its continuous wetlands in areas where it has not been developed. 
Smaller regions outside the study area include the Southern Slope, the Gulf Coastal 
Lowlands, the Florida Bay Mangrove Islands, the High Coral Keys, the Reticulate Coast 
Swamps, the Green Ridge and the Osceola Plain (Figure 5). Except for the coastal or 
beach ridges, the study area is relatively flat with a gradual reduction in elevation 
southward. 


Topographic relief and the nature of surface water features affect the distribution 
of recharge and discharge within the study area. The Lower East Coast of Florida has low 
elevation, typically less than + 25 feet above sea level. The bottom of Lake Okeechobee 
is approximately at sea level with the water level in the lake ranging from + 11 to + 18 ft 
NGVD. The land immediately surrounding the lake has an elevation of about + 20 to 
+25 ft NGVD. 


The topographic data was collected from a number of sources, listed in Table 1, 
including Light Detection and Ranging (LIDAR) data for the coastal areas of the model, 
the USGS High-Accuracy Elevation Data, as well as data from several other Digital 
Elevation Models developed for smaller regions of the study area (Hinton 2004), as 
illustrated in Figure 6. The USGS High Accuracy Elevation Data was used in Everglades 
National Park and Water Conservation Area 1. This data was collected on approximate 
50 foot centers. Data for the Everglades Agricultural Area was collected using transects 
on 500 foot centers. In North Palm Beach County, a 5 foot Digital Elevation Model was 
developed from data collected using LIDAR. A digital elevation model was also 
developed using LIDAR for coastal Miami Dade County on 25 foot intervals. The 
remaining areas used small high resolution studies to fill small gaps or the USGS 24,000 
Quad points which consist of various elevations and benchmarks assembled by the 
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USGS. LIDAR data was also used along the coast of Broward County and Martin 
County. If not already converted, all topography was converted to NGVD. 
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Figure 4. Physiographic Regions in the Lower East Coast Subregional Model (Modified from 
Florida Department of Environmental Protection). 
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Figure 5. Data Sources Used to Develop The Topography for the Lower East Coast 
Subregional Model. 
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Table 1. Description of Sources Used to Develop the Topography. 
1992 Florida Game and Fish 
Management Area grid unreported 
U.S. Geological Survey 24K 
USGS24K Quad Points (Topography Five varies varies +/-2.5 ft NGVD29 
Foot USGS 24k Points ) 
Survey of Water Conservation 
Combined Structure and 
Operational Plan for MWD and 
C-111 LIDAR Surveys Digital prepared bai 
Elevation Model 
SFWMD Subregional 
NGVD29 
unreported NAVD88 


Publication Beésolüfióh Vertical 
Date Accuracy Vertical Datum 
1986 AeroMetric Corporation 
Survey of the 8.5 Square Mile 1986 300 ft +/-4 in NGVD29 
Area 
Commission Survey of .5-minute 
Holeyland Wildlife latitude/longitude 
1992 Florida Game and Fish 
Commission Survey of 
ROTENBERGER Rotenberger Wildlife dree 
Management Area rid unreported 
EAA-DEM Eve ETS 2003 500 ft unreported NGVD29 
Digital Elevation Model 
1992-3 Keith and Schnars 
WCA2A Survey of Water Conservation 0.0003-0.0329 m NAVD88 
Area 2A 
1999 EarthData International 
1999 5m fl-15 cm NAVD88 
Area 3A 
NPB-DEM NOTI) A A 2004 B ft unreported NG VD29/NAVD88 
Elevation Model 
Groundwater Modeling 
Topography 
Evaluation Data Collection 
+/-0.7 ft NAVD88 
NOAA Soundings from I I 
NOAA Everglades National Park Staff varies varies MLLW, NGVD29 





The study area has a tropical to subtropical climate (Pendleton, Dollar and Law 
1976). The summers are warm and relatively wet with mild, relatively dry winters. 
Temperatures are moderated by the close proximity of the Gulf Stream, which is located 
Just off the south Florida coastline. The Gulf Stream originates in the Gulf of Mexico and 
flows around the tip of Florida, heading north into the northeastern Atlantic Ocean. 
Maximum temperatures rarely exceed 100 degrees Fahrenheit (F) and minimum 
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temperatures are rarely below freezing. The maximum humidity occurs by dawn while 
the minimum humidity usually occurs in the afternoon. The relative humidity in Florida 
generally exceeds 50 percent all year. 


Data for rainfall and temperature were collected from the National Oceanic and 
Atmospheric Administration (NOAA) stations (Figure 7) and from the extensive 
SFWMD data base DBHYDRO. Rainfall stations were selected based on period of 
record, data quality, and spatial distribution. Stations with a long-term record were given 
preference. Published data from NOAA was preferred due to quality control. Twenty-six 
rainfall stations were selected for the study area. Similar criteria were used for selecting 
temperature stations which are used to determine reference evapotranspiration [ET] 
(SFWMD 2005). Unfortunately, long-term weather data are scarce in south Florida; most 
stations do not collect data for all the parameters required for a robust calculation of 
actual ET. 
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Figure 6. Locations of Rainfall and Temperature Stations. Temperature Stations are used to 
Develop Reference ET. 
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Rainfall 


Rainfall represents the largest input of water into the hydrologic system. The 
average annual rainfall throughout the study area is 54 inches (Figure 7). The annual 
variation is high. Historically, at the West Palm Beach Airport, the wettest year on record 
is 1947 with 108 inches and the driest year is 1954 with 37 inches recorded. Over the last 
twenty years, the wettest year, 1994 had 86 inches but the driest year, 1989 had only 39 
inches. Due to these extremes in annual rainfall, the region periodically undergoes 
flooding and prolonged drought. 


Rainfall was summarized based on twenty-six rainfall stations over the years, 
1965-2000. Rainfall, on an average annual basis is relatively consistent throughout south 
Florida; however, annual rainfall tends to generally decrease inland. Separate 
microclimates develop mainly associated with the tropical sea-breeze front, Everglades 
System, and Lake Okeechobee. Using 1992 as a typical year, average annual rainfall for 
coastal stations is 62 inches, which is above the annual rainfall of 58 inches for the inland 
stations. 





— 
= 
— 
8 
E 
e 
a 
E 
3 
E 
c 
< 
v 
O) 
G 
L-I 
> 
< 








Figure 7. Average Annual Rainfall (in) from 1965 -2000. 
On a daily or even weekly basis, the distribution of rainfall can vary dramatically 


across the region. This is especially true during the summer months when warming tends 
to create scattered afternoon thunderstorms referred to as convective rainfall. 
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Approximately 75 percent of the rain falls during the wet season (May through October) 
(Figure 9). Rainfall is provided by thunderstorms, tropical storms, and hurricanes in the 
summertime and passing cold fronts during the winter months. 
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Figure 8. Average Annual Rainfall (in) from 1965 -2000 by Month. 


Evapotranspiration 


Temperature stations (Figure 7) were used by the SFWMD to develop a regional 
reference ET data set based on the Simple Method (Abtew 1996). This method was 
selected to provide estimates of long-term historical (1965-2000) wet marsh potential ET 
for predictive simulations (SFWMD 2005) and will be discussed in Chapter 3. Reference 
ET is defined as the rate of ET from a hypothetical crop with an assumed height, a fixed 
canopy resistance, and albedo, which would resemble evapotranspiration from an 
extensive surface of green grass cover of uniform height, actively growing, completely 
shading the ground and not short of water (FAO 1990). The reference crop should be 
taken as a hypothetical crop with fixed parameters and resistance coefficients. These crop 
coefficients are complex empirical factors derived from experimental data and encompass 
all characteristics of the crop which differ from those of the reference crop. These 
coefficients quantify how soil and crop conditions affect actual ET. 


Evapotranspiration represents the largest water loss from the system and was 
summarized based on seventeen NOAA stations over the years, 1965-2000. In general, 
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reference evapotranspiration is on the order of 55 in to 59 inches per year (Figure 10) for 
a wet marsh crop. In south Florida, it is recommended that grass be used as the reference 
crop, since previous studies have been conducted utilizing grass as a reference crop. 
Consequently, all the crop coefficients are based on grass. However, to ensure 
consistency between the regional and subregional modeling efforts, this study will use the 
data set that applies the wet marsh crop (herein referred to as reference crop) developed 
and used by the regional models. 
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Figure 9. Average Annual Evapotranspiration (in) for Wet Marsh from 1965 -2000. 


Actual ET depends on the availability of water to the crop and most often is less 
than the potential ET which has no constraints on available water. Although not as varied 
as rainfall, actual evapotranspiration rates in south Florida can be greater than 40 in/yr 
and can approach and even exceed rainfall rates during drought events (USGS 1996). 


The monthly distribution of reference ET (Figure 11) follows a bell-shaped 
curve. Over the thirty-six year period of record, the month with maximum ET rates is 
May. Reference ET is lowest during the month of December. In general, the distribution 
of reference ET follows the solar radiation curve. 
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Figure 10. Average Annual Evapotranspiration (in) for Wet Marsh from 1965 -2000 by Month. 


LAND USE 


At the time of data collection, available land use and land cover data included 
1988, 1995, 2000, and future data sets. The 1988 land use data were produced by photo 
interpretation of color infrared photography from the late 1980's and the accuracy of 
parameters has not been assessed. The 1995 land use data 19 a composite of vegetation 
types from a variety of sources and was linked to crosswalk tables that associate land use 
codes with the vegetation. The 2000 and future land use analysis used photo 
interpretation and Palm Beach, Broward, and Miami-Dade Counties’ comprehensive 
water management plans to estimate the developable areas. 


Land use was summarized herein using 1995 data. Land use classification 
includes several levels which give details about the type of use. Figure 12 shows a 
partially aggregated level 3 classification (see Table 2). Wetlands cover half of the area 
and include large, expansive systems like the Everglades. Northern Palm Beach County 
contains several wetlands including those in Jonathan Dickinson State Park, Palmar, 
Corbett Wildlife Management Area, and Loxahatchee Slough. The second largest land 
use consists of urban and built up areas. The majority of these areas are generally within 
the boundaries of the LEC Service Areas; however, there are urbanized areas in the Lake 
Okeechobee Service Area. The largest percent of agriculture occurs in the EAA and is 
made up of mostly sugarcane. 
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Table 2. Percent of Total Area by Land Use Types in the Study Area. 


Land Use Type______| Acres | Percent | 
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Figure 11. 1995 Land Use / Land Cover. 
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SURFACE WATER CHARACTERISTICS 


Surface waters within the model area include Lake Okeechobee, numerous 
smaller lakes, retention ponds, limestone mining pits, a large wetland ecosystem 
(Everglades), numerous smaller wetland systems, the Wild and Scenic Loxahatchee River 
and its tributaries, numerous canals, coastal estuaries including Lake Worth Lagoon, 
Biscayne Bay and Florida Bay, and the Atlantic Ocean. 


Surface Water Management Systems 


There are three types of canal systems (or networks) — primary, secondary, and 
tertiary - in the SFWMD. Primary canals consist of canals and/or natural water courses 
providing final conveyance of overall drainage basin flows to the ocean or major inland 
water bodies. These canals are managed and operated by the SFWMD and, in certain 
cases, the U.S. Army Corps of Engineers (USACE). Secondary canals (e.g., drainage and 
water control district canals) are generally designed to control surface and groundwater 
elevations and maintain the quantity and quality of developed area runoff at pre- 
development levels. These canals typically discharge into the primary canal system or 
other natural waterways. Tertiary canals include swales, ditches, and retention ponds and 
are designed to remove stormwater in areas that are subject to inundation or provide 
irrigation to agricultural operations. 


Water management in south Florida is controlled by the Central and Southern 
Florida Project originally constructed by USACE and maintained and mainly managed by 
the SFWMD. The project is divided into four regions which are hydrologically 
interconnected and include: 1) the Kissimmee River — Istokpoga Basin; 2) Lake 
Okeechobee and the Everglades Agricultural Area; 3) the Water Conservation Areas and 
Everglades National Park and 4) the East Coast Canals (USACE 1999). The Kissimmee 
River-Istokpoga Basin lies north of Lake Okeechobee and provides inflow to Lake 
Okeechobee; however, this basin lies outside the model area so a detailed discussion is 
not included. 


Lake Okeechobee and the Everglades Agricultural Area 


Lake Okeechobee is the largest freshwater lake in the State of Florida, and is one 
of the largest freshwater lakes in the United States. It is regulated to meet numerous uses 
including protection of fish and wildlife, flood control for surrounding agricultural and 
urban lands, water supply for the Everglades Agricultural Area, the environment and the 
coastal users, protection of the dikes from hurricane damage, minimum flows and levels 
demands from coastal estuaries and environmental systems, protection of coastal aquifers 
from salt water intrusion, recreational users, and waterway navigation across the lake. 


Lake Okeechobee is regulated through a series of pumps and canals and is 


surrounded by approximately 1,000 miles of levees which range in height from 
approximately + 32 ft NGVD to + 45 ft NGVD. These pumps and canals maintain stages 
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in Lake Okeechobee within a seasonally changing band of high and low levels. Water 
levels in the lake are drawn down at the beginning of June to prepare for the hurricane 
and wet seasons. The target for this date is + 13.5 ft NGVD. Releases are made if the 
water level is higher than + 13.5 ft NGVD. Water levels lower than + 10.5 ft NGVD will 
trigger water restrictions from the lake. Releases are designed to protect the levees 
during a hurricane from wind-driven waves and tides, and also increases flood control for 
lands which discharge into the lake. Starting in October, water levels in the lake are 
allowed to rise to provide increased storage for water supply demands from users during 
the dry season. Releases will be made if water levels exceed the target of +15.5 ft NGVD. 
Water levels lower than +13.0 ft NGVD will trigger water restrictions from the lake 
during this time. 


Demands on Lake Okeechobee to supply water come from a variety of sources 
and occur at any time of year but demands are highest during the dry season. The 
Everglades Agricultural Area located to the east and south of the lake is a large 
agricultural area which relies on the lake almost exclusively for its source of water. The 
communities that surround the lake also rely on the lake for public water supply; 
however, several of these communities are looking to groundwater from the Floridan 
aquifer system to meet their future demands. In addition, coastal canals are maintained 
with lake water, when other sources are depleted, to prevent saltwater intrusion, provide 
water to coastal agricultural operations and recharge coastal wellfields and lake systems. 
The lake also supplies water to environmental features, when local supplies are depleted, 
including Everglades National Park and coastal estuaries. 


The Water Conservation Areas and Everglades National Park 


The Water Conservation Areas are a vast expanse of saw grass and island 
hammocks which are divided into five compartments which are surrounded by levees 
with the exception of a small gap on the west side of Water Conservation Area 3A (see 
Figure 4). The southernmost compartment, Water Conservation Area 3A, is adjacent to 
the Everglades National Park, which is the only subtropical National Park in the United 
States. Similar to Lake Okeechobee, the Water Conservation Areas are regulated to meet 
numerous uses including protection of fish and wildlife, water supply for the Everglades 
National Park, water supply for coastal users and wellfields, protection of coastal aquifers 
from saltwater intrusion and recreational uses. 


The Water Conservation Areas are regulated through a series of pumps, weirs and 
canals. These pumps and canals maintain stages in the Water Conservation Areas within 
a seasonally changing band of high and low levels. Water levels are drawn down at the 
beginning of June, which is the start of the wet season. During the fall, water levels are 
allowed to rise to provide increased storage for water supply demands from users during 
the dry season. Inflows to the WCAs are from local rainfall and runoff from the 
Everglades Agricultural Area to the north which is starting to be filtered (to remove 
nutrients) through several storm water treatment systems. During times of extreme 
drought, when the WCAs are below there minimum level, water is brought into the 
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WCAs from Lake Okeechobee, if available, and delivered to Everglades National Park or 
the coastal users. 


East Coast Canals 


The East Coast Canals are a series of canals which stretch approximately 170 
miles from St. Lucie County to the north, southward through Martin, Palm Beach, 
Broward and Miami Dade Counties ending at Florida Bay. This area consists of 
numerous canals, spillways, structures, culverts and pump stations. Water levels are 
maintained in the canals at specific levels to provide water supply and improve flood 
protection for coastal urban and agricultural interests, but are also controlled at high 
enough levels to stabilize the saltwater interface in the Surficial and Biscayne aquifers. 
During times of drought, water is brought into these canals from either the Water 
Conservation Areas or Lake Okeechobee to maintain the minimum levels, recharge local 
wellfields and irrigate landscaping and agricultural operations by moving water from the 
primary to the secondary canal systems. 


The East Coast Canals tie into secondary canal systems. The secondary canal 
systems were not constructed as part of the Central and Southern Florida Project. Runoff 
from the urban areas is discharged through the secondary canals into the primary canals, 
which ultimately release the water to tide for flood protection. During the dry season, or 
in some cases, during periods of low rainfall, secondary canal systems pump water from 
the WCAs or East Coast Canals to maintain minimum levels. Other times, secondary 
systems close off the connection to the East Coast Canals to avoid free canal flow from 
the secondary canals. 


Figure 13 shows the spatial variation in the primary and secondary canal systems 
due to different rules of operations across the study area. The control levels in the 
secondary canal system are primarily based on dry and wet season criteria and have a 
large impact on the Surficial Aquifer System by trying to maintain the control levels. 
The control levels for these secondary canal systems are determined by the local entity 
operating the system. Therefore, there are several hundred different control elevations 
for these canal systems. In general, the control elevations decrease from the north to the 
south and towards the coast. 
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Figure 12. Spatial Variation of Control Levels in Primary, Secondary, and Certain Tertiary Canal 
Systems. Each Color is Associated with an Operational Rule and/or a Structure. 
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Regional Flow System Operations 


Lake Okeechobee is the primary storage center of the hydrologic system in the 
LEC; however, its origin is really in the upper Kissimmee Chain of Lakes area in central 
Florida. Lakes Myrtle, Alligator, Mary Jane, Gentry, East Tohopekaliga, Tohopekaliga, 
and Kissimmee are the principal sources of inflow to Lake Okeechobee in the upper 
Kissimmee watershed (Figure 14). These lakes flow into the Kissimmee River through 
the S-65 structure. The Kissimmee River contributes about 69 percent of the inflow to 
Lake Okeechobee through the S-65E structure at the north end of the lake (Abtew et al. 
2002). Other major inflow sources to Lake Okeechobee are Lake Istokpoga (through S- 
68), direct rainfall, Fisheating Creek, the Taylor Creek-Nubbin Slough Basin, reverse 
flow from the Caloosahatchee River, the St. Lucie Canal, and back pumping from the 
Everglades Agricultural Area. 


Prior to the completion of the North New River, Hillsboro, Miami, and West Palm 
Beach canals in the early part of the 20th century, water would overflow from Lake 
Okeechobee and slowly move south through the Everglades wetland system (Sklar et al. 
2002). Discharges from Lake Okeechobee are now regulated by SFWMD by means of a 
series of canals and control structures. Outflows from Lake Okeechobee into the LEC 
occur primarily through structures on the south and southeast side of the lake. Lake levels 
are managed seasonally to provide flood protection and to meet urban and agricultural 
water supply demands. During wet periods, seepage from the Everglades is generally 
more than adequate to maintain water levels in the coastal aquifers; however, releases 
through coastal canals may be required to maintain regulation schedules in natural 
storage areas, such as Lake Okeechobee and the WCAs, and to provide flood protection 
(Figure 15). Alternatively, during dry periods additional water may be released from the 
regional system through the coastal canals to help recharge the Surficial Aquifer System 
in the coastal basins (Figure 16); these water supply releases can be triggered by a 
decline in canal water levels or by movement of the saltwater front inland (SFWMD 
2000). 


The four major canals connected to Lake Okeechobee are, from north to south, 
the West Palm Beach, Hillsboro, North New River, and Miami canals. They route water 
from Lake Okeechobee through the Everglades Agricultural Area to the south and into 
Water Conservation Areas (WCAs) 1, 2, and 3 (Figures 15 and 16). The WCAs are 
remnants of the historical Everglades, compartmentalized by canals, levees, and control 
structures into five separate reservoirs — WCA-1 (known as the Arthur R. Marshall 
Loxahatchee National Wildlife Refuge), WCA-2A. WCA-2B, WCA-3A, and WCA-3B. 
The WCAs retain excess floodwaters from Lake Okeechobee, the EAA, and the 
Everglades, as well as maintaining SAS groundwater levels, reducing saltwater intrusion 
into coastal basins, and benefiting fish and wildlife (SFWMD 2000). 
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Figure 13. Lakes and Rivers that are Inflow Sources to Lake Okeechobee. 
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Surface water enters Everglades National Park through the S-12 structures and 
other structures at the south end of WCA 3 (along the Tamiami Trail), eventually 
discharging to Florida Bay. A series of canals and structures on the eastern side of the 
WCAs maintain water levels in the coastal watershed. Coastal canals are regulated to 
provide flood control, to maintain groundwater levels for municipal wellfields, and to 
prevent saltwater intrusion. 


The opening and closing of the gates on the majority of the operating structures 
on the canals are controlled remotely from the South Florida Water Management District 
headquarters. By incorporating the ability to remotely operate the structures with 
weather radar and satellite images, the District has the ability to quickly drop the water 
levels in the canals, which improves the existing flood protection for that area. 
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Figure 14. Water Conveyance in the LEC Regional System during Wet Periods. Arrows Indicate 
Direction of Pumpage or Flow. 
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Figure 15. Water Conveyance in the LEC Regional System during Dry Periods. Arrows Indicate 
Direction of Pumpage or Flow. 
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Surface Water and Groundwater Interaction 


Groundwater levels in the LEC are closely related to surface water levels due to 
the high hydraulic conductivity of the SAS. The interconnection of surface water and 
groundwater levels has also been demonstrated from analysis of stable isotope data, 
which shows that the predominant source of groundwater in a study done in the north- 
central Everglades (i.e., northern WCAs) was from surface water recharge due to water 
control operations, historical changes in the water table and land subsidence (Harvey ef 
al. 2002). Recharge refers to flow from surface water to groundwater. Conversely, flow 
from groundwater to surface water is known as discharge. Harvey et al. 2004 found that 
both processes occur in the Everglades when water is exchanged across the wetland 
ground surface. 


Because the  predrainage Everglades ecosystem covered a large, 
uncompartmentalized area with free-flowing water, and with a small change in 
topography from north to south, the recharge and discharge rates were small (Harvey et 
al. 2004). Currently, the compartmentalization of the WCAs impounds the water and 
causes larger surface water-groundwater gradients than in the predrainage system. 
Surface water levels in WCA1 are higher than those in WCA2. Near the levees, 
groundwater from WCAI discharges into WCA2 due to the hydraulic gradient. Likewise, 
groundwater discharges from WCA2 into WCA3 as discussed above. Seepage from the 
WC AS can also occur west to east as underflow from the WCAs to the protective levees 
that separate the natural and urban areas. These discharges are intercepted by wellfields 
located in the eastern, urbanized areas. Harvey et al. 2004 recognized that there is little 
understanding of the source areas, flow paths, and travel times required for Everglades 
surface water to reach the eastern wellfields and also of how groundwater discharges 
sustain wetlands during drought. Moreover, the residence times of particles traveling 
through the shallow, peat sediments will be shorter (i.e., hours to days) compared to 
particles that travel deeper into the aquifer and emerge at a much later date (i.e., decades 
to millennia). 


Additionally, the SAS and canal system are well connected in the LEC due to the 
proximity of the canals to the water table and the permeability of the SAS, especially in 
the Biscayne aquifer. The main purpose of the canal system is to provide flood protection 
and prevent saltwater intrusion near the coast. Historically, increased pumping lowered 
water levels and groundwater recharged the canal system. However, saltwater intrusion is 
mitigated by holding surface water levels higher near the coast, consequently raising 
adjacent groundwater levels. Initially developed for flood control, the secondary canal 
structures have allowed the canal network to be used increasingly to regulate ground- 
water levels by draining the land during wet events and raising canal levels with water 
from the regional system during droughts. This secondary use has become as important as 
flood control in water-management practices (Miller 1998). Tertiary canals will provide 
storage and conveyance for local basin runoff. These canals can be a main source of 
local recharge into the aquifer system. 
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HYDROSTRATIGRAPHY AND HYDROGEOLOGY OF SOUTH 
FLORIDA 


Geologic Setting 


The basement complex of south Florida is composed of rhyolites and basalts 
believed to have formed in the vicinity of a Mesozoic hot spot (Heatherington and 
Mueller 1997). Overlying the volcanic basement rocks is an approximately 20,000-foot 
thick sequence of carbonates, evaporites and clastic sediments deposited continuously 
and ranging in age from the Triassic period to the present. Parker et al. (1955) and 
Randazzo (1997) have presented a generalized discussion of the pre-Pleistocene deposits 
of south Florida. This thick sequence of clastic and carbonate deposits is generally devoid 
of significant faulting or folding in southern Florida, suggesting a stable tectonic 
platform. Perkins (1977) has shown that the paleotopographic features of south Florida 
were inherited from older structural features with no indication of tectonic activity in the 
region throughout the Pleistocene epoch. In addition, there is no evidence that south 
Florida was subject to glacial loading during the Pleistocene. Thus, the deposition of the 
marine units in south Florida must be regulated by eustatic sea level fluctuations 
associated with glacial and interglacial stages of the Pleistocene rather than tectonic 
activity or crustal rebound. Figure 17 illustrates the lithologic units encountered in the 
study area. 
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Figure 16. Lithostratigraphic and Geohydrologic Units of the Surficial Aquifer System in 
Southeast Florida. Source: Adapted from Reese and Cunningham 2000; Perkins 
1977. 
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Quaternary Deposits 


Holocene Deposits 


The Holocene deposits are the most recently laid sediments and overlay the 
Anastasia Formation and the uppermost member of the Fort Thompson Formation. These 
recent deposits of muck, windblown sand, and peat (Parker and Cooke 1944) cover most 
of the model area from Lake Okeechobee south through the Everglades. Underlying these 
sediments is the Lake Flirt Marl Formation which is widely distributed throughout the 
model area. It consists of gray calcareous mud with shells of freshwater snails (Helisoma) 
and ranges in thickness from a few inches to 3 feet. Lake Flirt Marl underlies the peaty 
accumulations of the Everglades and overlies the Pamlico Sand formation and other 
Pleistocene or Pliocene deposits (Cooke 1945). 


Pleistocene Deposits 


Pleistocene deposits in south Florida are dominated by shallow-water marine 
carbonates and clastics. Interbedded within these units are indications of sub aerial 
exposure that include laminated crusts and freshwater limestones (Perkins 1977). Rapid 
facies changes across the region have resulted in the establishment of several distinct 
formations even though they appear to have been deposited contemporaneously. The 
stratigraphy of the Pleistocene sediments utilized in this report follow Parker et al. 
(1955), Brooks (1974) and Perkins (1977). The principal Pleistocene deposits within the 
study area include the Pamlico Sand, Miami Limestone, Fort Thompson Formation, 
Anastasia Formation, Caloosahatchee Formation and Key Largo Formation. Each of 
these deposits may be overlain by recent Holocene deposits. 


Pamlico Sand 


The Pamlico Sand formation was named by Stephenson (1912) to describe the 
fine sandy loams, sands and clays from Pamlico Sound in North Carolina. Parker and 
Cooke (1944) extended the Pamlico formation into Florida to include all marine 
Pleistocene deposits younger than the Anastasia Formation. 


The Pamlico Sand Formation unconformably overlies the Miami Limestone, 
Anastasia Formation, and Fort Thompson Formation and covers most of the region south 
of Lake Okeechobee. The formation does not usually extend above the 25-foot contour, 
which was the approximate shore line location of the third interglacial stage (Cooke 
1952). The formation is primarily composed of quartz sand with some local bodies of 
clay. The Lake Flirt Marl and recent deposits of muck, peat, and windblown sand overlie 
the Pamlico Sand (Parker and Cooke 1944). 
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Miami Limestone 


The Miami Oolite was studied in the mid to late 1800s by Tuomey (1851), 
L. Agassiz (1852), Shaler (1890), A. Agassiz (1895, 1896) and Griswold (1896). In 1909, 
Sanford proposed the name Miami Oolite for the oolitic deposits in the vicinity of Miami. 
However, Sanford (1909) separated the mainland oolite from the oolite found in the 
lower Keys. Later, Cooke and Mossom (1929) redefined the Miami Oolite to include all 
oolitic deposits in south Florida, including those found in the Florida Bay and the Florida 
Keys. Hoffmeister et al. (1967) suggested that the name be changed to the Miami 
Limestone based upon two distinct members, an upper oolitic facies and a lower 
Schizolle floridana bryozoan facies. 


The Miami Limestone unconformably overlies the Fort Thompson Formation 
(Perkins 1969) and is generally less than 20-feet thick but can exceed 40-feet thick along 
the coast in Miami. The oolites were deposited in a shallow sea similar to present 
conditions along the Bahama Bank, and this deposit has been dated at 125,000 years BP 
by Osmond ef al. (1965) and Broecker and Thurber (1965). The unit encompasses a 
broad arc shape from Boca Raton southward through Miami and terminates in the 
vicinity of the Key West/Dry Tortugas area. 


Fort Thompson Formation 


As originally proposed by Sellards (1919), the Fort Thompson Formation consists 
of alternating beds of marine, brackish and freshwater limestones located at the type 
location along the Caloosahatchee River. Although Sellards (1919) originally called this 
outcrop the Fort Thompson beds, Cooke and Mossom (1929) gave it formational status. 
The unit overlies the Caloosahatchee Marl in the northwest portion of the study area but 
lies directly above the Tamiami Formation through the vast majority of the study area., 
and 1s Pleistocene in age (Cooke and Parker 1944; Mitterer 1975). The formation inter- 
fingers with the Anastasia and Key Largo Formations on the eastern and southern 
boundaries. 


The Formation may be subdivided into four distinct members. The uppermost 
unit, termed the Coffee Mill Hammock member, is a marine unit characterized by 
abundant Chione cancellata mollusks. The unit grades into clastic beach and bar deposits 
of the Anastasia Formation eastward and into the carbonate oolitic shoal deposits of the 
Miami Limestone southward. Beneath the Coffee Mill Hammock member is Unit 1 of 
Brooks (1974). This unit is characterized by mollusk packstones and wackestones, 
including Chione cancellata. The clastic content decreases southward, transitioning from 
a shallow marine bay into an open marine platform environment (Perkins 1977). The 
lower middle unit, Unit 2 of Brooks (1974), is similar to Unit 1 but contains a very 
pronounced discontinuity surface between the units. Arenaceous, mollusk-fragment 
packstones characterize this unit, suggesting a restricted marine bay environment 
(Perkins 1977). The lowermost unit of the Fort Thompson Formation, Unit 3 of Brooks, 
has an increased clastic content compared to the overlying units, with mollusk-fragment 
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packstones and quartz sandstones common (Perkins 1977). The depositional environment 
appears to suggest an extensive somewhat restrictive marine bay (Perkins 1977). 


The depositional environment of the Fort Thompson Formation can be linked to 
late Quaternary sea level fluctuations. Parker and Cooke (1944) identified five separate 
marine invasions in this formation punctuated by freshwater deposition or subaerial 
exposure between the marine transgressions. Amino acid racemization work conducted 
by Mitterer (1975) revealed that three of the beds are late Pleistocene in age, ranging 
from 134,000 years before present (BP) to 324,000 BP. The age of the lowermost unit 
could not be determined using Mitterer's (1975) methodology. 


The Fort Thompson Formation covers the largest geographic expanse of all 
Quaternary formations in south Florida. The unit is generally less than 10-feet thick at the 
Caloosahatchee River but increases to over 50-feet thick in eastern Miami-Dade, 
Broward and Palm Beach Counties, where it makes up the highly productive zone of the 
Biscayne aquifer (Fish and Stewart 1989). This formation can be found in all areas of 
south Florida with the possible exception of the coast in Palm Beach County where it 
inter-fingers with the lower members of the Anastasia Formation. 


Anastasia Formation 


Originally named by Sellards (1912) for a series of coquina outcrops on Anastasia 
Island off the east coast of north Florida, the Anastasia Formation was extended by 
Cooke and Mossom (1929) to include all marine deposits located along the coast from 
Palm Beach County northward. The deposit consists of alternating offshore bar, beach 
ridge and dune system environments and may be at least 50 feet thick along the coast 
(Perkins 1977). The age of the formation is estimated to be late Pleistocene and 19 
considered to be deposited contemporaneously with the Miami Limestone and Fort 
Thompson Formation (Parker et al. 1955). The formation general lies directly upon the 
Tamiami Formation. 


The Anastasia Formation can be divided into two distinct facies, a coquina facies 
and a shell rock facies (Lovejoy 1983). The coquina facies represents a high energy 
environment typical of an offshore bar complex and is generally aligned with the present 
coastline. Behind this bar complex is the shell rock facies characterized by a diverse 
molluscan fauna with minimal damage to the fossils in a fine-grained quartz matrix, 
suggesting a shallow bay origin (Lovejoy 1983). 


Caloosahatchee Formation 


The Caloosahatchee beds or marls were considered the lower shell beds along the 
upper portion of the Caloosahatchee River by Dall (1897). The formation was officially 
named the Caloosahatchee Marl by Matson and Clap (1909). The age of the formation is 
believed to be late Pliocene to early Pleistocene. The formation occurs from Saint 
Petersburg on the west coast, southeast past the Caloosahatchee River, to the east coast 
(Randazzo, 1997). The formation generally thickens toward the east coast (Parker and 
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Cooke, 1944). Within Palm Beach County, however, the unit is restricted to the western 
part of the model area where it may directly underlie the Fort Thompson, Anastasia or 
Holocene sediments. The Formation is restricted to the northwestern portion of the study 
area in the vicinity of Lake Okeechobee. It consists of sandy, shelly marls, silt and clay 
with occasional stringers of well-consolidated sandy limestone. 


Key Largo Limestone 


Sanford (1909) proposed the name Key Largo Limestone for a coralline limestone 
exposed on Key Largo, Florida. This coralline limestone composes the upper Florida 
Keys from Soldier Key southward to Bahia Honda. The Key Largo Limestone also 
interfingers with the Miami Limestone along the western edge of Biscayne Bay. The 
formation is dominated by the corals Monsterea and Diplora, indicating an inner reef 
depositional environment (Hoffmeister 1975). The age of upper unit has been dated as 
approximately 125,000 years BP (Osmond et al. 1965; Broeker and Thurber 1965). 


Tertiary Deposits 
Tamiami Formation 


The Tamiami Formation, originally named by Mansfield (1930), 1s the principal 
Tertiary deposit in the LECsR model area. The formation covers most of south Florida 
from Lake Okeechobee southward and stretching across the entire State from the Gulf of 
Mexico to the Atlantic Ocean. It outcrops in the general vicinity of Naples and the Ten 
Thousand Islands and dips eastward and northward. The formation thickens to the east 
and the southeast. It inter-fingers with the Caloosahatchee Formation (Parker and Cooke, 
1944). The two primary members of the Formation are the Pinecrest Sand Member and 
the Ochopee Limestone Member. It is primarily composed of sandy limestone, 
calcareous sandstone, quartz sand, clay and marl, and often contains solution holes filled 
with soil (Parker and Cooke, 1944). The formation was deposited during the Pliocene to 
late Miocene epoch. The upper Tamiami Pliocene deposits 


Pliocene Deposits 


Pinecrest Sand Member 


The Pinecrest sand members originally named Pinecrest beds by Olsson(1964). 
Pinecrest sand is well sorted with abundant well preserved shells (Scott, 1992). The 
Pinecrest sand member has three lithofacies: 1) Quartz and sand, 2) rudstone and 
floatstone and 3) mudstone. The quartz and sand is the most characteristic of the 
Pinecrest sand member (Reese and Cunningham, 2000). The Pinecrest sand member is 
thickest in central and south-central Miami-Dade County. 
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Ochopee Limestone Member 


The Ochopee Limestone Member is made up of limestone composed of very fine 
to fine quartz grains (Reese and Cunningham, 2000). The Ochopee Limestone is well 
indurated, phosphatic and variably sandy (Scott, 1992) The Ochopee Limestone is 
thickest in southwest Palm Beach County and northwest Broward (Reese and 
Cunningham, 2000). In Miami-Dade County the Ochopee Limestone pinches out and 
transitions to the overlying Pinecrest Sand (Reese and Cunningham, 2000). 


Miocene Deposits 


Unnamed Formation 


The boundary between the Tamiami Formations Ochopee Limestone and the 
Peace River Formation is composed of a combination of quartz, gravel, sand, silt and 
carbonate rocks. At bottom of the formation is clean sand (Reese and Cunningham, 
2000). 


Peace River Formation 


The Hawthorn Group is composed of an upper and lower formation. The top of 
the Upper Hawthorn is the Peace River Formation. The Peace River Formation is 
principally carbonic rocks with siliciclastic sediment (Scott, 1990, Reese and 
Cunningham, 2000) interbedded with quartz, sand and clay. This represents the base of 
the SAS and concludes our description of the Hawthron Group. 


Chronostratigraphy 


Chronostratigraphy is the study of the ages of strata and may include the 
correlation of strata using relative or absolute methods of age determination. Large areas 
of south Florida are marine carbonate deposits. Emergence of these deposits by a 
reduction in sea level subjects the carbonate deposits to atmospheric conditions. The 
correlation of these marine units across the study area provides an indication of areas 
where potential reduced vertical movement of water may occur in the groundwater 
system. 


The development of advanced groundwater flow models to accurately simulate 
subsurface conditions in south Florida requires the vertical discretization of the Surficial 
Aquifer System (SAS) into an appropriate layering scheme. Historically, discretization of 
the SAS has generally fallen along hydrologic lines with minimal consideration for 
stratigraphic boundaries (Restrepo etal. 1992; Merritt 1998). Although previous 
approaches resulted in an adequate calibration for the model simulation periods, the layer 
schemes tended to cross distinct depositional boundaries and occasionally even 
formational boundaries. The approach utilized herein delineates the layering scheme 


Draft - 38 


LEC subRegional MODFLOW Model Documentation DRAFT 


along the stratigraphic boundaries originally proposed by Perkins (1977) while preserving 
the hydrologic properties of the aquifer system. 


The sediments deposited during the Quaternary Period are important to the 
interpretation of the stratigraphy of the SAS. The Quaternary sediments appear to have 
been deposited during the middle to late Quaternary Period. Within the Quaternary 
sediments, there are several episodes of subaerial exposure that resulted from the low sea 
level stands associated with Pleistocene glacial advances. Depositional indicators of these 
subaerial exposure periods include the development of dense caliche-type crusts with 
localized deposits of freshwater limestones on the underlying marine limestones. These 
caliche-type crusts generally have moderate to low hydraulic conductivities, and 
represent significant departures from the extremely permeable marine limestones 
generally located above and below these units. These subaerial exposure zones tend to 
retard vertical movement of water. 


Identification of these subaerial exposure surfaces becomes a critical component 
for the stratigraphic subdivision of the Quaternary sediments within the study area. 
Perkins (1977) indicates that these zones may only be deposited within a 6-inch to 1-foot 
layer. Therefore, geologic logs based on cuttings collected at 5-foot intervals or where the 
geologist was not specifically documenting subaerial exposure surfaces are not very 
useful for this type of analysis. Detailed continuous geologic logs are required to identify 
these subaerial exposure surfaces. 


In addition to the lithologic descriptions of the strata, age estimates are required to 
correlate the marine beds within a formation to world-wide eustatic sea level fluctuations. 
Parker and Cooke (1944) attempted to correlate a series of marine deposits located along 
the river bank of the Caloosahatchee River in south Florida to interglacial periods of the 
Pleistocene. Parker et al. (1955) refined the formational breaks in south Florida based, in 
part, on the earlier division proposed by Parker and Cooke (1944). Brooks (1968) 
subsequently refined the number of interglacial marine deposits of Pleistocene age. 
Building on the findings of these and other authors, Perkins (1977) differentiated the 
various Quaternary deposits of southeastern Florida into five distinct marine units 
punctuated by episodes of subaerial exposure (Q1, Q2, Q3, Q4, and Q5 in Figure 25). 
The estimated age of the uppermost marine unit was assigned based on radiometric 
dating of the Key Largo Limestone and Miami Limestone (Osmond 67 al. 1965; Broecker 
and Thurber 1965). Dates for the remaining four Quaternary units were estimated from 
amino acid racemization results from the Fort Thompson Formation (Mitterer 1975). 


Correlation of marine beds to oxygen isotope stages provides an understanding of 
the length of time the upper surfaces may have been exposed to terrestrial processes and 
also allows for an understanding of the facies distribution and paleotopography across the 
study area during a single marine transgression (Perkins 1977). 


Facies distribution and secondary solutioning also help control the hydraulic 


properties of the SAS. For example, buried in-situ reef systems may exist adjacent to 
sandy shallow marine flats. Although both facies may have been deposited during a 
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single high sea level stand, their hydraulic properties may be significantly different. 
Identification of facies distributions and general aquifer parameters associated with the 
individual facies need to be recognized prior to discretization of the model to allow for 
accurate simulation of the SAS. 


The tops and bottoms of the Holocene, the 5 Q Units, and the two Members of the 
Tamiami Formation were determined from a variety of sources including Caursas (1985), 
Caursas (1987), Reese and Cunningham (2002), Harvey et al (2000), and Giddings 
(1999). However, Perkins (1971) was used as the primary source. 


Figure 18 represents the location of five cross sections across the study area 
showing only the thickness of the Holocene sediments, Q Units 1 through 5, and the 
Pinecrest and Ochopee members of the Tamiami Formation. Only thickness of the units 
is shown to identify potential tops where sub-aerial exposure surfaces may occur. 
Figures 19, 20 and 21 are oriented in a north-south direction. No noticeable change 
occurs regarding the base of the SAS from north to south, although there appears to be a 
slight thinning northward. The Holocene sediments are extremely thin when compared to 
the Quaternary and Tertiary Sediments. The Tamiami Formation is the dominant 
formation along the west portions of the study area. The Quaternary sediments begin to 
thicken in Figure 20 and eventually are thickest along the coast (Figure 21). Deposition 
within the Quaternary sediments 1s relatively consistent although Q Units 2 and 3 appear 
to have had thicker depositional periods in many wells especially along the coast. 


Figures 22 and 23 are cross sections with an east-west orientation across the 
study area. These two Figures illustrate the drastic increase in thickness of the 
Quaternary sediments, and the thinning of the Tamiami Formation from east to west 
across the study area. The Quaternary Period appears to include the additional eastward 
expansion of the southeastern Florida peninsula. 
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Figure 17. Cross-section Base Map. 
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Figure 18. Cross-section A-A’. 
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Figure 19. Cross-section B-B’. 
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Figure 20. Cross-section C-C’. 
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Figure 21. Cross-section D-D'. 
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Figure 22. Cross-section E-E’. 


Groundwater Hydrology of the Surficial Aquifer System 


The Surficial Aquifer System (SAS) 1s an unconfined to semi-confined aquifer 
system, and includes all saturated sediments from the ground surface down to the clay 
and marls of the impermeable Hawthorn Group. The system is dominated by carbonate 
and clastic sediments ranging in age from the Holocene through the Miocene. The 
limestone beds constitute the major water-producing component of two aquifers: the 
Biscayne aquifer and the Gray Limestone aquifer. These two aquifers can grade into one 
another, especially along the coast. The SAS includes a composite of the Biscayne 
aquifer, the Gray Limestone aquifer along the coast, and the Water Table aquifer. 


The Biscayne aquifer generally comprises all Quaternary sediments including, the 
Fort Thompson Formation, Anastasia Formation and Miami Limestone (Fish 1988). The 
Gray Limestone aquifer is a permeable zone within the Tamiami Formation of Pliocene 
to Miocene age (Fish and Stewart 1991). In western Broward and Miami Dade Counties, 
the Biscayne aquifer is separated from the lower Gray Limestone aquifer by a 
semiconfining unit of the upper Tamiami Formation, with the two aquifers merging into a 
single hydrologic unit along the coast. The remaining portion of the SAS, the Water 
Table aquifer, has much lower yields than the Biscayne and Gray Limestone aquifers, 
and is located in western and northern Palm Beach County and all of Martin County. 
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Biscayne Aquifer 


The principal groundwater resource for the southeast coast of Florida is the 
Biscayne aquifer. This groundwater system provides most of the fresh water for public 
water supply and irrigation demands within the model area, and is also used to maintain 
canal and lake systems as well as provide groundwater base flows to the Everglades and 
important estuaries including Biscayne Bay and Florida Bay. The Biscayne aquifer is an 
unconfined aquifer and a component of the SAS underlying most of southeast Florida as 
shown in Figures 3 and 24. 
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Figure 23. Geologic Cross-Section of the Surficial Aquifer System in Southeast Florida. 


The major geologic deposits that comprise the Biscayne aquifer include the 
Miami Limestone, Fort Thompson Formation, Key Largo Formation and the Anastasia 
Formation. Information about these deposits is summarized in Figure 25. The base of the 
Biscayne aquifer is generally the contact between the Quaternary Fort Thompson 
Formation and the Pliocene-Miocene Tamiami Formation. 


The Biscayne aquifer is composed of interbedded, unconsolidated sands and shell 
units along (in combination) with varying thicknesses of consolidated, highly solutioned 
limestones and sandstones (Shine ef al. 1989). In general, the Biscayne aquifer contains a 
smaller percentage of sand and a greater percentage of solutioned limestone than the 
other portions of the SAS. The Biscayne aquifer is one of the most permeable aquifers in 
the world, with estimated transmissivities in excess of 7 million gallons per day, per foot 
of drawdown (Parker et al. 1955). 
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Figure 24. Lithostratigraphic and Geohydrologic Units of the Surficial Aquifer System in 
Southeast Florida. Source: Adapted from Reese and Cunningham 2000; Perkins 
1977. 






Gray Limestone Aquifer 


Underlying the Biscayne aquifer and composing the base of the SAS is the 
Tamiami Formation. This formation consists of an upper clastic unit with low to 
moderate permeability. Beneath the upper clastic unit is an intermediate limestone unit 
locally referred to as the Gray Limestone aquifer or an eastern extension of the Lower 
Tamiami aquifer which is characterized by moderate to high permeability. The Gray 
Limestone aquifer is located in the western portion of the study area where it 1s semi- 
confined to unconfined. In the eastern portion of the study area it interfingers with the 
Biscayne and Water Table aquifers to form the SAS. The base of the SAS is also 
generally considered to be the base of the Tamiami Formation, which consists of a lower 
clastic unit of moderate to low permeability (Fish and Stewart 1991). 


The Gray Limestone aquifer includes the Ochopee Limestone Member of the 
Tamiami Formation. It is overlain and underlain by upper and lower confining to 
semiconfining units throughout most of the Lower East Coast of Florida. These confining 
units are usually composed of siliclastics of low to very low hydraulic conductivity 
(Reese and Cunningham 2000). 


Hydrogeologic Parameters 


Aquifers transmit water at a rate based upon hydraulic conductivity, storage 
properties, and the hydraulic gradient. Henry Darcy first introduced an equation in 1856 
to express factors controlling the movement of groundwater. This equation, known as 
Darcy’s Law, is expressed as follows: 
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Oh 
= Ka e 
e ol 


Where Q is the quantity of water, K is the hydraulic conductivity, A is the cross- 
sectional area at a perpendicular to the flow direction, and the hydraulic gradient is the 
change in head per change in length, dh/dl. 


The primary aquifer characteristics required for simulating flow in an unconfined 
aquifer are transmissivity and specific yield. Transmissivity is the capacity of an aquifer 
to transmit water and is determined by multiplying the hydraulic conductivity of the 
aquifer by the thickness of the aquifer. Darcy’s Law can be rewritten in terms of 
transmissivity by expressing the cross-sectional area component of the equation as the 
product of the width, w and the aquifer thickness, b. The revised equation is as follows: 


ah 
siy 
L-Iw 


Therefore, in order to solve the groundwater flow equation, the hydraulic 
conductivity, the aquifer thickness or the transmissivity of the aquifer will need to be 
identified. These parameters are generally determined from doing some form of aquifer 
performance test (APT) on wells with known aquifer and screen interval depths and 
pumping rates. Laboratory analysis can also help determine some estimated properties 
of the aquifer as well. In the study area, there are three types of data that will be used to 
determine the hydraulic conductivity and aquifer thicknesses. They are geologic control 
wells with continuous cores, APT”s, and specific capacity tests on existing production 
wells. Locations of the three data types used to estimate hydraulic conductivity are 
shown in Figure 26. The hydraulic conductivities derived from the Holocene, five 
Quaternary, and Tertiary units are depicted in Figures 27 to 36. 
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Figure 25. Locations of Control Points Used to Estimate Hydraulic Conductivity 
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The lack of APTs in Miami-Dade and Broward Counties in Figure 28 is 
misleading. Fish and Stewart (1991) and Fish (1988) did a very detailed analysis of the 
hydraulic conductivity of the geologic control wells in those two counties. So the 
hydraulic conductivity for each section of each geologic control well has inherently 
incorporated into it the large number of APT’s that have been conducted in south Florida 
over the years. Fish (1988) and Fish and Stewart (1991) used Causaras’ (1985, 1987) 
detailed descriptions of continuous cores drilled in Miami-Dade and Broward Counties. 
Through a complex process, they determined hydraulic conductivity ranges continuously 
down each geologic control well. The information utilized in developing this data 
included geologic descriptions, lithology, previous APT’s, new APT’s, flow rates while 
drilling the test wells, hydrologic inferences from samples, sieve analysis, and the degree 
of sorting and grain size for clastic sediments (Fish and Stewart 1991). A similar process 
was utilized for the geologic control wells in Palm Beach and Martin Counties. 


The geologic control wells allowed for a continuous profile of the hydraulic 
conductivity from ground surface to the base of the SAS. By overlaying our thickness 
maps for each of the five Q units and the members of the Tamiami Formation discussed 
in the Chronostratigraphy Section, estimated hydraulic conductivities for each unit were 
developed; they were based upon the weighted average for the hydraulic conductivities 
associated with that unit. This approach was used for every geologic control well in the 
study region. 


Fish (1988) and Fish and Stewart (1991) also estimated the transmissivity of the 
production zones from specific capacity tests for a large number of municipal public 
water supply production wells in Miami-Dade and Broward Counties (Figure 27). 
Specific capacity tests are general short term pumping tests where only water levels in the 
production well are recorded for both drawdown and recovery. Utilizing the following 
empirical equation, they were able to estimate the transmissivity of the production zone 
for each well. 


i= 270€ 
S 

Where T is the transmissivity, Q is the well discharge in gallons/minute, and s is 
the drawdown, in feet. 


Hydraulic conductivities were determined from the screened or open-hole interval 
for each production well and by the unit(s) that it intercepted as determined from the 
chronostratigraphic maps. Since the production wells generally are screened across the 
main production zones of the aquifer, and the specific capacity tests were of short 
duration, it was assumed that the screened or open-hole interval was the main source of 
water and the transmissivity represented that zone of the aquifer. The hydraulic 
conductivity of the Chronostratigraphic unit(s) was then determined by dividing the 
transmissivity by the screened interval of the production well and a weighted average was 
then given to that unit(s). 
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APTs were used to help develop the geologic control wells in Palm Beach and 
Martin counties. In addition, they were used to provide point locations for the composite 
transmissivity of the aquifer as well as an estimate of the hydraulic conductivity of the 
production zone. Typical methods were used to analyze the results of the APT including 
Neuman, Hantush-Jacob, Boulton, and Cooper. In addition to providing transmissivity 
information, the APT’s allowed for an estimation of the storativity of the aquifer. 
Storativity values general ranged between 2x10-1 to 5x10-5 with an average of 
approximately 1x10-3. 


For the unconfined section of the aquifer, the specific yield needs to be known to 
ultimately be able to determine transient flow through the system. Merritt (1996) utilized 
rainfall and the corresponding rise in water levels in monitoring wells to estimate the 
porosity of the Miami Limestone. Short-term-duration rain events and predevelopment 
data were utilized to minimize the impact of runoff and rapid removal of water from 
secondary canal systems. It was determined that a specific yield of between 0.20 and 
0.25 was reasonable for the Miami Limestone. Considering the similarities of the Miami 
Limestone with the other limestone units in the area, a specific yield of 0.2 for the entire 
SAS appears reasonable. 


The vertical hydraulic conductivity is also a necessary component when 
evaluating flow through an aquifer system. The vertical hydraulic conductivity allows 
for movement vertically through the stratigraphic layers and is generally a fraction of the 
horizontal conductivity which is the preferred flow path. The vertical anisotropy ratio is 
the ratio of vertical to horizontal conductivity. Although numerous authors have 
estimated this ratio, the works of Geotrans (1986) and Langevin (2002) probably provide 
the best estimate of this ratio. Both of these works address movement of the salt water 
interface in south Florida in a three dimensional manner. The salt water interface in the 
study area is evaluated in a cross sectional manner and takes on an “S” curvature. The 
“S” curvature is where the toe of the saltwater interface is further inland than at the top of 
the aquifer. Geotrans (1986) reported that the degree of this “S” curvature is sensitive to 
the anisotropic ratio and suggested utilizing a ratio of 10:1. Langevin (2002) utilized a 
value of 100:1, but through sensitivity analysis reported no significant effect when this 
number is reduced. 
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Figure 26. Horizontal Hydraulic Conductivity (ft/day) from Specific Capacity Tests. 
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Figure 27. Horizontal Hydraulic Conductivity (ft/day) from Aquifer Performance Tests. 
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Figure 28. Horizontal Hydraulic Conductivity (ft/day) of the Holocene Sediments. 
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Figure 29. Horizontal Hydraulic Conductivity (ft/day) of the Q5 Unit. 
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Figure 30. Horizontal Hydraulic Conductivity (ft/day) of Q4 Unit. 
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Figure 31. Horizontal Hydraulic Conductivity (ft/day) of Q3 Unit. 
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Figure 32. Horizontal Hydraulic Conductivity (ft/day) of Q2 Unit. 
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Figure 33. Horizontal Hydraulic Conductivity (ft/day) of Q1 Unit. 
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Figure 34. Horizontal Hydraulic Conductivity (ft/day) of the T2 Unit. 
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Figure 35. Horizontal Hydraulic Conductivity (ft/day) of the T1 Unit. 
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Saltwater Intrusion 


A primary issue of concern for the Biscayne aquifer is contamination by a 
westward-moving saltwater front (SFWMD 2000). Several factors have contributed to 
the landward migration of the interface. Draining the Everglades and increasing wellfield 
pumping lowered the water table elevations and also increased peat and muck soil loss 
due to oxidation and compaction (Renken et al. 2005). As a result of extensive previous 
studies, the position of the current saltwater interface is a result of the following 
mechanisms: lateral landward movement when seawater moves from the Atlantic Ocean 
into the SAS, seepage from tidal canals containing saline water, and upconing when 
wellfield pumping causes relict seawater to move upward (Renken er al. 2005). 


The 250 milligrams per liter (mg/L) saltwater intrusion line (Figure 20) at the 
base of the SAS is used in south Florida to evaluate potential saltwater contamination 
when wellfield withdrawals are increased. Additionally, the SFWMD has implemented 
minimum levels at several surface water structures to prevent landward migration of the 
saltwater front (SFWMD 2000). 


One of the primary concerns associated with the Biscayne aquifer along the 
southeast coast of Florida is the threat of salt water intrusion. In order to minimize the 
threat of the inland migration of the saline interface, a sufficient fresh water head should 
be maintained landward of the saline interface. If a sufficient fresh water head is not 
maintained, then salt water intrusion can occur. Such was the case in 1939 when more 
than 10,000 water supply wells went salty in south Florida, including the partial loss of 
five major wellfields (Parker et al. 1955). 


It is generally believed that one of the major reasons that salt water 
intrusion has historically occurred in south Florida is the loss of the fresh water mound 
behind the coastal ridge system (Parker et al. 1955; Fish and Stewart 1991). Prior to 
development, a large fresh water mound would develop behind the Atlantic Coastal ridge 
during the rainy season. Ground water flows would be so great seaward that springs 
would boil up off the coast providing freshwater to early mariners. As south Florida 
became developed, a series of canals and drainage ditches were constructed which 
drained this mound. The result was a significant decline in groundwater flow towards the 
ocean and, consequently, the inland migration of the saline interface. In addition to the 
canal network, large wellfields located adjacent to the coast have also been responsible 
for localized salt water intrusion problems. 


In 1981, a severe drought in south Florida resulted in widespread movement of 
the saline interface. By 1987, the City of Hallandale had to permanently reduce pumpage 
by 50 percent and shut down their primary wellfield due to salt water intrusion (SFWMD 
1987). Over the next couple of years, several golf courses in the Hollywood/Hallandale 
area were also required to curtail or eliminate groundwater withdrawals due to salt water 
intrusion. Koszalka (1994) reported that the saline interface continued to move inland in 
Broward County between 1980 and 1990 due to the lowering of regional groundwater 
levels and wellfield pumpage. 
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In the late 1800's, two scientists working independently observed that salt water 
was not found at sea level in an aquifer but at a specified depth depending upon the 
height of the fresh water in an aquifer. The depth to the salt water was approximately 40 
times the height of the fresh water above sea level. They concluded that the distribution 
was attributed to a hydrostatic equilibrium existing between the two fluids of different 
densities (Todd, 1980). In general, they observed that for each 40 feet of aquifer 
thickness, a one foot fresh water head will need to be maintained to stabilize the salt 
water front under static conditions. The equation derived to explain their observations is 
referred to as the Ghyben-Herzberg relationship. The aquifer system along the southeast 
coast of Florida ranges in thickness from approximately 100 to 300 feet. Therefore, in 
order to ensure that salt water intrusion does not occur within these aquifer systems, 
freshwater heads of between 2.5 and 7.5 feet will need to be maintained based upon the 
Ghyben-Herzberg relationship. 


Work done at Cutler Ridge in south Dade County indicates that the salt water 
front is not static as assumed by the Ghyben-Herzberg relationship (Kohout 1960). In 
addition, the observed position of the saline interface is several miles seaward of the 
position calculated using the Ghyben-Herzberg relationship. Kohout (1960) observed 
that as saltwater moved inland, a significant portion of the diluted sea water was 
circulated back toward the sea along the zone of diffusion. It is estimated that up to 20 
percent of the salt water intruded into the aquifer is returned seawater with the remaining 
80 percent being retained in the aquifer (Kohout 1960). This cyclic flow acts, in part, as 
a deterrent to salt water intrusion due to a percentage of the salt water being returned to 
the sea. It also indicates that salt water intrusion is dynamic and not static as assumed by 
the Ghyben-Herzberg relationship. However, Kohout’s (1960) results suggest that the 
use of Ghyben-Herzberg levels to control the position of the salt water interface may be 
adequate because it tends to over predict the levels needed to maintain the position of the 
saline interface. 


The City of Hallandale in southeastern Broward County, Florida has been an area 
susceptible to salt water intrusion for a number of years. A series of monitor wells 
located perpendicular to the coast have recorded the inland migration of the saline 
interface for over 35 years. Evaluation of the data suggests the saltwater front has 
consistently migrated inland at a rate of approximately 80 feet per year. Anderson et al. 
(1988) conducted a detailed evaluation of the salt water interface in the vicinity of 
Hallandale utilizing a coupled flow/solute transport three dimensional finite element 
model. They evaluated several potential causes for the continued intrusion including 
wellfield pumpage, rainfall deficiency and lowering of inland canal stages due to 
urbanization. Although their model could not localize the problem causing the salt water 
intrusion, the results clearly demonstrate the sensitivity of the saline interface to 
maintained ground water stages. Lowering the inland canal stages by even several tenths 
of a foot could result in widespread movement of the saline interface. In addition, 
Anderson et al. (1989), observed from their model simulations that a significant lag time 
existing between the lowering of hydraulic heads and the subsequent movement of the 
saline interface. 
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Because of these issues concerning movement of the saline interface, the 
SFWMD developed a series of Minimum Levels for the primary coastal canal structures 
based upon both analytical and numerical analysis (SFWMD 2000). Based upon the 
result of the analysis, it was recommended that canal stages at key coastal primary 
discharge structures be maintained a minimum level for a specific period of time during 
the year to begin to stabilize existing or future movement of the saline interface. Table 3 
is the recommended minimum canal operational levels for the Biscayne aquifer. 


Table 3. Minimum Canal Operational Levels (ft NGVD) for the Biscayne Aquifer. 


Wet 


Season Average ET 
Canal/Structure C Management 
ontrol Canal Level C 
ontrol Level 
Level 
C-51/S-155 8.50 8.12 7.80 
C-16/S-41 8.20 8.23 7.80 
C-15/S-40 8.20 8.39 7.80 
Hillsboro/G-56 7.70 7.43 6.75 
C-14/S-37 7.20 6.82 6.50 
C-13/S-36 5.60 4.43 4.00? 
NNR/G-54 4.00 3.68 3.50 
C-9/S-29 3.00 2.16 1.80 
C-6/S-26 4.40 2.55 2.50° 
C-4/S-25B 4.40 2.55 2.50° 
C-2/S-22 3.50 2.86 2.50° 


a. Duration Criterion — water levels within the above canals may fall below 
the proposed minimum canal level for a period of no more than 180 days per 
year. 


b. These levels will be maintained if sufficient water is available. 


c. Minimum Canal Operational Levels needed to protect against MFL 
violations during drought conditions. 
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Figure 36. Location of the Saltwater Interface. Source: After Renken et al. 2005 
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WATER USE 


Within the study area, ground-water and surface-water withdrawals for public 
water supply, industrial (including power plants), landscaping and agricultural uses were 
estimated from the South Florida Water Management District permit data base. In 
Miami-Dade, Broward, Palm Beach and Martin counties there are approximately 7,000 
water use permits issued with a total daily permitted demand of approximately 3.2 billion 
gallons per day in the year 2004. The main agricultural users of surface water are the 
citrus groves in the northern portion of the study area and the Everglades Agricultural 
Area. The public water supply systems in the study primarily use the Surficial Aquifer 
System with the exception of some Floridan aquifer use in the northern portion of the 
study area. Surface water as a public water supply source is restricted to the small 
communities surrounding Lake Okeechobee and the City of West Palm Beach. 


Of the top five largest permitted users of water in the study area, three are for 
Public Water Supply, including Miami-Dade Water and Sewer Department and Palm 
Beach County Water Utilities Department, while the other two are for diversion of 
surface water from the regional system into secondary canal systems and include the 
Lake Worth Drainage District located in central and southern Palm Beach County and 
Old Plantation Water Control District in central Broward County. These diversion 
permits accomplish multiple uses including re-hydration of urban wetland systems, 
provide water for agricultural and landscaping purposes, and maintain the vast system of 
secondary canals to recharge local wellfields and help prevent saltwater intrusion. 


The two largest use classes of water in the region are agriculture and public water 
supply, which account for approximately two-thirds of the entire demand. The main 
industrial use is primarily associated with limestone rock mining. Palm Beach County 
has the greatest amount of permitted water use in the region with over 1 billion gallons 
per day permitted. Domestic self-supplied use is not included in the analysis because, 
with the exception of some rural areas, the majority of the population is supplied from a 
public water supply system. The main domestic use is through irrigation of home-site 
landscaping with private wells, which even then 18 prohibited in some communities. 


Table 4. Permitted Water Use in South Florida, Expressed in Million Gallons Per 


Day (MGD). 

Mi 

am 
Broward i Palm Beach Martin Totals 

Da 

de 
Agriculture 23.8 761.8 211.8 1034.2 
Dewatering 76.4 274.5 18.4 426.5 
Diversion 96.2 0.0 189.6 0.0 285.8 
Golf Course 18.6 Pc 48.1 8.7 85.7 
Industrial a5 ES 32.2 8.0 144.1 
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Landscaping 18. 
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Figure 37. Water Use by Class in 1994. 
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Figure 38. Water Use by County in 1994. 
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Figure 39. Public Water Supply Wells in Southeast Florida. 
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Figure 40. Consumptive Use Locations Other than Public Water Supply Categorized by Type of 
Use in Southeast Florida. 
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CHAPTER 3 


Simulation of the Flow System 


The process of model design for the LECsR truly began with the first subregional 
groundwater models at the SFWMD in the early 1990’s. As stated in Chapter 1, the 
knowledge base in south Florida is growing and there are many research efforts 
contributing to this wealth of knowledge. Development of LECsR relied on the 
knowledge gained from the previous subregional models and experienced staff 
throughout the SFWMD, as well as the most current, best available data. 


Model design for the LECsR Model has been an iterative process. The first 
version of this model was initiated three years ago and resulted in a calibrated model 
which was further refined due to a request by management to extend the calibration 
period by ten years. During this time, the model evolved by incorporating new data and 
modeling capabilities within MODFLOW-96. These events required revisiting the 
modeling protocol. It is through this iterative process of re-formulating the conceptual 
model, re-designing and re-calibrating that has produced the current version of the 
LECsR Model. In the future, the model will evolve and mature using this iterative 
process as needed to keep up to date with the best available data without invalidating the 
peer review (i.e., a step-wise post-audit). 


COMPUTER CODE SELECTION 


Once modeling objectives have been established and the predominant hydrologic 
processes within the area of interest have been determined, a model code that can meet 
the model development and application objectives is selected. MODFLOW-96, a code 
created by the U.S. Geological Survey (USGS) (Harbaugh and McDonald 1996; 
McDonald and Harbaugh 1988), was selected for this purpose for the following primary 
reasons: 


e It has been widely accepted in the groundwater modeling profession 
for over 20 years. 


° The code is well documented and within the public domain. 


° The code is readily adaptable to a variety of groundwater flow 
systems. 


e The code is modular and easily facilitates any modifications required 
to enable its application to the types of unique groundwater flow 
problems encountered in south Florida. 


e MODFLOW was used to develop existing groundwater flow models 
located within the LEC that could be upgraded to meet the current 
objectives. 
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e [n the future, MODFLOW models may be updated to the variable 
density (with solute transport) transient flow model, SEAWAT 
(Langevin et al. 2003), which is fully compatible with MODFLOW. 


MODFLOW with District Source Code 


MODFLOW simulates groundwater flow in aquifer systems using the finite- 
difference method. The aquifer system is divided into rectangular or quasi-rectangular 
blocks by a grid (Figure 41). The grid of blocks is organized by rows, columns, and 
layers, and each block is commonly called a cell. 


ie 000 


X III DASS 





XR ef, > ALLAN 





Cells inside of 
aquifer system 


Figure 41. Example of Model Grid for Simulating Three-Dimensional Groundwater Flow. 


For each cell within the aquifer system, the user must specify aquifer properties. 
Also, the user specifies information relating to wells, canals, and other hydrologic 
features for the cells corresponding to the locations of the features. For example, if the 
interaction between a canal and an aquifer system is simulated, then for each cell 
traversed by the canal, the required input information includes layer, row, and column 
indices; canal stage; and hydraulic properties of the channel bed. Also, MODFLOW 
allows the user to specify which cells within the grid are part of the groundwater flow 
system and which cells are inactive (i.e., outside of the groundwater flow system). 


The MODFLOW model code consists of a main program and a series of 
independent subroutines called modules. The modules, in turn, have been grouped into 
packages which each deal with a particular hydrologic process or solution algorithm. The 
packages used for the LECsR Model development and application, including those 
developed or enhanced by SFWMD staff and contractors, are shown in Table 5. The 
enhanced packages have been added to the MODFLOW-96 computer source code over 
the last decade. Therefore, the SFWMD has not migrated to MODFLOW-00, but intends 
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to do so in the near future. The planned migration will include updating MODFLOW- 
2000/SEAWAT-2000 with the add-on (or enhanced) packages. 


Table 5. MODFLOW Packages Used in the LECsR Model. 


Basic and Output Defines stress periods, time Handles the primary 
Control steps, starting heads, grid administrative tasks 
(BAS and OC) specifications, units, and output associated with a 
specifications. simulation. 


Block-Centered Specifies steady state vs. Derived primarily from 
Flow (BCF) transient flag, cell sizes, geologic data used to 
anisotropy, layer types, and construct the model. 
hydrogeologic data for each 
layer. 


Utility Generation Creates input files during model Generates input for Well, 
(UGEN) execution by linking static (time General Head Boundary, 
invariant) data with time series Diversion, River, 
(variant) data. Reinjection Drainflow, and 
Drain Packages. 


Surface Water Stresses and Processes 


Recharge (RCH) Simulates aerially distributed Preprocessed using an 
recharge to a water table during Agricultural Field-Scale 
each stress period. Irrigation Requirements 
Simulation (AFSIRS) 
based ET- Recharge 
model. 


Preprocessed using an 
AFSIRS based ET- 
Recharge model; 
Saturated ET rate 
diminishes with increasing 
water table depth. 
Canal stages are usually 


based on measured stages 
or control elevations. 














Simulates removal of water from 
the water table via transpiration 
and direct evaporation. 






Evapotranspiration 























(EVT) 
River (RIV) Simulates groundwater 
interchanges with canals that 
can either recharge or drain the 
aquifer. 
Drain (DRN) Essentially the same as the 
River package except that 
canals can only drain the aquifer 
and water removed by the drains 
is removed permanently from the 
model. 


Simulates the effects of water 
control structures (either 
pumping stations or gravity flow 
drains) on water levels. 









Canal stages are usually 
based on weir elevations 
and topography. 















Allows simulation of 
operational rules at 
structures and diversion of 
water, including external 
sources like locks and 
external sinks like the 
ocean. 


Allows simulation of 
seepage control systems 
where water is diverted out 
of the seepage canal 
based on a control level. 







Diversion (DIV) 









Reinjection Essentially the same as the 
Drainflow (RDF) Drain package except that it 
allows water to be redirected to 


another location in the model 
instead of being permanently 
removed from the model. 


Horizontal Flow Simulates thin vertical barriers to Used for siurry walls. 
Barrier (HFB) horizontal flow. 
Wetland (WTL) Simulates the overland flow in Enhanced to also simulate 
wetlands using the uppermost barriers to flow like levees. 
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model layer. OF 


General Head Simulates groundwater Boundarv stages are 
Boundarv (GHB) exchange between selected usually based on 
cells and a specified boundary measured stages. 
as a function of water level 
difference. 


Water Supply and Management 


Well/Multiple Wells Simulates withdrawals from Includes Public Water 
(WEL) wells. Enhanced by the District Supply (PWS), irrigation, 
to read multiple input files and Aquifer Storage and 

Recovery (ASR) wells. 


Trigger (TRG) Simulates wellfield withdrawal Cutback zones are based 
cutbacks as a function of water on calibration of model to 
level in trigger cells and in Lake historical water shortages 

Okeechobee; simulates LEC including the historical 
water shortage policy associated Lake Okeechobee 
with saltwater intrusion cutbacks. 


Solution Algorithms 


Strongly Implicit A mathematical solution Enhanced by District to 
Procedure (SIP) algorithm internal to the model improve model stability by 
allowing maximum 
iterations to be exceeded 
and the closure criteria to 
be loosened for a short 
number of stress periods. 


Preconditioned A mathematical solution PCG2 uses the 
Conjugate Gradient algorithm internal to the model preconditioned conjugate- 
(PCG) gradient method to solve 
the equations produced by 
the model for hydraulic 
head. 


Output Management 


Output Control Sums cell-by-cell flows to reduce Set up in the output 
Summation output size control. 
(OC) 


Multibud (BUD) Outputs an internal water budget Frequency of flow output is 
for a set of specified cells at a user-specified (e.g., daily 
given frequency to reduce output or monthly). 
size. 





Add-On District Packages 


The modular structure of MODFLOW readily allows for modifications or the 
creation of additional packages. The SFWMD has taken advantage of this feature and 
developed several additional packages and made enhancements to the original 
MODFLOW-96 code which are described in this section. Additional information about 
these packages may be obtained by contacting the SFWMD. 


Wetland Package 


The Wetland (WTL) Package was developed by the SFWMD and the Center for 
Hydrology and Water Resources at Florida Atlantic University (FAU) (Restrepo et al. 
1998). The current version has been revised and improved by FAU in association with 
the SFWMD, as the package has been applied to subregional models in south Florida, to 
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fulfill the need for a sound, physical-based representation of wetlands and surface water- 
groundwater interaction. 


Prior to the development of the Wetland Package used in this model, wetlands 
were often simulated as constant head or general head boundaries in groundwater models 
(Merritt 1995). These approaches have limitations, especially with the periodic drying 
and rewetting of wetland cells that occur when water levels fluctuate below and above 
land surface (Wilsnack et al. 2001). 


The Wetland Package is incorporated into the MODFLOW groundwater model 
code (Harbaugh and McDonald 1996) and enables the top layer of the grid system to 
contain overland or groundwater flow. The wetland module can account for vegetation 
characteristics, simulation of sheet flow, sloughs, levees, and barriers, and 
evapotranspiration as shown in Figure 42. This figure also illustrates the interaction 
between the channel flow (e.g., slough) and aquifer and between sheet flow and the 
aquifer. In the model, the top layer can include sheet flow together with aquifer and the 
muck/peat layer. This add-on package makes it possible to simulate the areal expansion 
and contraction of wetland systems and the associated water routing (horizontal and 
vertical) in response to different hydrologic conditions (Restrepo ef al. 1998). 
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Figure 42. Schematic of Wetland System Representation. Source: After Restrepo ef al. 1998. 


The Wetland Package uses a semi-empirical Manning-type equation (Kadlec 
1990) rather than Darcy's Law to represent surface water movement through dense 
vegetation. The Kadlec equation is given by: 


= K hr S 
d ii f Equation 1 
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where q is the discharge per unit width (LT); Kw is the hydraulic conductance 
coefficient for overland flow (L//T/L'); h is the flow depth (L); 6 is an exponent related 
to microtopography and the stem density-depth distribution; Sr is the hydraulic gradient; 
@ is an exponent that reflects the degree of laminar or turbulent flow conditions. An 


explanation of the mathematical formulations for the package may be found in Restrepo 
et al. (1998). 


This package is suitable for modeling the wetting and drying of the wetlands by 
combining vegetation and soils (under the ponded water in Figure 42) as part of model 
layer 1. This approach allows the flow equations to remain valid when the water surface 
falls below the soil surface (1.e., Darcy flow) as shown by the level 2 line in Figure 43. 
When the model is applied to a water ponding case (level 1), the effect of the value of 
soil transmissivity is insignificant relative to the corresponding value for the surface 
water, which is governed by Kadlec. However, the soil layer plays an important role in 
the vertical direction, as is expected. 
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Figure 43. Surface and Groundwater Representation in Weiland Package. Source: After 
Restrepo ef al. 1998. 


The Wetland Package simulates preferential flow paths and vertical impermeable 
barriers to both horizontal overland and groundwater flow. These features will either 
increase or decrease the flow rate. Levees and sloughs are indirectly defined as a 
hydraulic characteristic by using the cell faces for the interblock transmissivity to 
designate preferential flow orientations. In the case of fully or partially penetrating 
barriers, the package has the option of completely restricting flow through the cell face or 
allowing some seepage through the face. The establishment of barriers to flow requires 
the user to specify which cell faces in the modeled feature are perpendicular to the cell 
flow in order to block the flow. Specifying preferential flow paths requires the user to 
select the cell faces in the modeled feature that will follow the primary spatial extent of 
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the slough to allow water to move freely down a path. This methodology works best 
when the principal direction corresponds to the grid orientation. 


Diversion Package 


The Diversion (DIV) Package was developed by the SFWMD and the Center for 
Hydrology and Water Resources at Florida Atlantic University (FAU) (Restrepo et al. 
1998). The current version has been revised and improved by FAU in association with 
the District, as the package has been applied to subregional models in south Florida. The 
package was created to simulate the effects of man-made structures such as pumping 
stations or drains on groundwater and surface water levels. Originally designed as a part 
of the Wetlands Package (Restrepo et al. 1998), it was subsequently developed as an 
independent package in order to utilize its features in a wider range of situations. 





The Diversion Package defines water control structures as those structures that 
regulate groundwater or surface water levels on a seasonal basis. The operational rules of 
these structures will vary seasonally and in response to extreme weather events such as 
droughts or hurricanes. In general, during wet periods, water control structures remove 
water from a basin to prevent flooding. In dry periods, water is retained in basins in order 
to maintain groundwater levels for crops or to deter saltwater intrusion in coastal areas. 


The Diversion Package routes water according to operational rules or physical 
constraints of the system by conserving mass balance (keeping the water in the system). 
The water is supplied by a group of source cells, and is distributed to a group of sink cells 
associated with each structure. Sources and sinks can be external or internal cells or 
areas outside of the model domain. External cells are located outside of the active model 
domain; internal cells are always located inside the active model area. Sources may be 
located outside the active domain (e.g. external surface storage which makes deliveries 
within the active area). Sinks may be located inside the active domain (e.g., a swale 
which receives water from a source). Source and sink sets should be defined based on 
particle travel times of one to two times the length of the time discretization. 


The Diversion Package implemented flows due to pumping by adding the source 
and destination (layer, row, and column) of the cell to which the drain flow is directed. If 
no destination cell is specified or if the destination cell is inactive, the package functions 
exactly like the Drain Package. In the case of pumping, the withdrawal rate is calculated 
as the product of the pumping flow (given by the user) and efficiency factor. This 
package can be used, for example, to represent the removal of water from a reservoir and 
dispersal onto a nearby wetland or swale; or to an injection into an ASR well. 


Structures have an optimum level at the head (upstream) and at the tail 
(downstream). The package redistributes the ground water flows from source to sink 
cells, according to the operational rules of the structure which are based on control levels. 
The case of diversion flow by a pumping station is depicted in Figure 44 below. 
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Cross-Section Case of Diversion Flow by Pumping 
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Figure 44. Cross-Section Diagram of Diversion Flow at a Pumping Station. Source: After 
SFWMD 2005. 


The flow rate variable, q. is supplied by the user for each stress period and 
represents the movement of water from source to sink (for qç > 0) or from sink to source 
(q. < 0). The term on level source refers to the lowest water level limit in a source at 
which water can be released to the sink. Off level sink is the highest water level limit in 
the sink cells at which they can receive water from the source. In the case of positive flow 
(q. > 0), the water level in the source cell (hydraulic head at source = hsource) will be 
greater than or equal to the on level source limit and the water level in the sink (hydraulic 
head at sink = hsinx) will be less than or equal to the off level sink limit. Off level source 
1s the highest water level at which source cells at the structure can receive water from the 
sink. On level sink is the lowest level at which water in sink cells can release water 
through a structure to source cells. When flow is negative (qe < 0), hsource Will be less than 
or equal to the off level source limit, and baak will be greater than or equal to the on level 
sink limit. The on/off source/sink levels function as operational rules for the structure, 
and are specified by the user for both wet and dry seasons. 


The package is also used in a wider range of situations such as modeling 
reservoirs and impoundments, spreader swales and flow ways. Current enhancements 
allow routing water as series of “cascading” sub-basins and representation of several 
rules at same location. Flows can be routed from upstream to downstream basins. The 
basins can include natural and urban areas, but must be specified as active Wetland 
Package cells. The Wetland Package, in combination with the Diversion Package, will 
facilitate both channel or sheet flow to move the water to the downstream outlet. In this 
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case, basin-scale rules of operation or natural drainage characteristics (e.g., change in 
topography within the source) can be used to control the amount of flow leaving the 
basin. 





Flow (q. > 0) can be predetermined from estimated or historical basin runoff. The 
model estimates the quantity of runoff that will be removed from the sub-basin and 
subsequently routed downstream according to the drain elevation. This quantity may be 
reduced due to the number of cells that actually meet the criteria which are defined by the 
rules of operation. When a structure does not exist (e.g., tributary to a river), the drain 
elevation is derived from the topography. The discharge rate released by this structure is 
constrained by a percentage defined by the number of cells that meet the optimum level 
in the source and/or sink set. The package redistributes the ground water levels/ surface 
water levels from source to sink cells that meet this condition. 


Reinjection Drainflow Package 


The Reinjection Drainflow (RDF) Package was created in 1999 for the SFWMD 
by McDonald Morrissey Associates, Inc. to simulate a drain from which flow could be 
directed into another cell rather than being permanently removed from the groundwater 
flow system, as is the case in the original Drain Package (Jones 1999). The reinjection 
option is invoked by adding the destination (layer, row, and column) of the cell to which 
the drain flow is directed. The user specifies both a source destination and a sink 
destination. The destination in the RDF is similar to the source and sink concept from the 
Diversion Package (Restrepo et al. 1998); however, the source and sink destinations in 
the RDF have a 1:1 relationship while the Diversion Package source and sink sets have a 
many to many relationship. If no destination cell is specified or if the destination cell is 
inactive, the RDF Package functions exactly like the Drain Package. The reinjection rate 
is calculated as the product of the drain conductance and the elevation of the head in the 
cell about the elevation of the drain. The RDF Package can be used, for example, to 
represent the removal of water from a reservoir and dispersal onto a nearby wetland or 
injection into an ASR well. 


Enhancements for the RDF Package were made by FAU (Restrepo 2003) to 
incorporate a second level for the sink destination cells. The purpose of this second level 
is to restrict water moving into the sink cells by setting a maximum level. If the heads in 
the sink cells are below the maximum level, the discharge will be allowed only of the 
source cells are above the minimum. The benefit of using this approach is that the user 
can vary the operations on a daily basis, as in an operational schedule, by varying the 
control levels. 


Horizontal Flow Barrier Package 


The Horizontal Flow Barrier (HFB) Package was developed by the USGS in order 
to simulate thin, vertical low-permeability barriers to horizontal groundwater flow. These 
barriers may be natural geologic features, such as vertical faults; or may be man-made, 
such as slurry walls, sheet piles, or levees. Prior to the development of this package, 
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vertical low permeability features were modeled by reduced model grid spacing, or by the 
introduction of variable grid spacing in the region of concern. The HFB Package allows 
these features to be represented without increasing the number of model cells, thereby 
improving efficiency (Hsieh and Freckleton 1993). 


Figure 45a shows a plan view of a thin, vertical low-permeability geologic 
feature in a model layer. This feature is represented in Figure 45b as a series of seven 
horizontal flow barriers located on the boundary between adjacent cells in the finite- 
difference grid. 


ow 


.] permeability 
feature ` 





A B 


Figure 45. Representation of a low permeability feature (A) and its simulation in the MODFLOW 
HFB Package as a series of horizontal flow barriers within the model layer. Source: 
Hsieh and Freckleton 1993. 


The HFB Package makes the assumption that the width of the barrier is negligibly 
small in comparison with the horizontal model cell dimensions. Barrier width is not 
explicitly defined, but is included implicitly as a hydraulic characteristic either as (1) 
barrier transmissivity divided by barrier width in a layer with constant transmissivity or 
(2) barrier hydraulic conductivity divided by barrier width in a layer with variable 
transmissivity (Hsieh and Freckleton 1993). The HFB Package was not used in model 
calibration runs for this model; however, it is used in predictive simulations, such as in 
modeling the L31N Seepage Barrier Pilot Project (Model Application Group 2003). 


Trigger Package 


The Trigger (TRG) Package was developed in 1992 by the SFWMD and Water 
Resources Management, Inc. in order to initiate cutbacks in well pumpage in response to 
head declines in specified “trigger cells”. It was originally created for use in the SFWMM 
and was modified later in 1999 to work with MODFLOW-96 (Randall 1992). 
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The Trigger Package employs two different methods to reduce pumpage: the cell 
mechanism and the time series mechanism. The cell mechanism triggers cutbacks when 
water levels fall below user-specified limits at selected cells, while the time series 
mechanism triggers cutbacks at user-specified time periods, regardless of water levels. 
Trigger cells are located in areas of concern, and are associated with larger, rectangular 
zones within the model. A combination of both the cell and the time series mechanisms 
determines which zones will be subjected to cutbacks. 


In the cell mechanism, a violation is signaled when the head in a trigger cell falls 
below the specified trigger value. Cutbacks in pumpage are calculated and applied in 
phases based upon the severity of the violation. The user may require that the violation 
must occur for a certain number of consecutive stress periods (called the “delay period”) 
before a cutback is applied in order to avoid overreacting to short-term dips in head. 


The time series mechanism applies cutbacks according to specified time periods 
such as the dry season, or periods when low Lake Okeechobee levels are planned. These 
cutbacks remain in effect until the end of the dry season, or until a more severe cutback is 
triggered. Violations occur in phases, with Phase 1 being the least severe by convention. 
Dry season cutbacks without associated low trigger cell levels are always Phase 1 level 
cutbacks. If there are trigger cells in a zone having different levels of violation, the cell 
with the most severe violation determines the cutback phase for that zone. 


Pumping reductions are applied to different wells according to their usage type. 
Some examples of usage types are urban, agricultural, golf, and “zero” type, which 
means that the well is exempt from cutbacks. When pumpage reductions are triggered, 
the program outputs two files: one showing where the triggers were set and what phase 
was assigned, and another showing which wells will have reduced pumpage, and by what 
amount. 


In 2005, source code modifications to the Trigger Package were requested by the 
Water Use Division at the SFWMD to better reflect existing SFWMD water shortage 
policies during drought periods. It was necessary to modify the source code by adding a 
constraint to the logic of the cell mechanism. In south Florida, the wet season lasts from 
June Ist through October 31st. Historically, the SFWMD's Governing Board has not 
made a water shortage declaration at the end of the dry season (i.e., April and May), since 
the wet season would be approaching to bring the needed rainfall. Therefore, an 
additional user-defined interval or “no cutback window” was programmed to specify 
times when no cutbacks would occur. It was recommended that the “no cutback window” 
specification (in the input data) begin on April 1st and end on October 31st. 


The Trigger Package was not used for calibration runs because pumpage 
reductions were represented by actual pumpage records obtained from the SFWMD's 
Water Use Regulation records. This package is used for predictive scenarios. If the model 
is re-calibrated, the trigger package must then be re-calibrated as well. 
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Well Package Modifications and Additions 


Rather than modifying the original MODFLOW Well (WEL) Package, a special 
Trigger Well Package was created by Water Resources Management, Inc. to be used with 
the Trigger Package. The Trigger Well Package has an added field called “tipe” (the 
name “type” could not be used, as it is a FORTRAN keyword) which associates a usage 
type with each well (e.g. urban, agricultural, etc.). It also contains a temporary storage 
place for input pumpage. This allows the reuse of the previous stress period’s pumpage if 
that pumpage had been reduced by the trigger mechanism. 


Both the original Well Package and the Trigger Well Package were modified by 
the District staff in 1999 to allow wells to be read from multiple files. This is useful when 
changes are made frequently to certain types of wells (1.e. public supply wells) while 
other well data remains fairly static. The primary well file will allow up to two additional 
unit numbers to be included. These changes allow the reuse flag (-1) to be invoked 
separately for each file. For example, the first file may have 500 wells, the second 20, and 
the third 10 for the first stress period. In the second stress period, one might decide to 
reuse all the wells in the first two files, but specify 40 new wells in the third file. 


Utility Generation Package 


The Utility Generation Package (UGEN) was developed by the Center for 
Hydrology and Water Resources at Florida Atlantic University, as a tool to generate time- 
related MODFLOW input during model execution. Traditional MODFLOW input files 
must be built prior to execution for each package and for the entire simulation period. 
UGEN creates the input files “on the fly” by linking the static input parameters, such as 
physical location, with the dynamic temporal data, much like a relational database 
(Restrepo et al. 2003). The static and dynamic data are linked by “location name” 
identifiers. The UGEN package can be used to generate these MODFLOW packages: 
River, Drain, Well, General Head Boundary, Retnjection Drainflow, and Diversion. It 
can also be used to correct heads due to the presence of the salt-water interface and to 
calculate the hydraulic conductance for the River, General Head Boundary, and Drain 
packages. 


When MODFLOW is run with UGEN, three input files are needed. The first, the 
UGEN file, defines the number of stations for stages and/or flows and any flags. The 
second type of input file needed is the observation file. UGEN observation files contain 
stage and/or flow values for each stress period for all stations. Each record in the 
observation file represents one stress period, as shown in Table 6. 
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Table 6. Example of UGEN Observation File Displaying Time-Varying Data 


ID Year Month Day Station-1 Station-2 Station-3 
1 1990 1 1 14.3 7.9 11.1 
2 1990 1 2 14.5 8.1 11.5 
3 1990 1 3 14.5 8.0 10.8 
4 1990 1 4 14.3 8.3 11.0 


The third type of input file is the modified MODFLOW input package file to be 
used with UGEN (e.g. RIV, DRN, WEL, GHB, DIV, or RDF). This input file will 
provide the same information that MODFLOW requires. However, instead of having 
stage or flow values, it will have the location name identifier. This represents the “static” 
portion of the UGEN input files that the utility links with the “dynamic” information 
provided in the observation files. These modified MODFLOW input files may be kept 
constant during all simulations, or may be changed if necessary. 


UGEN is completely compatible with MODFLOW in both its modular structure 
and its programming language (FORTRAN). When activated within MODFLOW, 
UGEN saves both storage space and execution time. The placement of the static data and 
time-series data into separate files also reduces model setup time. Table 7 is an example 
of a portion of aMODFLOW River Package input file. The first two lines have the same 
information a normal MODFLOW River Package input file contains. The third line 
contains the static information that is in the usual input file; however, the stage data is 
replaced by the station name. 


Table 7. Example MODFLOW River Package Input File Formatted for UGEN. 


53 0 
1 
1 75 44 500 1.3 Station-1 
1 121 39 500 -1.5 Station-2 
1 14 208 500 2.4 Station-3 
-1 
-1 Indicates data being reused for the last three stress periods 


-1 
Strongly Implicit Procedure Package Enhancements 


Two alternative enhancements were developed by the SFWMD in 1998 for the 
Strongly Implicit Procedure (SIP) Package in order to improve or maintain model 
stability. Both alternatives have added optional variables. If the variables are not used, the 
SIP package will function normally. 


In Alternative 1, two optional variables are added to the second line of the SIP 
input file. These are HCLOSEMAX and NOSTOP. When the maximum number of 
iterations 18 reached and the maximum head change in a cell is less than HCLOSEMAX, 
SIP continues to the next time step rather than aborting the simulation. This allows a tight 
closure criterion (via the original HCLOSE term) for most of the simulation, while 
tolerating a few problem stress periods. When NOSTOP is included and set equal to 1, 
the program will not terminate if HCLOSEMAX is violated. Instead, the problem cells 
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are reset to their values at the end of the last time step and a warning message 19 written 
to the output file. This is helpful in trying to improve a model with stability problems. 


In Alternative 2, four optional variables were added to the SIP input file. These 
are MNITER and NITERSL on the first input line, and HCLOSEMAX and DACCL on 
the second input line. MNITER is the minimum number of iterations. NITERSL is the 
minimum number of iterations before deceleration is allowed. DACCL specifies the 
fraction by which the simulation will decelerate. HCLOSEMAX is the same as described 
in Alternative 1. HCLOSEMAX and HCLOSE together serve as an upper and lower 
bound. Deceleration allows the model to iterate slower, thereby helping maintain 
stability. The simulation will terminate if the closure criterion exceeds HCLOSEMAX. 


Output Control - Summation Enhancement 


The MODFLOW summation utility was developed by the SFWMD as an 
enhancement to the output control file. Because of the size and long simulation periods 
used in SFWMD models, and the use of daily stress periods, looking at cell-by-cell flows 
was an overwhelming task. The summation utility allows the user to sum the cell-by-cell 
flows after any number of stress periods so that the output data is more manageable. 


The enhancement was made by adding the variable INUMSUM to the first line of 
the output control file, and by modifying the variable ICBCFL. INUMSUM tells the 
program how many packages there are to be summed, and is specified once per 
simulation. ICBCFL is the cell-by-cell flow flag. Originally if this variable was 0, flows 
were not saved, and if it was 1, flows were written to disk. These two options are still 
available. The enhanced version, however, allows the user the use the number 4 to add 
the current flows to those in an accumulator, and to use the number 5 in the stress period 
when the final flows are added to the accumulator and the summed flows are written to 
disk. Currently the SFWMD practice is to sum the cell-by-cell flows at then end of every 
month. 


Multibud Enhancements 


Multibud was originally developed by the SFWMD as a post-processing tool for 
the subregional models. Due to the Output Control Summation Enhancements, flows are 
summed to a monthly rate (to save disk space). After model execution, Multibud extracts 
the flow (or water) budget from pre-defined cells or regions. For example, 1f Multibud 
was used to extract the flow budget from one cell, Multibud would provide a time-series 
of flows from that cell for the length of the simulation. If a region consisted of 50 cells, 
Multibud would sum the flows from the 50 cells and produce a time-series of the water 
budget. Multibud produces flow budgets for MODFLOW terms, including: storage, cell- 
by-cell flows from all six faces, constant head and general head flows, pumping rates, 
leakage from drains and rivers, ET and recharge rates, flows from water diversions, and 
reinjected drainflows. 
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In 2004, the SFWMD’s Processing Section proposed to incorporate Multibud to 
the MODFLOW-96 source code with the intention of producing daily flow budgets 
without writing MODFLOW’s cell-by-cell file for the entire matrix for the length of the 
simulation. In 2005, with the help of Ecology & Environment, Inc. this idea was 
implemented and verified for use in the MODFLOW-96 source code. 


MODEL DESIGN 


The MODFLOW-96 source code has been substantially altered by introducing 
many add-on packages. For that reason, customized pre- and post-processing tools have 
been developed to aid in the process of model design, as well as model applications. A 
combination of mostly FORTRAN, Arc Macro Language (.amls), and geoprocessing 
programs are used to construct model input and output. In addition, the SFWMD's 
Processing Section has developed a Performance Measure (PM) Viewer which manages 
and graphically displays input and output from multiple simulations. 


Spatial and Temporal Discretization 


The LECsR Model grid is uniformly discretized into 704 foot by 704 foot cells 
(each cell covering about 11 acres) with 1033 rows and 408 columns resulting in a full 
grid of 421,464 cells per layer (Figure 46). The grid limits, from the lower left corner 
node coordinates (Xmin, Ymin) are given as (680961.0, 318790.0) in U.S. State Plane, 
Florida East Zone, NAD 83. 


In the study area, the hydraulic gradient is from north to south from Lake 
Okeechobee to the Everglades; water also moves west to east (towards the coast). A 
larger nodal spacing could be used when the regional water table slopes gradually from 
about 20 ft, NGVD to 1 ft, NGVD. However, changes in the water table towards the 
coastline occur over a shorter distance, partly due to induced hydrologic stresses (e.g., 
pumping wells and canal drainage systems) and could require a smaller nodal spacing. 
For instance, there can be up to 8-10 feet in elevation change between the headwater and 
tailwater of a structure (e.g., S-80, S-155 and G-56). Selecting the 704 foot by 704 foot 
cells was a compromise between the regional (i.e., large nodal spacing) and local (1.e., 
small nodal spacing) hydraulic gradients. Moreover, the cell size of 704 foot by 704 foot 
was selected in order to interface with the SEWMM's 2 mile by 2 mile cells for 
predictive runs. Consequently, there are 225 LECsR Model cells within 1 SFWMM cell. 


The model will be run under transient conditions. A steady-state version of the 
model will not be developed, since the south Florida hydrologic system is very dynamic 
(e.g., wellfield withdrawal rates vary, operational changes over time). The temporal 
discretization was chosen to accurately reflect the hydrologic system changes (e.g., 
rainfall and canal stages) over a recent wet and dry cycle representative of 1-in-10 year 
events. Daily input data were available to construct the hydrologic packages (1.e., ET, 
Recharge, River, Drain, General Head Boundary, and Operations) and excluded pumping 
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stresses. Average monthly pumping rates (fluxes) were used to generate daily 
withdrawals because the pumping records are compiled monthly. 


A temporal discretization of one day was chosen. This discretization will be 
applied as a one day time step in a one day stress period. A daily stress period was 
selected for the following reasons: 


* Groundwater levels have been observed to vary due to rainfall on a 
daily basis. 


* Groundwater levels have been observed to vary considerably (e.g., 
6 inches) due to the individual or combined effects of canal drainage, 
overland flow, and wellfield withdrawals in time periods of less than a 
week, 


° Considerable portions of the model are impacted by drainage systems 
in which the canal stages may change in a single day. 


* The natural areas (e.g., J.W. Corbett Wildlife Management Area) 
provide considerable surface water storage which can not be 
accurately modeled during wet periods (when the groundwater is 
above the ground surface) with a weekly or monthly average value. 


* Evaluation of structural (e.g. pump sizes) and operational rules (stage 
discharge relationships) are more easily and intuitively incorporated 
using a daily time step. 

° Wetland hydroperiods can be highly variable (e.g., inundation periods 
of less than one month). 


This uniform grid and daily stress period allow for an adequate level of 
subregional accuracy with manageable run and post-processing times. A SFWMD 
computer analyst estimated disk space and processing requirements for several model 
grid designs based upon the SFWMM domain. 


Table8. Initial Estimated Processing Requirements for the LECsR Model. 


WMM| CellsWMM| LECsR Total Bytes Total Bvtesj Approx Output) Cell 
Col 36 year daily 36 year monthly Size, Size 
e 
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(65 | 41 | 400 | 
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a 4264000] T 110,522,880,000 |278,609,760,000 


id | NEST 90,846,281,2001 44 433,171,200 | 84,279,452,400 960. 


At the time of model design, the Model Application Section was somewhat 
limited by working across the network and executing models on Sun workstations 
running Solaris 5.8, primarily SunBlade 150's. The standard configuration for these 
Workstations includes 512 MB of RAM, a 550 MHz UltraSparc processor, and they have 
local external disks providing 100-400 GB of storage on each. Four terabytes of network 
disk space became available (via a Storage Area Network (SAN)) months after the initial 


grid design. 





The MODFLOW-96 source code has been compiled using Sun's Forte 
FORTRAN compiler (version 6.2) under Solaris 5.8 and Digital Visual Fortran (version 
6.6) under Windows XP Professional. The compilation is done with an optimization that 
provides closest to the maximum performance possible for this architecture (level -O4 or 
—fast). 


At the present time, the Section is in the process of migrating over to two, faster 
Windows-2003 servers through Terminal Services. With hyper-threading activated, there 
are 8 processors (four virtual and four physical) per server. The total amount of disk 
storage between the two servers is four Terabytes. Due to large storage requirements of 
the model, the Section has been asked by the SFWMD centralized Information 
Technology Department to make every effort to reduce its disk storage requirements. 


Overall, the model grid was designed with parsimony by trying to choose a spatial 
and temporal discretization that would represent the hydrologic system well and keep the 
disk storage to a minimum. At the same time, it was critical to model the LEC region 
using a single model. The Section receives more modeling requests for this area than any 
other within the SFWMD boundary. Often, the modeling requests involve system-wide 
changes in the LEC. The best way to fulfill the requests is to have one model, even 
though the computer disk and processing requirements are significant. 


Draft - 86 


LEC subRegional MODFLOW Model Documentation DRAFT 









5 EZ x Z= 
^ Ve 나 
D EE L 
~ , 
I 
em) un im mus mu. Mm: en RA ` 
1 Ç 
——-4 —— 1-2 ; < 
Y 
4 ㆍ w A. 


l 
kj—t 
= 


E Sa Nees Every 10th Model 
eee U 0 Row 
ae EE . gu — ——— Every 10th Model 
KE ECHOBEE pt A Wm t Column 



















——À49——— — $—— dono o c bre A q 


PALM BEACH... 





— — — — Se A L 

"-— SWAT 진주 a — — — 998 29 E — +> 
=) ac = ome = = w. LI 
| — 1. ba 


A pe jew = —— P ——— + —— emma 





바보 y +Y 





D E R T T LS. Se L TL LOL. LLLI LIA b LLL 
AA ALA --+-BROWARD ` SSS SS SS SSS SSS eee 
bach l e Ee Ee E IL 


— + —— — —+—— s - py oad 





mtm ptt meme n nan em E EE 





———————————————————————————————————————— l) 






e -~ e EE --— ~- <e e- jamrah t———-—— + -<-—* — - . 
8, R E. NEN! - + zmn. 


mI s sss e "mm — T MU ` ven J i — 





r————— M P—— X ta —— As dw 





la 2 oe bdo O ME — — ——n 


E 


e ' 4 , 
| L l 
Y U + 
| m man am em em "mm 
~… ~ E wi ie B ui E os en 100 sess 1 
y I 1 t i 1 ㅣ 4 i U U i 


b 41 5 E 4 
H » = $ + 
9 - ` bag on 
Ue Ft Fas p pr e nn qon b DPS RR qa 


œ ` $ tn $ 
P "f acci Xi 
5 h ve | - > 
Tei le n n beten ST fu - 
"M. “ark du ` KT fale ^ ] < 45 1 l ex (DA ID Ñ 
000 k ETATIS M hey RUM Mg tes KK 1 Mi 
Y^ ZA Ae iN EE CN vw GK OË, et ite T^ ㆍ ` 





RTE? E E ET EE ep en op tt 


i 
LA 









A i $ d $ SS i i 1 d 
I o pod <> 

} $. l i? Í i 

+ e jo h e ¿—+ L L 









Cee ENSE RS RN ei (n MESI a UNE EET NHAT 
MONROE pM J H H ++ m. po joo qe mn ———- ge E: : L p Je š i U : ¡A A yd 
$ 1 i j d i i 42 V 1 i ` OE: > $ 
TIT ak LIITTTTTT +) ft. EF j L ㅣ f 해 E FRK I ES A 
Ai ġe a Oe O eS 4 TO O OS oe ee 10 oe H ` ċ Gë : 











Miles | Lower East Coast 
E subRegional Model 








Prepared by : MAPU 
Date: 03/04/2006 
Map Doc.:modelmesg_06.mxd 






Figure 46. Model Mesh (Displaying Every 10" Row and Every 10" Column). 
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Groundwater Flow System 


The LECsR Model contains three layers comprising the entire Surficial Aquifer 
System (SAS), providing good representation of the hydrogeologic zones within the 
aquifer system as well as the partial penetration of canals and wells. The layering scheme 
incorporates the two principal permeable zones targeted by production wells (i.e., a very 
high transmissive zone representing the prolific Biscayne aquifer and a more laterally 
extensive moderately transmissive production zone representing the Gray Limestone 
aquifer). 


Methodology for Model layers 


Initial Layer Design 


Based on the conceptual model, the LECsR Model was vertically discretized into 
three layers. With three, active layers the LECsR Model allocates space for roughly 1.2 
million cells; about 30 percent of these cells are inactive. Appendix A contains geologic 
wells with x and y coordinates, land surface elevations, hydraulic conductivities for each 
unit, lithology and thicknesses for the Holocene, Quaternary era ( Q5, Q,4, Q3, Q2, Q1), 
and Tertiary era (T2, T1). Figure 26 displays the locations of the wells. The wells are 
distributed over the model area with gaps mainly seen in the Water Conservation areas 
(WCA’s) and in the Everglades Agricultural Area (EAA). Units were grouped based on 
similar hydrostratigraphic properties (e.g. depositional environments) Model layers 
consist of the following units: 


* Layer 1 - Holocene+Q5+Q4, 
* Laver 2 - Q234 Q24 Ql, and 
° Layer3- T2+T1. 


The thickness at each well was calculated for each model layer. To obtain the top 
of layer 2 the thickness of layer 1 was subtracted from the ground elevation. The top of 
layer 3 was calculated by subtracting the thickness of layer 2 from the top of layer 2. The 
bottom of layer 3 was obtained by subtracted the thickness of layer 3 from the top of 
layer 3. Inverse distance weighting interpolation was used to create the top and bottom 
surfaces. 


Fence diagrams showing the layering scheme are based on the cross-sections from 
Figure 18. The layers thin out in the western portions of the model (Figures 47 and 49). 
Thinning out was not a problem in the north-south direction (Figure 48). Figures 47 to 
49 show the thickness based on the IDW interpolation. In some areas the layer thickness 
needed to be increased for model stability. Near L-47 and L-65 the interpolated top of 
layer 2 was higher the ground surface (top of layer 1) which is outside the active model 
area, in these areas the elevation of top of layer 2 need to be lowered. The changes to the 
model layers are discussed below. 
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Figure 48. North-South Cross-section B-B' with Model Layers 1 to 3. 
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Figure 49. East-West Cross-section D-D' with Model Layers 1 to 3. 


Maximum Layer Elevations 


In some locations interpolation by IDW generated values for the top of layer 2 
that were higher then the ground surface (these areas were outside the active model area). 
In other areas layers 1 and 2 thin out, as discussed above. For model stability, maximum 
elevations were applied to the model layers. 


Although the layer thickness was not adjusted to reach a better calibration, 
modifications to the layering scheme were necessary to prevent drying of cells. In the 
southern and central portions of the model, layers 1 and 2 become extremely thin. As a 
result, during periods of drought, cells along the western boundary tended to go dry. In 
order to overcome this problem and produce a more stable model, the maximum value for 
the top of layer 2 was set to 0.0 feet NGVD. The 0.0 NGVD value was chosen because 
cells in the extreme southwestern portions of the model are subject to tidal boundary 
conditions. In addition, the top of layer 3 was set to -10.0 feet NGVD because numerous 
cells had the bottom of layer 2 above 0.0 feet NGVD in the conceptual model. The 
bottom of layer 3 was not corrected; it was determined from the contact of the Hawthorn 
Group 


The topographic surface generated from the 100-foot DEM was developed after 
the model layering scheme. The thickness of layer 1 was originally determined from the 
geologic control wells. The values at each of the geologic well points should be similar in 
both the topographic surface and top of layer 1 surface. Figures 50 to 52 show the model 
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layers as applied and Figures 53 to 55 illustrate the thickness of each layer. The model 
layers were examined to assure that none of the layers crossed (or overlapped). 


Layer Properties 


Layers 1 to 3 represent the whole SAS; there are no unrepresented confining 
units. Layer 1 1s set to “unconfined” in the BCF package, and contains all river, drain, 
wetland, diversion, reinjected drainflow, recharge, and evapotranspiration cells. Layer 1 
extends from land surface to an elevation of -60 feet NGVD incorporating the 
unsaturated/saturated zone contact. The bottom of layer 1 was set sufficiently deep to 
prevent drying of cells especially in areas with variable hydroperiods. 


Layer 2 is set to “confined/unconfined” in the BCF package and represents the 
lower permeable sediments within the upper SAS and the higher permeable limestone of 
the Biscayne aquifer. This layer extends from about -10 feet to -142 feet NGVD. The 
bottom of this layer resides in the Biscayne aquifer. 


Layer 3 1s set to “confined/unconfined” in the BCF package and represents the 
Gray Limestone/Lower Tamiami aquifer within the lower SAS. This layer extends from 
about -67 feet to -246 feet NGVD. The bottom of this layer resides at the top of the Peace 
River Formation. 
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Figure 50. Top of Layer 2 in feet NGVD1929. 
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Figure 51. Top of Layer 3 in feet NGVD1929 
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Figure 52. Bottom of Layer 3 in feet NGVD 1929 
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Figure 53. Thickness Layer 1 in feet NGVD 1929 
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Figure 54. Thickness of layer 2 in feet. NGVD 1929 
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Figure 55. Thickness of Layer 3 in feet NGVD 1929 
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Hydraulic Conductivity 


The hydraulic conductivities (Hy) for each well in each layer were calculated by 
first computing the transmissivities (T) for each unit (Holocene, Q5 to QI and T2 to T1). 
The transmissivity is calculated by multiplying the hydraulic conductivity for each unit 
by the thickness. The hydraulic conductivities are determined by the sum of the 
transmissivities divided by the total layer thickness. 


° Layer 1 Hy = (T Holocene + T Q5 + T (04) / (Thickness of layer 1) 
* Layer 2 Hy = (T Q3 + T Q2* T Q1) / (Thickness of layer 2) 
° Layer3 Hy = (T T2 + T T1) / (Thickness of layer 3) 


The hydraulic conductivity surfaces were created in the same manner as the 
geologic surfaces, using Inverse Distance Weighting on the points. The hydraulic 
conductivity for model 1 (composite of Holocene, Q5 and Q4) is shown in Figure 56. 
The hydraulic conductivity values for the Holocene are low (below 50 fee/day) for the 
whole LECsR area. In Q5 and Q4 test values of above 10,000 feet /day were measured 
for Miami — Dade County. The hydraulic conductivity for layer 1 ranges from 20 to 
18,500 ft/day. Low hydraulic conductivities are observed in the northern counties and 
much higher hydraulic conductivities in Miami-Dade. 


The methodology for computing the hydraulic conductivity in layer 2 varied 
slightly from the other layers as data from two additional sources were available. The 
hydraulic conductivities for the units Q3-Q1 were derived solely from the geologic 
control wells as displayed in Figures 32 to 34. The locations for the geologic control 
wells, specific capacity tests and APT's used in layer 2 are shown in Figure 26. The 
additional data for hydraulic conductivity was based on specific capacity tests and aquifer 
performance tests (APT) from the main production zone. Figure 57 shows the hydraulic 
conductivity values for layer 2, which ranges from 1 to 75,000 ft/day. The data from the 
APT and Specific Capacity increased the number of data points along the coast 


The hydraulic conductivity for model 3 (composite of Holocene, Q5 and (24) are 
shown in Figure 56. The hydraulic conductivity values for periods T2 and TI are 
displayed in Figures 35 to 36. The hydraulic conductivity values for Layer 3 are much 
lower then layers 1 and 2. Layer 3 1s the thickest of the three layers at most of the well 
locations seen in Figure 55. 
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Figure 56. Hydraulic Conductivity of Layer 1 in feet/day 
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Figure 57. Hydraulic Conductivity of Layer 2 in feet/day 
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Storativity 


The specific storage is the volume of water released from storage per unit change 
in head per unit change of aquifer for a saturated formation (Fetter 1988). The specific 
storage is generally very small. A specific storage of 1x 10? was used as determined from 
APT's. Storativity or storage coefficient is estimated as: 


S=5.010°5 


Where b is the thickness of the layer. The storativity values for layers 2 and 3 
were calculated based on the layer thickness. The storativity ranges from 0.85x107 to 
0.65x10° in layer 2 and 0.14x10° to 0.85x10° in layer 3. The distribution of the 
coefficient can be observed in the thickness maps (Figures 55 and 56). 


Specific yield generally ranges from 0.1 to 0.4 (Anderson and Woessner 1992) 
and can be derived from APT’s. The specific yield in layer 1 was set to the constant, 
0.20. Due to the small range of this parameter, it was varied during sensitivity analysis. 


Vertical Conductance Coefficient 


In MODFLOW vertical flow between layers is controlled by the vertical 
conductance coefficient (Vcont), which is a composite term expressed in units of 1/day. 
Vcont 18 an expression of the vertical conductivity in confining unit and the thickness of 
confining unit in that model cell (McDonald and Harbaugh 1988). It is calculated for the 
two nodes located at vertically adjacent hydrogeologic units (i.e., layers) using the 
equation (McDonald and Harbaugh 1988): 


1 

ji A 

om Az,/2 Az,/2 
K K, 


ZH 


Equation 2 





where z, and zare the thickness of the upper and lower layers (ft), Kzu and Kz are 
the horizontal hydraulic conductivities for the upper and lower layers (ft/day) and K,, is 
the hydraulic conductivity for the confining unit. 
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Figure 59. Vertical conductivity coefficient Vcowr 12 (1/day) between layers 1 and 2 
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Figure 60. Vertical conductivity coefficient Vcowr 23 (1/day) between layers 2 and 3. 
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In previous studies, vertical anisotropical factor ranged from 0.5 to 0.1 (Restrepo 
et al. 2001, Nair et al. 2001 and Wilsnack ef al. 2000). The vertical anisotropical factor 
was set to a value of 0.05 due to the extremely high horizontal hydraulic conductivities in 
the Biscayne aquifer. This represents a 1:20 relationship between the vertical and 
horizontal conductivity for each cell. 


Surface Water Flow System 


The surface water flow system in the LECsR Model is simulated using a 
combination of standard and add-on MODFLOW packages. This approach introduces 
complexity into the model and distinguishes it from being a purely groundwater flow 
model. The primary packages that allow the user to represent general operational 
mechanisms in the south Florida surface water system are Diversion and Reinjection 
Drainflow. These packages complement the Wetland Package by introducing non-linear 
flow (1.e., Kadlec) in the upper model layer. 


Topography 


Input to the Wetland and ET Packages are derived from topography. Input to the 
ET Package includes the ET surface elevation - the water table elevation at which 
maximum ET occurs. The Wetland Package identifies the boundary between the surface 
and groundwater systems which corresponds to the top elevation of muck/peat sediments. 
Though topography is not explicitly entered as input, it does function as a surface at 
which fluxes occur; therefore, developing accurate, high resolution topography is 
beneficial to modeling fluxes. 


Data sources shown in Table 1 were incorporated into the model topography by 
assigning each model cell with an average value derived from a 100-foot Digital 
Elevation Model (DEM) constructed for the LEC. The number of points used to calculate 
model cell averages varied depending on data resolution. The quality of data varied also; 
however, this data set is considered to be the best available data in the LEC (Hinton 
2004). The result was a grid with topography for the 704 ft by 704 ft cells. The next step 
included overlaying line features of levees. Levee top elevations were identified from the 
C&SF Project as-built cross-sections for primary canals. Additional levees were 
identified and added after reviewing satellite imagery and consulting with SFWMD staff. 
Since model cells are 704 ft wide, most often the levee and canal are within the same cell. 
Figure 62 illustrates the land surface elevations and associated levee system within the 
LECsR Model. 
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Figure 61. Land Surface Elevations (ft) in the LECsR Model. 
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Wetlands 


Flow within wetlands was modeled using the Wetland Package (Restrepo et al. 
1998). Wetland areas were identified from 1995 Land Use / Land Cover as shown Figure 
63. Certain areas that are classified as non-wetlands were selected to maintain contiguous 
wetland boundaries. This process required some human interaction with the land use data. 
In the Everglades and Water Conservation Areas, the upland forests (e.g., tree islands) 
were incorporated in the active wetland boundary. Conversely, small (less than 50 acres) 
isolated wetlands were dissolved out of the active wetland boundary. 





Welland-Aquifer Interaction 


Model layer 1 consists of the wetland flow system (or overland flow regime), 
underlying soil/muck layer, and all underlying geological strata down to 0 feet NGVD. 
By incorporating these strata into the wetland flow system, the wetland cells will not 
become dry and require rewetting. The wetland layer includes ponded water with a 
specific yield of 0.9; specific yield of the muck/peat is mainly 0.3. 
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Figure 62. Locations of Modeled Wetlands and Levees. 
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Wetland-Aquifer Parameters 


The data requirements for the wetland model, which are needed to simulate 
overland flow in wetlands include the following: boundary of wetland areas; elevation of 
the wetland soil top surface (e.g. derived from land surface); hydraulic conductivity of 
the muck and aquifer (as a composite), anisotropy ratio, top elevation of the muck, 
Kadlec coefficient, specific yield of the muck and surface water body; depth of the soil 
capillary fringe; Kadlec conductance coefficient for overland flow; specific yield of the 
wetland water body; and coefficients a and B for the Kadlec equation (e.g., as derived 
using the Manning equation). 


The Kadlec coefficient is assigned to each cell based on the land use or vegetation 
type. The suggested values were not derived from field data; they were obtained during 
model design. Numerical difficulties may result when a K larger than 1 million is used. 
However, the model was assigned high K’s where appropriate (e.g., open water, sloughs). 
The highest K in the model is 2.0x10% ft2/day/ftf and was applied to wetland areas with 
sparsest vegetation (1.e., along the southern tidal boundary). A lower bound was 
established based on the highest estimated K (2.0x10° ft2/day/ftB).The lower bound is 20 
percent of the maximum K and was applied to wetland areas covered with the most dense 
vegetation. The upper bound describes flow exhibiting less friction and a higher velocity 
while the lower bound describes flow of a higher resistance. All other types of wetland 
classes were assigned a value of hydraulic conductance between the two extremes. An 
exponential function was used to compute K to interpolate values from the upper and 
lower bound. 


Surface Water Management 


A combination of packages was used to manage or handle surface water flows in 
the model. In combination, the Wetland, Diversion, and Reinjection Flow Packages can 
work together to simulate non-Darcy flow, operational schedules, management rules, and 
simple, lumped routing of sub-basin flow downstream. 


Rivers and Drains 


MODFLOW requires a canal conductance to compute the flux to or from the 
aquifer. Hydraulic canal properties (1.e., cross-sections) were used to estimate river and 
drain conductances. Primary and secondary canal networks were derived from the 
previous subregional models for North Palm Beach, South Palm Beach, Broward, and 
South Miami-Dade. Additional missing or misplaced canals were identified from satellite 
and aerial photos. Canal dimensions were obtained from typical cross-sections in 
permits. As in several cases where no or suspect data existed, field transect 
measurements were made across the canals to obtain the required canal dimensions. 
Controls for some agricultural areas had to be estimated from topography and known 
operational procedures for the crop being grown (ie canals are maintained two feet below 
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ground surface for optimal growing conditions). The final canal network used for model 
design is shown in Figures 63 to 67. 


The three main categories used to represent the canal network include canals as 
rivers, drains, and flows. Canals were assigned as rivers when the control elevations 
change over time and water is delivered from the regional system (e, external source). 
Most of the canals in the model act like rivers. Canals were specified as drains when 
controls do not change over time (e.g., weir). The last category assigns canals to the RDF 
and Diversion Packages in order to route canal flows downstream. 
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Figure 63. EAA Canals 
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Figure 64. North Palm Beach Canals 
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Figure 65. South Palm Beach Canals 
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Figure 66. Broward Canals 
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Figure 67. South Miami-Dade Canals 
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Water Deliveries in LECsR 


The Water Conservation Areas operated by the SFWMD are, for all practical 
purposes, above ground impoundments. They are completed surround by a levee system 
and water is moved into them, through them and out of them to meet both urban and 
environmental demands as well as a place to provide storage during wet periods. Each 
compartment is regulated with an upper limit to protect the levees from breaching. 
Lower limits are also addressed in some areas for environmental concerns. However, 
they do occasionally dry out and wildfires are not an uncommon occurrence. Another 
above ground impoundment is the West Palm Beach Water Catchment Area, or Grassy 
Waters Preserve. This impoundment stores water for the City of West Palm Beach and 
utilizes surface water for its water supply. Most of the time, the Catchment Area does not 
receive water supply deliveries from Lake Okeechobee. 


Simulation of above ground impoundments is accomplished primarily through a 
combination of the Wetland package and the Diversion package. Grassy Waters Preserve 
will be used to illustrate how these facilities are simulated. The cells that Grassy Waters 
Preserve are initially identified and assigned as wetland cells. The levee system (Oe, top 
elevations) that surrounds the impoundment is added to the topography in the wetlands. 
The higher topography associated with the levees ensures that water does not overland 
flow to the adjacent cell that are also wetlands but not part of the impoundment. 


The Diversions Package is required to move water into and out of the 
impoundment. In the case of Grassy Water Preserve, water is brought into the 
impoundment from Lake Okeechobee which is outside of the model domain. The 
structure component of the Diversion package allows for the assignment of an outside 
source to ensure that mass balance is preserved. Water is moved into the facility, up to 
the volume specified. The criteria in this case are that the water level must be below 18.9 
feet NGVD at each cell. Water is then removed from the impoundment utilizing the 
standard well package at historical water withdrawal rates. 


Inflow from Lake Okeechobee to Water Conservation Areas 1A, 2A and 3A are 
done similar to Grassy Water Preserve however, outflows to the urban areas and 
movement through the various compartments is accomplished in a slightly different 
manner. The net flow for a conservation area is equal to the Lake Okeechobee inflows 
minus the urban outflows. If the Lake Okeechobee inflows are greater than the urban 
area outflows the net volume for the Conservation area will be positive indicating that 
water is entering the WCA, if it is negative then water is being removed from the 
impoundment. This is done for each individual day of the simulation period. Urban 
outflows from Water Conservation Area | are primarily to the Lake Worth Drainage 
District, Broward County and the City of Boca Raton. 
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Figure 68. Modeled Surface Water Management Features. 


Boundary and Initial Conditions 


Model boundary definition is an important consideration during model design, 
since boundaries can affect flow patterns. In order to minimize errors in model 
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simulation, two principles were followed to define the type and placement of boundaries: 
1) the observed physical system dictated the type of boundary conditions employed, and 
2) where possible, placement of model boundaries was made at great distances from the 
areas of interest in order to reduce the boundary effects upon the simulation results. This 
resulted in a series of general head boundaries (apart from the no-flow inactive cells) 
being applied throughout all model layers as depicted in Figure 70. 


Outer Boundary Conditions 


According to McDonald and Harbaugh (1988), a general head boundary consists 
of a water source outside the modeled area that supplies or removes water to a model cell 
at a rate proportional to the head difference between the source and the active cell. The 
rate at which water is supplied to a cell is given by: 





O=C(H, —h.) Equation 3 


where Q is the flow rate (f /day) to or from the cell from the boundary and C is 
the constant of proportionality boundary. H, is the average head at the source boundary 
and h, is the average head in the cell. 


The General Head Boundary (GHB) Package is applied at all of the cells located 
along the LECsR Model boundary. As described in Chapter 2, the study area was 
delineated by hydrologic and hydrogeologic boundaries. Surface water stages define the 
northern and western boundaries. The eastern and southern boundaries are defined by 
modeled tidal elevations. The GHB Package was chosen so that the conductance could be 
calibrated along the saltwater interface and boundary conditions could vary with time. 


Boundary conditions were applied to the model by identifying several stations 
(Figure 70). Historical or modeled water levels were applied to a section of cells. For 
instance, historical water levels at station, S-80 HEAD were assigned to over 100 cells 
representing the upstream canal stages in C-44. Modeled tidal elevations were derived 
from the hydrodynamic model, WNAT (Hagen 2005). Modeled data was preferred over 
the limited historical, tidal data from NOAA. The modeled tidal data characterize the best 
representative daily value while minimizing the number of stations used in the GHB 
Package. 


Surface water stages include those along the C-44, L-24, and L-28 Canals. All of 
these stages are applied to model layers 1, 2 and 3. Since the aquifer is unconfined, there 
is a good connection between the aquifer and canal stages. In the case of C-44, this canal 
acts like a groundwater divide in which groundwater flow from the north and south is 
intercepted by the canal. The other two canals, along the western boundary, are located in 
a region with low hydrologic stresses (e.g., no pumping wells). This region is also near to 
a regional groundwater divide. 
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The eastern and southern boundaries of the LECsR Model Area follow the 
Intracoastal Waterway, Biscayne Bay, and Florida Bay. Modeled tidal elevations are 
defined by four sections in this area. Initially, astronomic tides were generated for 134 
locations. Individual locations were grouped into sections based on similarities in tidal 
range, phase, and composition (Hagen 2005). 


Sections 1 to 4 are within close proximity to the saltwater interface. The density 
correction option in the UGEN Package 19 activated in order to compute equivalent fresh 
water heads for these cells during model execution. The following formula was used to 
calculate the vertical pressure distribution along the boundary of the freshwater heads: 


7, g 
Hu (H LCD. Equation 4 
J 


Where He, is the equivalent freshwater head at the boundary; H, is the tidal stage; 
Le is the elevation within the aquifer where the equivalent freshwater head is to be 
applied; ys is the specific weight of salt water and yr 1s the specific weight of freshwater. 
If Le is -320 ft, the equivalent freshwater head for a tidal elevation, H, of 2 ft is 2.3 ft 
when Le is-10 ft and y, and yr are 1.021 and 1.0, respectively. The equivalent freshwater 
head for the same tidal elevation at a depth, L. of -70 ft is 3.5 ft (holding all above 
assumptions true). When Le is -220 ft, NGVD, the equivalent freshwater head is 6.6 ft. 
This case (Figure 69) illustrates that when equivalent freshwater head is computed at 
depth there is an upward vertical gradient, which is expected along the saltwater 
interface. 


Prior to model execution, conductances are computed as input to the GHB 
Package. In order to avoid large vertical flows at the boundary while mimicking the 
natural processes at this location, the following scaling factor was applied to the 
conductance values along the tidal boundary: 


2 
Sc = gon — d x3 Equation 5 
(H, iu Ly) 


where Sc is the scaling factor and La is the elevation of the Surficial aquifer 
system base at the node center. The concept of the scaling factor was applied by Restrepo 
et al. 2001, Wilsnack et al. 2000, and Nair et al. 2001 in previous mulit-layered 
groundwater models. The factor is applied to all GHB cells in layers 1 to 3 to 
approximate reduced flows from the saltwater interface into the aquifer at depth. In the 
upper layers of the model, the factor facilitates the upward movement of flow from the 
deeper layers. The movement of freshwater upward along the saline interface toward the 
near surface strata is consistent with previous studies ((Kohout 1964) and Geotrans 
1988)). Since MODFLOW assumes a constant density of water, the scaling factor reflects 
the assumption that the saltwater interface is fixed in time and space. Table 9 shows the 
average conductance for each station and for each layer. 
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Table 9. Representative GHB Conductance for Modeled Tidal Stations. 
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Figure 69. Diagram of Model Layer Used to Compute Head at Depth. 
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Figure 70. Location of Active Model Area and General Head Boundary Stations. 
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Initial Conditions 


Transient models require initial conditions which are boundaries in time. These 
conditions are given at all locations in the model prior to model execution. In this case, 
initial conditions are provided for the three model layers and are representative of the 
water table. 


The SFWMD version of MODFLOW-96 was developed for transient conditions; 
therefore, the user is prevented from executing a steady-state run. During initial model 
design, initial water table elevations were approximated by reducing the topography by 2 
ft in each model cell. The same initial conditions were used in every layer. The 
justification is that the aquifer is unconfined and the water table is shallow and close to 
land surface in many places. It was observed that in higher relief areas (e.g., coastal 
ridges) the initial conditions were too high and in lower relief areas (e.g., marshes) the 
initial conditions were too low. 


The initial conditions were then improved upon by applying pseudo-steady-state 
conditions. A pseudo steady-state run is one that achieves steady-state conditions by 
repeating (daily) stresses for a length of time. The initial conditions based on the lowered 
topography were used as input to a pseudo-steady-state run. All packages with time- 
varying heads or flows were set up to repeat the first stress period. For example, if a well 
pumped 200 cfs in the first stress period, this value was repeated for the next 10 stress 
periods. Pseudo-steady-state conditions were applied for 10 stress periods. The resulting 
heads for layers 1 to 3 were then applied as initial conditions to a second run. This 
process was repeated once more with those results being used in the model at the present 
time (Figure 71). 
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Figure 71. Initial Conditions (ft NGVD) in Model Layer 1. Stretched to 1 Standard Deviation. 
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Evapotranspiration and Recharge 


The calculation of maximum potential evapotranspiration (ET) rate in the LECsR 
is based on reference crop ET which is adjusted according to crop type, available soil 
moisture content, and location of the water table. Algorithms used to calculate actual 
evapotranspiration vary geographically because of different data availability, and varying 
physical and operational characteristics of different areas within the model domain. 


Recharge and maximum potential evapotranspiration (ET) rate time series are 
computed using an ET-recharge model (Restrepo and Giddings 1994). The methodology 
is based on a daily unsaturated/saturated water balance approach on each unique 
combination of land cover and soil type. This model is an extension of the Agricultural 
Field-Scale Irrigation Requirements Simulation (AFSIRS) Program (Smajstrla 1990). 


ET-Recharge Model 


In the LEC and Everglades areas, irrigation supply, unsaturated zone ET and 
recharge into the Surficial Aquifer System are preprocessed, i.e., pre-calculated 
quantities, and used in the unsaturated zone moisture accounting. These quantities, 
among others, were output from the ET-Recharge model (Giddings and Restrepo 1995). 


The unsaturated zone is treated as a separate control volume where infiltration, 
percolation, evapotranspiration and changes in soil moisture are accounted for. The 
reasons for the unsaturated zone accounting are the need to quantify LEC irrigation 
applied to the unsaturated zone and to more accurately assess changes in irrigation 
requirements associated with changes in land use. 


The ET-Recharge model is based on Land Use/ Land Cover, soil, and a reference 
table for each unique land use and soil type area in relating land use classification to the 
following: runoff coefficients; crop type; growing season; percent pervious area; switch 
indicating if area is irrigated or not; and water use type classification. As outlined in 
Restrepo and Giddings (1994), a daily water balance was conducted on each unique 
combination of land cover and soil type. A corresponding maximum evapotranspiration 
rate was assigned to each of these combinations. The daily water balance analysis yields 
the amount of recharge to the water table that is input into the model. The ET and 
recharge values that were developed this way were held relatively constant in the model 
throughout the calibrations process. Figure 72 shows the distribution of recharge in the 
model area. 
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Figure 72. Average Annual Recharge Rate in inches/year (1986-2000) 


Draft 125 


DRAFT LEC subRegional MODFLOW Model Documentation 


Daily rainfall and wet marsh potential ET are defined as inputs to AFSIRS. 
AFSIRS calculates the crop potential evapotranspiration (ET,) using the formula: 


ET, = KeopET, Equation 6 


where Kao is the crop coefficient that varies with crop type and crop growth 
stage, and ET, is the wet marsh potential ET (see wet marsh potential ET section for more 
details). 


The rate at which water is returned from the soil to the atmosphere by 
evapotranspiration is controlled by two factors: atmospheric demand and soil-water 
availability (Jensen e£ al. 1990). For a pervious area, AFSIRS is called to perform crop 
root zone water balance on a daily basis: 


ASTO = RAIN + NIRR — DRAIN — ET Equation 5 


Where ASTO is the change in root zone soil water storage, (in); RAIN is the 
rainfall (in); NIRR is the net irrigation requirement or irrigation supply (in.); DRAIN is 
the drainage and surface runoff (in.); and ET is the evapotranspiration (in). 


AFSIRS calculates irrigation requirements and crop evapotranspiration rates as a 
function of crop type, soil type, irrigation system, growing season, and climatic 
conditions. It assumes that crop requirements are met from the unsaturated zone through 
rainfall or supplemental irrigation. An irrigation management option within AFSIRS was 
selected such that the exact amount and timing of the irrigation is to be used to restore the 
root zone to field capacity (i.e., maximum yield and thus, maximum or potential ET is 
always maintained). 


Irrigation deliveries calculated from the ET-Recharge model are treated as target 
irrigation demands in the LECsR. These irrigation demands can be met from various 
sources (mainly the water table but also from wastewater reuse and public water supply) 
and are the basis for implementing the LEC trigger and cutback modules. 


The calculation of maximum potential evapotranspiration (ET) rate in the LEC 


and Everglades areas can be partitioned into non-irrigated and irrigated areas as explained 
below. Figure 73 shows landscape irrigated and non-irrigated areas in the LECsR. 
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Figure 73. Irrigated and Non-Irrigated Areas in the LECsR Model. 
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ET in Non-irrigated Areas 


Vegetation in the non-irrigated vegetated areas of the LEC utilizes the available 
water in the unsaturated zone and supplements demands from the water table when 
possible. This implies that as water is used up in the unsaturated zone by the plant, the 
plant develops a greater dependency from the saturated zone. The ET in the unsaturated 
zone approaches zero during a deficient rainfall period due to reduction in soil storage. 


For non-irrigated areas in the Water Conservation Areas, Everglades National 
Park and portions of Big Cypress National Preserve, the following assumptions are made: 
(1) moisture content between land surface and water table does not change; (2) ET comes 
only from the saturated zone and/or ponding; and (3) infiltration equals percolation. 


In the saturated zone, the maximum saturated ET is equal to the supplemental 
crop demand. The supplemental demand is equivalent to what the crop needs from the 
water table once the available soil water storage in the unsaturated zone has been 
depleted. In groundwater models (ep. MODFLOW), the maximum saturated zone ET 
rate is a linear function of the distance from the ET surface to the water table and on the 
extinction depth applied for a specific crop. 


ET in irrigated Areas 


The ET which occurs in the irrigated LEC areas is computed essentially equal to 
the ET for the crop because the crop is assumed to be well irrigated. This implies that 
water is removed from the saturated zone at a rate equal to the ET for the crop. If the ET 
rate exceeds the water available in storage in the unsaturated zone, then, water is 
provided to the plant via the irrigation system or rainfall. 


In the saturated zone, the maximum saturated ET is computed by the difference 
between the potential crop ET and the unsaturated zone. The saturated zone ET is equal 
to zero for well irrigated areas because the unsaturated zone can meet all the crop 
demands due to irrigation. 


Deep Percolation Rates 


The deep percolation zone is the zone in which rainfall or excess irrigation 
reaches the aquifer and provides direct recharge to the groundwater system. It is 
calculated as the residual function of the water balance equation within the root zone. 
Due to excessive recharge rates, runoff coefficients are utilized to reduce percolation 
rates for irrigated and non-irrigated pervious areas. 
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Wet Marsh Potential Evapotranspiration 


The SFWMD Simple Method (Abtew 2003) was used to provide estimates of 
long-term historical (1965-2000) wet marsh potential ET for long-term hydrogeological 
modeling as described in Irizarry-Ortiz (2003). The SFWMD Simple Method was 
developed as a modification to the Penman-Monteith formula due to the lack of a 
comprehensive meteorological database for south Florida. The LECsR Model used the 
estimated wet marsh potential ET for the years, 1986-2000. 


Wet marsh potential ET, ET, (mm d-1), is computed as: 
ET, = Kı Rs/A Equation 1 


Where Ki, is a coefficient; Rs is the solar radiation received at the land surface (MJ 
m-2 d-1); and A is the latent head evaporation (MJ kg-1). 


Potential ET is defined as “the rate at which water, if available, would be removed 
from saturated soil in the form of latent head per unit area or the equivalent depth of 
water" (Restrepo and Giddings 1994). Due to the difference in roughness characteristics 
between marsh and reference grass surfaces, the crop coefficients developed with respect 
to a grass reference ET needed to be modified for use with wet marsh potential ET. 


Due to the scarcity of solar radiation and cloud cover data, solar radiation (Rs) 
was estimated as a function of atmospheric transmissivity, temperature, and 
extraterrestrial solar radiation (Irizarry-Ortiz 2003). 


Extraterrestrial solar radiation is calculated from latitude and time of year by 
integrating the instantaneous radiation intensity at the outer atmosphere from sunrise to 
sunset. The potential ET for wetland marsh was calculated and applied at 17 NOAA 
stations with long-term daily temperature data in order to provide estimates of Rs for 
hydrogeologic modeling. Annual time series and summary statistics of wet marsh 
potential evapotranspiration estimated at 17 NOAA stations are presented in Appendix 
A. The Inverse Distance Weighting method was used for spatially-interpolating the wet 
marsh potential ET across the LECSR grid. 


Evapotranspiration Package in MODFLOW 


The ET Package uses the following input data: 


e An ET surface array depicting the elevations at which 
evapotranspiration from the water table occurs at a maximum rate. 


° Anarravof maximum ET rates. 


ㆍ An array of extinction depths that represent the water depths below the 
ET surface where evapotranspiration rates from the water table 
become negligible. | 
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The generalized form of the ET function in MODFLOW is (McDonald and 
Harbaugh 1988): 


ET sat = K fact PT eas Equation 7 


where ET, is the saturated ET rate (inches/day), Kg is an adjustment factor, and 
ET sat-max 19 the maximum potential saturated ET (inches/day). 


Ktact or extinction depth varies linearly from ET surface array (ground surface) to 
the bottom of the shallow and deep root zone. This factor is used to simulate the 
diminishing ability of vegetation to use water and is function of vegetation/crop type and 
water table. 


The ET sat-max rate is computed by the difference between the potential crop ET, 
ET», and the unsaturated zone (ET ynsat). These quantities, among others, are output from 
the ET-Recharge model (Giddings and Restrepo 1995). The ET Recharge Model outputs 
the daily ETsat-max rate. Figure 74 displays the annual average ET sat-max rate. Rates in the 
EAA and other agricultural areas appear lower than expected. The reason for this is that 
AFSIRS is accounting for the water removed from the unsaturated zone due to crop water 
requirements. 
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Figure 74. Average Annual Maximum Evapotranspiration Rate (in/year) (1986-2000) 
Groundwater Use. 
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Water Use 


Data for the Well Package were extracted from the individual water use permits 
issued by the South Florida Water Management District and actual usage reported to the 
SFWMD and USGS. The individual permits were divided into two categories: public 
water supply and non-public water supply. 


Public Water Supply 


Annual allocations (in million gallons per year) for public water supply demands 
were obtained from monthly pumpage reports supplied by utility companies. Individual 
well withdrawals were estimated by creating a wellfield distribution percentage and a 
monthly seasonal component for each permit. The wellfield distributions were changed 
three times (1.6., 1986, 1989 and 1992) during the simulation period when there were 
significant shifts in demands. 


The annual allocations vary in the model according to the historical pumping 
records for that year. If the utility does not have a wellfield cap due to a restriction (e.g. 
wetland impacts, saltwater intrusion), the percentage of the annual allocation is evenly 
distributed among the total primary number of wells per permit. Most of the PWS wells 
in the study area are unrestricted. For example, Permit X contains 2 wells; consequently, 
50 percent of the pumpage is assigned to each well. 


A monthly seasonal component was developed using historical pumpage from 
January to December of each year. These values represent the fraction of pumpage 
distributed throughout the year - monthly seasonality. The fractions add to 1.0. If a 
permit’s annual allocation is 120 mgy and the utility reported that 10 mgy was pumped 
each month, the fractions are 0.0.83 for each month. Monthly pumpages are then 
converted to daily values by dividing by the number of days in each month. The daily 
pumping rate was held constant for the month. 


Non-Public Water Supply 


Non-public water use classifications include agriculture, industry, golf course, 
nursery, and recreation areas. The historical (or actual) pumpage was applied when 
available. When the pumpage data were not available, the modified Blaney-Criddle 
equation was used to calculate the crop water requirements and then to estimate the 
irrigation demands based on rainfall, crop type, and irrigation efficiencies. In reality, not 
all of the water pumped for irrigation is used, some of it is returned to the aquifer as 
recharge. 
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CHAPTER 4 
Model Calibration 


Calibration is the procedure by which input parameters are adjusted within the 
model until a reasonable representation of the physical system is achieved. Model 
calibration is considered to be accomplished when the model is capable of simulating a 
set of field measurements within specified tolerances. The goal of the calibration process 
is to change input parameters, such as boundary conductance or horizontal hydraulic 
conductivities, within a predetermined range in an attempt to produce simulated heads 
and fluxes that match historical values. The predetermined ranges for input values are set 
by a combination of site specific data (e.g. pumping tests) or typical values (e.g. specific 
yield of fine sand which typically ranges from 0.05 to 0.3). The calibration period chosen 
for the LECsR model was January 1986 to September 1999, which has a wide range of 
hydrologic conditions ranging from very dry to average to very wet hydrologic periods. 
The September, 1999 time period was chosen as the end of the calibration period due to a 
significant change in the surface water operations of the overall everglades system which 
occurred in the fall of 1999. The transient calibration comprises 5000 stress periods. Each 
stress period is a day and has one time step. 


The model was calibrated to only to transient conditions as discussed in Chapter 
3. The total number of observation wells used in the calibration is 195 which have 
continuous recorders on them. Figure 80 shows the location of the observation wells. 
With the notable exception of the use of monthly demands (which were converted to 
daily averages) for wellfield withdrawals for utility and agricultural demands, daily stress 
periods were used for the hydrological inputs including rainfall, evapotranspiration, canal 
stages, structure flows, and general head boundary conditions (tidal and groundwater 
levels at the boundaries of the model). 


Calibration was achieved primarily by adjusting parameters within pre-specified 
ranges to better match computed water levels and structure flows with the observed 
historical record. The model was initially run utilizing the original model design. The 
more sensitive parameter where then varied to there maximum and minimum values and 
evaluated for statistically significant trends. Max and min values were selected based on 
physical meaning. For example, horizonal hydraulic conductivity was not less than zero 
ft/day and did not exceed 75,000 ft/day (which was the highest value in the Biscayne 
aquifer from pump test). If any noticeable trends were observed by varying any 
parameter, regional changes were applied to the model to improve calibration. The 
model was then ran with the new parameters and compared to the previous “calibrated” 
model to determine the extent of improvement. If significant improve was observed, the 
model was permanently modified to include these changes. This process was repeated 
multiple times by slowly decreasing the minimum and maximum tolerances of the 
parameters extremes until a satisfactory calibration was achieved. At no time were the 
minimum and maximum tolerance values exceed during the process. The final results of 
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the model runs were subsequently evaluated both qualitatively and quantitatively. No 
automated calibration tools were utilized in the calibration process. 


CALIBRATION CRITERIA AND TARGETS 


The mean error, the mean absolute error and the root-mean square error are the 
statistics that are commonly used to evaluate the comparison between simulated and 
observed heads (Anderson and Woessner 1992). Langevin (2001) suggest that the mean 
absolute error “may be the most useful statistic because it provides a true measure of the 
average difference between observed and simulated heads”. However, other authors have 
used either the mean error or root-mean square error as their primary statistical indicator 
(Anderson and Woessner 1992). In addition, the South Florida Water Management 
District has historically used statistical criteria of plus or minus 1.0 feet for 75 percent of 
the simulation period (Restrepo et al. 1992) for the Lower East Coast of Florida. The 
limitation of this criteria is that it is independent of the variability of well. For this 
analysis, statistical indicators +/- 1.0, ME, MAE and RMSE will be considered the 
primary indicators with STD, MIN/MAX and RES as the secondary indicators. 


In order to determine if the LECsR Model is simulating hydrologic and 
hydrogeologic system accurately, a series of calibration targets were developed prior to 
model execution. A series of seven statistical measures were evaluated to assist in 
determining the degree of calibration for the model. These statistical measures were 
applied to the overall model and each individual well. 


For the entire model area, the following seven statistical criteria are calculated: 


e Residual (RES global) = Average over all 195 observation wells of 


RES 
° Standard deviation (STD global) = Average over all 195 observation 
~~ wells of STD 
e MIN/MAX global = Average over all 195 observation wells of 
MIN/MAX 


e +/- 1.0 global = Average over all 195 observation wells of +/- 1.0 


° Mean error (ME global) = Average over all 195 observation wells of 
ME 


e Mean absolute error (MAE global) = Average over all 195 observation 
wells of MAE 


* Root mean square error (RMSE global) = Average over all 195 
observation wells of RMSE 
For each observation well, the following seven statistical criteria are calculated: 


* RES = Percentage of time where the absolute value of the simulated 
residual minus the standard deviation of the residual is less than 25% 
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of the difference between the maximum and minimum observed heads 
for each observation well. 


STD = Percentage of time where the simulated head lies within the 
range of the observed head for that day plus or minus the observed 
head’s standard deviation for the period of record if the standard 
deviation is greater than 1.0 feet for each observation well. If it is less 
than 1.0 feet, than 1.0 feet is utilized. This is a modification of +/- 1.0 
(see below) to allow for those wells which may experience extreme 
fluctuations. 


MIN/MAX = Percentage of time where simulated head lies within the 
range of the observed heads minimum and maximum value for each 
observation well. 


+/- 1.0 Foot = Percentage of time where simulated head lies within a 
plus or minus one foot band of the observed head for each observation 
well. 


ME = the mean error (ME) is the mean difference between measured 
and simulated heads for each observation well. 


MAE = the mean absolute error (MAE) is the mean of the absolute 
value of the differences in measured and simulated heads for each 
observation well. 


RMSE = the root mean squared error (RMSE) is the average of the 
squared differences in measured and simulated heads for each 
observation well. 


DRAFT 


The seven global model criteria should not be used to ensure a satisfactory 


calibration in the model. However, they do represent objective numerical criteria which 
normally are indicative of calibration accuracy. Considering this model is a transient 
model with 5000 daily stress periods, it is not reasonable to expect that ever criteria can 
be met 100 percent of the time for those calibration targets which record if the simulated 
head fell within the predetermined range on a daily basis. Therefore, the calibration 
targets include a percentage, or an absolute number, depending upon which statistical 
criteria that is being looked at. 
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Table 10. Ranges and Targets for statistical calibration criteria. 


Statist Propo 
l sed 
ical taraġ 

criteri g 

a or 
Target 
Global model criteria 
RES 80 % 
STD 80 % 

EH 80 % 

+/- 1.0 n 
Foot 80% 

+/- 
= 0.75 
MAE 0.75 ft 
RMSE 0.75 ft 


Individual observation well criteria 


RES 75% 

STD 75% 
big 75 % 
+/- 1.0 

Foot 75 6 

ME +/- 1.0 

MAE 1.0 
RMSE 1.0 


The statistical criteria listed in Table 10 should be regarded as calibration targets, 
i.e., if the input data, the number of observation wells and the individual observations 
time series allow it, the calibration should be continued until the criteria are met. In 
general the model should describe the average level and the dynamics of the groundwater 
table fairly accurately. The criteria may serve to check if the deviations between 
simulation and observation have been sufficiently reduced. 


‘Soft’ calibration references 


Apart from field measurements, the model may be evaluated from a more general 
view. The ‘soft’ calibration references could include the general shape of the 
potentiometric surface, regional flow paths, aerial photography and ground truthing when 
evaluating inundated areas, and wetland community types and the environments that they 
thrive in. 


The model results are evaluated from the general knowledge and understanding of 
the model area. As no ‘hard’ data in terms of measurements exist, the comparison 
between simulations and the field conditions is qualitative. An example is where a 
known wetland community exists which has a distinct hydroperiod which may indicate 
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that water is above ground for 3 months out of the year. An average year would be 
evaluated to see if the simulated water levels are reasonable for the type of wetland 
community. Another example for soft calibration includes analyzing the typical regional 
flow pattern in the calibration runs. 


Calibration Process 


The calibration process was conducted in a uniform manner. First the model was 
originally run with the conceptual model and evaluated utilizing the calibration criteria 
established. The initial model run was reasonably calibrated considering the model was a 
partial combination of pre-existing calibrated models. However, these previous models 
did not have calibration periods remotely close to the 5000 stress periods utilized in this 
model. Therefore, additional calibration was necessary to achieve the desired goals. The 
first step in the calibration process was to vary all model parameters to there minimum 
and maximum values. Each of these runs were then compared to the original conceptual 
model for a determination of any improvements or degradation to the model calibration. 
For these first set of model simulations, every parameter was evaluated to see 1f the 
model was sensitive to that parameter. If the model was sensitive to that particular 
parameter in either a positive or negative way, 1t was retained for the calibration process. 
Table 11 shows the main parameters utilized in the calibration process and the tolerance 
levels imposed on each parameter. 


Table 11. Tolerance levels utilized during calibration process 


Model Minim Maxi 
Param um mum 
eter Value Value 
Genera Origin 
| Head RES al 
Bound ae Multipl 
ary "i ied by 
Condu alice 100 
ctance 
(ft2/day 109 
) 
Canal 0.1 100 
Hydrau ft/day fi/day 
lic 
Condu 
ctivity 
of the 


Canal 0.1 5 


thickne 
SS 

Specifi 0.2 

c Yield 

for the 0.8 
wetlan 

ds and 

underly 


ing 
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muck 
layer 
Hydrau 25 100,0 
lic fUday 00 
Condu ft/day 
ctivity 


d areas 
General Head Boundary Conductance 


The general head boundary was modified in portions of Palm Beach and Martin 
Counties. The conductance values were decreased by a factor of 100 along Lake Worth 
Lagoon and Jupiter Inlet. Noticeable improvements were observed in wells PB-565, M- 
1024, PB-1642, PB-809 and PB-1639. The reduction in the conductance value is at the 
extreme end of the range suggesting. The reduction in conductance values suggests that 
boundary conditions may influence water levels in this area of the model. Further 
reductions in conductance values in this area showed no noticeable improvement and 
worsen the calibration in several places. A total of six iterations were made to achieve 
the final general head conductance boundary. With the exception of central and northern 
Palm Beach and southern Martin County along the eastern boarder of the model, no 
changes were made to the general head conductance values elsewhere in the model. 


River and Drain Conductance 


Modifications to the river and drain conductance values, and to there assignment 
as rivers or drains, was the most intensive calibration effort undertaken for a number of 
reasons. In the conceptual model run, it became apparent that large numbers of canals 
were missing from the GIS coverages, and that data in the attribute fields for the existing 
canals was either absent or not correct. In addition, some of the existing canals were 
either in the wrong location or had an incorrect structure assigned to them. Initially, 
existing GIS coverages were obtained to fill in the missing canals. Then from an iterative 
process, each area of the model was evaluated further to investigate if additional canals 
were missing or mislabeled. This took numerous iterations to complete. Canal 
conductance values were modified by changing the thickness of the sediments in the 
canal. Initially, the sediment thickness was set to 1. Sediment thicknesses were then 
allow to vary between 0.1 feet and 5 feet which is within the range of sediment thickness 
found in the study area because some canals have near vertical slopes due to the rock they 
are cut into. The hydraulic conductivity of the sediments was initially set to 1 but was 
allowed to vary between 0.1 and 100 which still may be low considering the extremely 
porous nature of the Miami Limestone and Fort Thompson Formation which the canals 
cut into in Miami-Dade, Broward and southern Palm Beach County. In general, canal 
conductance tends to increase from north to south primarily due to Miami Limestone and 
Fort Thompson Formations which are outcropping in the southern portion of the model. 
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Specific Yield of the Wetlands Muck Layer 


The specific yield of the wetlands muck layer was originally set to 0.6. This 
value was chosen to try and simulate those areas where both wetlands and uplands were 
present in a single model cell. The reason for doing this was that the topo tended to favor 
the upland values and the isolated wetlands which are lower would not be accounted for. 
During the calibration process, the specific yield of the muck was allowed to vary 
between 0.2 and 0.8. The simulations determined that the 0.6 was too high for the 
specific yield of the muck and was reset to a more realistic value of 0.3 for the majority 
of the wetland areas in the model domain. 


Root Zone Extinction Depth and ET Surface Elevations 


The root zone extinction depth and the ET Surface elevation controls, in part, the 
amount of evapotranspiration that is removed from the saturated zone of the water table. 
The process for creating the water table ET and water table recharge attempts to account 
for water removed from the saturated zone. Therefore, in areas that are heavily irrigated, 
ET from the saturated zone in low because the potential ET is removed from the 
unsaturated zone because it is constantly irrigated. 


The ET function in the governing equation is a linear function and is dependent 
upon the depth of the root zone and the ET surface elevation relevant to the elevation of 
the water table. If the root zone depth is to shallow and the ET surface to high compared 
to the calculated head at that location, then little or no ET would be removed from the 
model. Conversely, larger amounts of ET would be removed if the ET land surface was 
low and the root zone deep when compared to the head calculated at that cell. 


Root zone depths were decreased in the Corbett Wildlife Management Area, the 
Tequesta area, and southeastern Miami-Dade County. The ET surface was also lowered 
in extreme southwestern Corbett Wildlife Management Area and the West Palm Beach 
Water Catchment Area. Known issues regarding the topography in northern Palm Beach 
County potentially being to high allowed some manipulation of the topography in these 
regions. The ET surface was slightly increased in Everglades National Park along the 
marl uplands to better simulate these features. 


Wetland Hydraulic Conductivity and Overland Flow Elevations 


The wetland hydraulic conductivity of the muck was originally set to values equal 
to layer one, as calculated from the hydrogeology of the units in layer one. These values 
generally ranged from approximately 25 ft/day to over 50,000 ft/day. 
The calibration process found that a value of 5,000 ft/day was generally sufficient to 
calibrate the wetland gages in most of the wetland systems. The one main exception to 
this was Shark River Slough. 


Shark River Slough is a deep slough system located in Everglades National Park. 
It is a broad southwesterly trending arc originating at the upper northeastern portion of 
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the park and terminating at Florida Bay in the southwestern portions of the park. It is a 
continuous deep wetland system interspersed with tree islands and is within a slight 
trough located within the underlying limestone. The system is managed by a series of 
pump stations and numerous culverts located along Tamiami Trail. Due to the low 
topographic relief in the area, the hydraulic conductivity of the muck and the overland 
flow elevations had to be modified significantly to achieve calibration. The hydraulic 
conductivity of layer one was increased to 500,000 ft/day and the topography lowered by 
a foot to help water move through the slough (or preferential flow path). 


Urban Recharge 


The final main changes made to the model during the calibration process was a 
modification to the et-recharge preprocessing program. This program calculates recharge 
to the water table dependent upon the amount of water retained in the unsaturated zone 
with some of the main input parameters including soil type, landuse, and crop type. The 
model originally utilized and standard curve number approach for estimating runoff for 
the various land uses in the study area. A slight improvement was noted in the majority 
of the urban areas when the curve number runoff estimates was not utilized. There are 
two apparent reasons for not utilizing runoff estimates in the model. The extreme amount 
of secondary canal systems, coupled with the high hydraulic conductivity of the Biscayne 
aquifer, rapidly removes water from the model, which, in effect is simulating the 
overland flow component. Second, south Florida has undergone tremendous 
development of the last 20 to 30 years. Developments that have recently been 
constructed are required to have onsite lake systems and surface water treatment facility 
to retain storm water onsite. So the runoff from these newer developments is limited to 
within the development itself and generally does not discharge to the primary canal 
systems except under heavy rainfall conditions. Even when they do discharge to the 
primary canal systems, there is a lag of hours to days where the water 1s retained on site 
prior to discharge. Considering the model is operating on a daily stress period and time 
step, this allows the river and drain package sufficient time to remove the excess runoff. 


Global Model Calibration 
During the calibration process, observations wells with continuous data recorders 
were analyzed. The location of the observation wells are shown in Figure 78. Table 12 


present the statistical results of the over all model calibration. 


Table 12. Calibration Results for the Entire Model Domain (195 observation wells) 


SR Propo Simul 
SCC sed ated 
criteria IUe i 

RES 100 
Global 85% % 
STD Š 80 % 
Global 85 % 
MIN/M 85 % 99 Yo 
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AX 
Global 
+/- 1.0 0 87 % 
Global 85 % 
ME +/- 0.00 
0.75 ft 
Global f 
MAE 0.75 0.55 
Global ft ft 
RMSE 0.75 0.70 
Global ft ft 


The global calibration results indicate that 6 out of the 7 criteria were met. The 
global model calibration was not within 1 standard deviation more than 85 percent of the 
time; however, the target was missed by only 5 percent. Based on these results, the model 
is calibrated under the global criteria. It should be noted that this result in upon itself is 
not indicative of a calibrated model but does provide a reasonable starting point in 
evaluating the effectiveness of the calibration. 


Figure 75 shows a histogram of the distribution of the mean error for each of the 
observation wells. For the entire model the mean error distribution appears to be 
symmetrical distribution with a slight bias to the positive side. This indicates that the 
overall model may be slightly over-predicting heads. The original calibration criteria 
target for was plus or minus 1.0 feet or less. The graph indicates that approximately 95 
percent of the wells meet this condition. 


Mean Error Distribution of Monitor Wells 
Lower East Coast sub-Regional Model 
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Figure 75. Mean Error Histogram for Model Calibration Period 
The histogram for the Mean Absolute Error is presented in Figure 76. It is not 


expected that the histogram should be normally distributed since the mean absolute error 
indicates no negative values with the theoretical goal of reaching zero for all values. The 
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graph suggests that over 90 percent of the wells have achieved the calibration goal of 1.0 
foot or less. 


Mean Absolute Error of Monitoring Wells 
Lower East Coast sub-Regional Model 
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Figure 76. Mean Absolute Error Histogram for Model Calibration Period 


The histogram for the Root Mean Square Error is presented in Figure 77. Similar 
to the Mean Absolute Error, it is not expected that the histogram should be normally 
distributed since the absolute indicates no negative values with the theoretical goal of 
reaching zero for all values. A slightly lower number of wells reach the calibration 
criteria of 1.0 foot with approximately 85 percent of the wells below the 1.0 foot target. 


Root Mean Square Error of Monitor Wells 
Lower East Coast sub-Regional Model 
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Figure 77. Root Mean Square Error Histogram for Model Calibration Period 


The final primary calibration statistic is plus or minus 1.0 foot which is that the 
simulated data should be within a band of plus or minus one foot for a minimum of 75 
percent of the time. The calibration results indicate that approximately 78 percent of the 
wells meet or exceed this condition. 
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Of the three secondary calibration statistics the objective was to achieve or exceed 
the desired result a minimum of 75 percent of the time. For RES, 99 percent of the wells 
achieve or exceed this condition. For the Min/Max and STD statistical criteria, they meet 
or exceed there condition 81 and 99 percent respectively. It should be noted that the S2 
secondary calibration statistic reports slightly better calibration results than the 54 
calibration statistic. The S4 calibration statistic is a plus or minus 1.0 foot band where 
the S2 calibration statistic keeps the plus or minus 1.0 foot band but allows the band to 
expand for those observation wells where the standard deviation exceeds 1.0 feet. The 
observation wells where this occurs are areas where fluctuation of the water table is more 
extreme. 


Individual Well Calibration 


The same observations wells with continuous data recorders were analyzed for the 
individual well calibration. Table 13 presents the statistical results of the individual well 
calibration. Figures 78 through 81 show the location of the monitoring stations. 


Table 13. Statistics for Calibration Period. 


100 
100 


0.33 
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89.63 
92.24 
94.66 
99.36 100 
| 62034 | 9936 | 7642 | 9994 | 913 | 006 | 045 | 06 ` 
| G2305 | 100 | 6735 | 100 | 3584 | -o66 | 167 | 208 | 
63259A | 100 | 9882 | 100 | 4794 | o | 104 | 117 | 
99.86 
| G3273 | 99.76 | 9054 | 9746 | 8796 | 007 | 052 oe | 
|.G3327 | 996 | 8181 | 100 | 9793 | 002 | 028 | 039 | 
| G3439 | 100 | 9125 | 999 | 9355 | aa 
| G3466 | 100 | 6708 | 100 | 6421 | -006 | 085 | 101 | 
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Well Std Dev MAE 
G347 0.41 
G3549 0.32 
G355 85.38 0.22 
0.33 
86.96 | 99.9 | 96.26 0.47 
G355 0.78 
: 1.04 
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G355 6188 | 100 | 851 | 06 | oe | 073 
0.32 
G355 0.37 
G3559 9634 | 9959 | 9969 | 009 | 0.14 | 023 
G356 92.96 0.56 
G3561 100 0.92 
G356 80.69 0.54 
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Std Dev 
94.72 
92.45 
96.84 
97.11 
97.31 
85.99 96.01 
79.9 
5.82 
83.16 
8762 | 9946 | 549 | -o9 | 111 | 14 | 
81.04 
99.29 
32.76 
6569 | 100 | 6937 | 042 | 08 | 109 | 
80.37 
89.38 
9099 | 9967 | 9641 | -0.09 | 03 | 04 | 
86.39 
96.88 100 
94.86 
89.46 
8284 | 9638 | 811 | -002 | 059 | 08 | 
81.18 
91.22 
77.24 
9984 | 9994 | 9996 | o | 008 | 012 | 
100 96.69 
80.24 
82.75 
9145 | 9916 | 916 | 002 | 048 | 066 | 
| NP67 | 999 | 8972 | 9984 | 9596 | am | 034 | 045 | 
83.34 
100 91.17 
100 
100 
6217 | 100 | 8769 | 042 | 06 | oó | 
| NPDO! | 100 | 8478 | 9989 | 8369 | ag | 06 | 076 | 
100 99.02 | 989 | 100 | ong | 011 | 014 | 
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Figure 78. Location of the Observation Network Showing Groundwater Observation Wells and 
Surface Water Gages. 
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Figure 79. Location of the Observation Network in the Northern Model Area. 
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Figure 80. Location of the Observation Network in the Central Model Area. 


Draft - 150 


LEC subRegional MODFLOW Model Documentation 


Groundwater 
Observation Well 


Surface Water 
Gage 


Active Wetland 
Boundary 





DRAFT 


G1 22562035 1226 “Y: 
G1 Y 


Lower East Coast 
subRegional Model 
"ys Fe GOY 
Prepared by : MAPU 
Date: 03/01/06 
Map Doc.: Obswells N C S 06.mxd 


Figure 81. Location of the Observation Network in the Southern Model Area. 
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In order to determine if any localized trends exist in the calibration results, the 
mean error was analyzed aerially across the model domain. The mean error was chosen 
because it would provide insight into localized areas where the model may be over- 
predicting or under-predicting heads. The location of the monitor wells and surface water 
gages used in the comparison of the mean error, and the mean error calculations are 
presented in Figure 82. Figure 83, Figure 84 and Figure 85 show the water levels for 
average (1992), wet (1995) and dry (1989) years Figures 86, 87 and 88 show the 
number of days when ponding occurs for average (1992), wet (1995) and dry (1989) 
years, in the southern portion of the model. 
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Figure 82. Mean Errors at Observation Sites for Average Daily Water Levels During Model 
Calibration Run. 
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Figure 83. Average (1992) water levels(average daily water levels) 
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Figure 84. Wet (1995) water level (average daily water levels) 
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Figure 85. Dry (1989) water level (average daily water levels) 
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Figure 86. Hydroperiod (pond 
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Figure 87. Hydroperiod (ponding) for Wet (1995) year 
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Figure 88. Hydroperiod (ponding) for Dry (1989) year. 
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Calibration and Simulation for Above Ground Impoundments 


Inflows to Water Conservation Areas 1A, 2A and 3A are a combination of water 
from Lake Okeechobee, which is an external source of water to the model, and outflows 
to the urban areas. In addition, movement of water through the various compartments is 
accomplished in a slightly different manner utilizing the redirected flow package. The 
net flow for a conservation area is equal to the Lake Okeechobee inflows minus the urban 
outflows. If the Lake Okeechobee inflows are greater than the urban area outflows the 
net volume for the Conservation area will be positive indicating that water is entering the 
WCA, if it is negative then water is being removed from the impoundment. This is done 
for each individual day of the simulation period. Figure XX compares the simulated 
versus historical net flows for Water Conservation Area 1 of Lake Okeechobee Inflows 
minus Urban outflows. Urban outflows from Water Conservation Area 1 are primarily to 
the Lake Worth Drainage District, Broward County and the City of Boca Raton. A 
similar approach is used to move water to Everglades National Park and the other water 
Conservation Areas. 


Flow from one compartment to another compartment within the WCA’s is 
accomplished utilizing the redirected flow package. The redirected flow package is 
needed to impose upper limits on the amount of water that can be within a water 
conservation area to protect the levees that surround it from failing. In addition, 
operational rules for moving water through the Water Conservation Areas to Everglades 
National Park can be simulated. Figures 89 and 90 show the simulated versus historical 
daily water levels, with a +/- 1 foot error band, for Water Conservation Area Number 1 
(Gage 1-9) and Everglades National Park in Northeast Shark River Slough (NESR2). 
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Figure 89. Stage Hydrographh for 1-9 
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Stage Hydrograph for NESRS2 
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Figure 90. Stage Hydrograph NESRS2 


Estimation of Flows to the Loxahatchee River 


Estimating flows to the Loxahatchee River is one of the most complex processes 
in the model as far as conceptualization, design formulation, and calibration. This process 
will be presented like a case study to fully understand the complexity of the natural and 
managed components of the system. The Loxahatchee River is a wild and scenic river 
located in the northeastern portion of the model in Palm Beach and Martin Counties. The 
contributing basins for the river include the Loxahatchee Slough, the C-18 west sub-basin 
and Jupiter Farms sub-basin. Flows from these three sub-basins enter the Loxahatchee 
River via the Lainhart Dam. The river then meanders northward until ultimately 
discharging to tide at the Jupiter Inlet. Figure 91 shows the main features associated 
with the Loxahatchee River. 
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Figure 91. Main Features of Loxahatchee River 
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The Loxahatchee River itself is simulated in the model as an internal boundary 
utilizing the Drain package because it is essentially brackish to marine for a significant 
portion of the river reach. Because of its direct tidal connection, and the need to keep 
freshwater flowing to the river so it doesn't complete become marine, flows over 
Lainhart Dam are consider crucial in maintaining the integrity of the river. A simplified 
approach was utilized to simulate flows over Lainhart Dam by determining flows from 
the three sub-basins which contribute to the flow. This was accomplished utilizing the 
wetland and diversion packages. The C-18 west and the Loxahatchee Slough sub-basins 
are already simulated utilizing the wetlands package. Therefore, there is no 
preprocessing of rainfall for these areas so recharge in the model is equal to actual 
rainfall and et is equal to the potential ET for wetlands. The Jupiter Farms sub-basin 1s a 
low density urban area so the ET and recharge for this area was modified and treated as 
wetlands during the recharge preprocessing program so full rainfall would be accounted 
for in the model in the Jupiter Farms area. The various canals that drain the three sub- 
basins were simulated using the diversions package instead of the conventional river or 
drain package. The diversion package allows the user to specify a specific volume of 
water to move through that reach assuming the upstream and downstream reaches meet 
the operational criteria specified. Taking the C-18 West weir as an example, water levels 
upstream would need to be above 17.6 feet NGVD, which is the weir crest, and the down 
stream reach must be lower than 17.6 for water to move from the C-18 West sub-basin to 
the Loxahatchee River sub-basin at pre-specified volume of water. The C-18 West sub- 
basin, because it has several different control elevations, higher than 17.6 feet NGVD, 
that actually discharge into the C-18 West canal was further sub-divided into 7 smaller 
sub-basins. The order of operations is as follows. When the water level rose above the 
control elevations, generally between 19.0 feet and 24.0 feet NGVD, for the upper sub- 
basins of the C-18 West sub-basin, water would move from these basins into the C-18 
West canal reach. If the C-18 West Canal reach rose above 17.6 ft NGVD, water would 
move into the C-18 Canal. If the water in the C-18 canal was above 13.0 feet NGVD, 
water would move into the C-14 Canal where it would flow over Lainhart Dam if the C- 
14 canal was above 11.0 feet. Jupiter Farms flows directly into the C-14 Canal and 
ultimately over Lainhart Dam. This simplified approach allows the user to move water 
from one section of the model to another utilizing the overland flow capability of the 
wetlands package and the operational components of the diversions package essentially 
creating a “cascade” type approach. 


Figures 92, 93, and 94 illustrate the simulated versus historical flow for the three 
control structures regulating flow to the Loxahatchee River. 
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Figure 92. C-18 West Weir Mean Monthly Flows. 
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Figure 93. G-92 Mean Monthly Flows. 
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Figure 94. Lainhart Dam Mean Monthly Flows. 


Miami-Dade Hialeah/Preston Wellfield Simulation 


The Miami-Dade Water and Sewer Authority’s Hialeah/Preston wellfield is one of the 
largest in the State of Florida. Calibration at this wellfield illustrates (Figure 95) some of 
the inherit problems with the model resolution, data collection and the daily stress period. 
The Hialeah/Preston water treatment plant supplies most of the people in northern 
Miami-Dade County. It receives water from five different wellfields including the 
Northwest wellfield (15 production wells), Upper Miami Springs wellfield (12 
production wells), the Lower Miami Springs wellfield (8 production wells), the Hialeah 
wellfield (3 production wells) and the Preston wellfield (7 production wells). Total water 
treated at the plant is approximately 160 million gallons per day (MGD). The Hialeah, 
Preston Upper Miami Springs and Lower Miami Springs wellfields are all in close 
proximity to each other and can be considered as a single wellfield from a regional 


perspective. 
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Figure 95. Location of Miami Dade Water and Sewer Authority Wellfields 


Historically, municipal utilities were required to submit total monthly raw water 
withdrawals for each individual treatment plant. In this case, five wellfields are 
supplying one water treatment plant. Daily data for each wellfield or each production 
well, is not available. In this case, the water treatment plant is producing approximately 
160 MGD but the distribution between wellfields or individual production wells is not 
known which is problematic. To further complicate things, it appears that a major shift in 
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wellfield withdrawals occurred midway through the calibration process as shown on 
Figure 96. 
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Figure 96. Water Levels in Miami-Dade County Wellfields. 


Figure 97 shows the location of the monitoring wells in the vicinity of the 
Hialeah/Preston and Miami Springs wellfields in Miami-Dade Counties. Observation 
wells G-3327, G-3328, S-68 and G-3467 are located in close proximity to the wellfields. 
These wells are reasonable well calibrated as shown in Figure 98. However, wells 
located within the wellfield itself including well S68 which are not as well calibrated as 
shown in Figure 99. One of the main reasons for this is the lack of wellfield and 
individual well data for the calibration period. The lack of actual daily pumpage data 
from the individual wells results in a poor calibration closer to the wellfield. This 
illustrates the point that for larger wellfields, additional data collection, potential 
rediscritation and partial recalibration of the model is needed before it can be utilized at 
the local scale. 
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Figure 97. Observation Wells Near X Wellfield. 
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Figure 98. Stage Hydrograph for G3327 
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Figure 99. Stage Hydrograph for S68 
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Volumetric Budget 


The volumetric budget for the model area shows the magnitude of the flow 
components in the active model domain. An analysis of the volumetric budget helps in 
developing better water management practices and plans. The volumetric budget is also a 
good indicator of whether the model results are reasonable. A transient state volumetric 
budget for the LECsR model at the end of the calibration simulation is shown in Table 
14. 


Table 14. LECsR Model Area Volumetric Budget at the End of the Simulation 


Cumulative Volume (ft3) % 10 % Out 
Inflow Outflow 
Storage 5.0128E+12 5.0541E+12 28.12 
28.35 
Constant 0.0 0.0 0.00 0.00 
Head 
Wells 1.0000E+9 5.9540E+11 0.01 
3.34 
Drains 0.0 5.4075E+11 0.00 3.03 
River 3.3327E+12 4.7444E+12 18.70 25.10 
Leakage 
ET 0.0 4.4268E+12 0.00 
24.83 
Head 6.1050E+11 1.5291E+12 3.42 8.58 
Dependent 
Boundary 
Recharge 6.9891E+12 0.0 39.21 
0.00 
Redirected 8.1845E+11 8.1845E+11 4.59 4.59 
Flow 
Source/Sink 1.0613E+12 3.8759E+11 5.95 
2.17 
Totals 1.7826E+13 1.7826E+13 100.00% 100.00% 


The largest inflow to the model is rainfall accounting for nearly 40 percent of 
inflows. Storage changes and river leakage represent the next two largest inflows 
accounting for 28 percent and 19 percent respectively. Source and sinks, redirected flow, 
head dependent boundaries and well reuse complete the remainder of inflows total 
approximately 14 percent. The primary outflows from the model are relatively evenly 
distributed between storage changes, canal drainage and evapotranspiration. These three 
account for over 80 percent of the outflows for the model. The remaining outflows are 
the head dependent boundaries, redirected flow, sources and sinks, and wells which are 
only 3 percent of the entire water budget. 


VERIFICATION 


Calibration is carried out for a period presenting different hydrological conditions. 
Verification serves to test whether the set of parameters selected during the calibrations 
are in fact suitable to represent a different period. In the LECsR Model project the period 
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1986-1999 has been selected as the calibration period and the period 1999-2000 is the 
verification period. In principle, the model calibration results should remain relatively 
unchanged for the two periods. 


The overall model statistics for the verification run do not change drastically from 
the calibration run as shown in Table 15.. The global standard deviation was less for the 
verification run than for the calibration run although the mean absolute error and the root 
mean squared error are slightly improved. 


Table 15. Calibration Results for the Entire Model Domain (195 observation wells) 
for the verification period (1995-2000) 


Verifi Calibr 
Statisti catio ation 
cal n Rang 
criteria Rang e 
e 
RES o 100 % 
Global 99% 
STD 80 % 
Global 71% 
MIN/M 99 % 
AX 96 % 
Global 
+/- 1.0 87 % 
Global 88 % 
ME 0.00 ft 
Global 0.11 ft 
MAE 0.55 ft 
Global 0.54 ft 
RMSE 0.70 ft 
Global 0.67 ft 


In general the wells that were consider calibrated after the calibration run 
remained calibrated for the verification run. No noticeable improvement in the wells 
consider not calibrated was observed in the verification scenario. Figure 100 shows the 
difference of the mean absolute error between the calibrated and verified simulation 
periods for each well. Approximately 95 percent of the wells are within a change ot +/- 
0.5 feet with approximately 85 percent of the wells changing less than 0.25 feet. 
Considering the verification period was one tenth in length as the calibration period, the 
results suggest that the model appears to be robust in its ability to simulate conditions 
adequately. 
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Figure 100. Difference Between Calibration and Verification Runs for the Mean Absolute 
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Figure 101. Simulated water levels and mean error for verification run. 
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SENSITIVITY ANALYSIS 


After a groundwater flow model has been calibrated, a final sensitivity analysis is 
performed. Sensitivity analysis is defined as “a quantitative evaluation of the impact of 
variability or uncertainty in model inputs on the degree of calibration of a model and on 
1ts results or conclusions” (ASTM, 2002). Sensitivity runs were conducted to test the 
magnitude of the model’s response on the range of simulated outputs to changes in 
stresses, aquifer parameters, and surface-groundwater interaction. 


The purpose of the sensitivity analysis is to vary input parameters within an 
acceptable range. Acceptable ranges of variation for input parameters were decided 
based on model output response and on expected error on input parameters. 


Parameters and Methodology 


The sensitivity of the output variables to variations in input parameters is 
estimated by the traditional approach of varying one parameter at a time by a constant 
factor. For each parameter, two model runs were completed by increasing and decreasing 
the value. The simulation period for the sensitivity analysis was within the calibration 
period, January 1986 to December 1995, which has a wide range of hydrologic 
conditions ranging from very dry to average to very wet hydrologic periods. Parameters 
tested included recharge, ET maximum rate, ET surface, ET extinction depth, horizontal 
hydraulic conductivity (HK) for all layers, vertical conductance (VCONT), specific yield 
(SY), river and drain conductance, general head boundary (GHB) conductance. The 
range of values tested for each parameter is shown in Table 17. 


Several parameters specific to the Wetlands package were also tested, including 
Kadlec coefficient, specific yield of the surface water body (SYWTL), specific yield of 
the muck/peat (SYMUC), horizontal hydraulic conductivity of the muck/aquifer layer 
underneath the water body (HYMUC) and bottom elevation of surface water (ZBOTT). 
Kadlec coefficient in wetland areas represents vegetative resistance to flow or drag forces 
on moving water as it flows through wetland vegetation. The soil/aquifer layer 
underneath the water body is simulated as part of the top layer (e.g., surface water body), 
as an independent confined/unconfined layer, or as part of the aquifer underneath. 
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Table 17. Sensitivity Analysis Parameters for LECsR Transient Simulation 







x 0.7 
Aquifer parameters SO 
x 0.2 












Surface-ground water interaction 
Kadlec coefficient x1.5 











ZBOTT 
GHB conductance 


x 0.01 
x 0.1 
x 0.1 








(1) Change from calibrated value: x or parameter multiplied by a factor of; 
+/- or parameter raised/reduced by 


Sensitivity Results 


For graphical data analysis for sensitivity model runs, dot diagram and frequency 
distribution table are included. These plots are built with residual values computed as the 
difference between simulated hydraulic heads for test run and calibrated model. The dot 
diagram plots show the distribution of the average, minimum, maximum and standard 
deviation of the head difference for analyzed parameters, and for each monitoring well 
throughout the simulation. They summarize visually the spatial dispersivity of the 
sensitivity for analyzed parameters and identify outliers. The boxplots show the 
distribution of the average, minimum, maximum and standard deviation of the head 
difference for the analyed parameters, and for the entire model active zone. 


Sensitivity results are grouped by urban and wetland areas. The histogram 
examines the frequency distribution of the mean daily residuals for observation wells, 
and for the entire model active cells with class marks of 0.5 ft. The dot diagram plots 
show the distribution of the average, minimum, maximum and standard deviation of the 
head difference for analyzed parameters, and for each monitoring well throughout the 
simulation. The boxplots summarize the median and the quartiles of the analized data. 
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Sensitivity Results for Observation Wells 


Frequency distribution of mean heads residuals (simulation tested minus 
calibration) based on 197 wells is presented in table 18. Tables 19 to 26 shows the effect 
of average, minimum, maximum and standard deviation head difference for each urban 
well and wetland in several parameters. 
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Table 18. Frequency Distribution of Mean Heads Residuals Based on 193 Wells 
CLASS OBSERVATION WELLS FREQUENCY 
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Table 19. Effect of Average Head Difference for Each Wetland Well in Several 
Parameters. 
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Table 20. Effect of Average Head Difference for Each Urban Well in Several 
Parameters. 
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Table 21. Effect of Minimum Head Difference for Each Urban Well in Several 
Parameters. 
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Table 22. 
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Table 23. Effect of Maximum Head Difference for Each Urban Well in Several 
Parameters. 
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Table 24. Effect of Maximum Head Difference for Wetland Well in Several 
Parameters. 
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Table 25. Effect of Standard Deviation Head Difference for Each Urban Well in 
Several Parameters. 
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Table 26. 
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Effect of Average Head Difference for Each Wetland Well in Several 
Parameters. 
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Stresses and Related Variables 


ET maximum rate 


The model is sensitive to increased ET for the urban and wetland environments 
(the range of mean residuals in -3 to 1.5 feet). The Wetlands areas are more sensitive to 
increases in the maximum ET rate. Gage 2B-Y in the wetland area is especially sensitive 
to variations in the rate with the mean residual values being reduced by up to 4 feet. 
Therefore the simulated water level values are lower than the calibration run when the 
rate is increased. Wetland wells are much more sensitive to changes in the maximum 
rate than non-wetland cells. Urban wells changed and average of between -1.5 to 0.5 ft 
from calibrated levels. 


ET surface elevation 


Decreasing ET Surface lowered the average mean residual by 0 to 1.5 ft, in 
wetlands. When the ET surface was increased the average residuals increased by an 
average of 0 to 1.0 ft. 


ET extinction depth 


Increasing extinction depth impacts the wetland wells more then the urban wells 
(an average maximum residual of -5 feet compared to -1). Decreasing the extinction 
depth, increased the mean residual by -0.5 foot for both urban and wetland wells As 
expected, when the extinction depth is deeper the water levels tend to drop, when it is 
shallower the water table increases. 


Recharge 


The model is sensitive to changes in recharge. When the recharge was increased 
the mean residual range for both Urban and Wetland area was between 0 to 1.5 ft. 
Frequency distribution table shows that the majority of the wells residual increased by 
0.25 feet. When the recharge values decreased the impact was greater in the wetland 
wells with the mean residual ranging from -1 to -2.0 ft. (The negative residual indicate a 
decrease in the water levels in the wetlands when recharge is decreased). The model 
response to recharge in non wetland areas is moderate due to the influence of canals 
which tend to remove the water quickly due to the high hydraulic conductivity of the 
aquifer. Based on frequency distribution table, water levels tended to decline in many 
wells by reducing recharge and a clear skew exist towards reduced water levels by up to 
2.0 foot. Only one well in urban areas was observed to increase significantly over 
calibrated values. This well is located near a drain. It is speculated that the drain is 
pulling the aquifer levels down under wet conditions but when reduces rainfall is applied 
the drain is off so it tends not to pull the water levels down. When recharge 1s increased, 
the exact opposite happens. Most wells show an increase in water levels. Frequency 
distribution table shows the skew towards increased water levels. 
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Aquifer Parameters 
Horizontal hydraulic conductivity in layer 1 (aquifer 


The model appears to be not very sensitive to changes in the hydraulic 
conductivity of layer 1. This is a result of several factors. In the wetland areas, changes 
to the Block Centered Flow (BCF) package of modflow has no effect on the wetland 
areas because the wetland package over-writes the hydraulic conductivity values 
originally read in from the BCF package. In general, layer | has a significantly lower 
hydraulic conductivity and thickness than the underlying layer 2 which is the primary 
production zone for most of the model domain. 


Horizontal hydraulic conductivity in layers 2 and 3 


Layers 2 and 3 were modified together rather than separately. The model showed 
relatively greater sensitivity to these parameter in both wetland and non-wetland wells. 
When the hydraulic conductivity was increased the Urban areas responded greatly with 
the range of mean residuals is -8 to 3 ft (with most of the wells falling 10 the -1 to 1.5 foot 
range). In the wetland areas, the range is -4 to 4 ft, with most of the wells in the -0.5 to 
0.5 foot range. When the hydraulic conductivity was decreased the mean residual in the 
wetland areas ranged from -3 to 2 ft. The residual range in the urban areas varies from -1 
to 1 ft. For the majority of the wells, the mean residual decreased by 0.25 ft when the 
hydraulic conductivity was raised and went up 0.25 ft when it was lowered. The 
histogram shows a less skew distribution when the hydraulic conductivity was increased. 
When the values were decreased the pattern is skewed. Most of the wells (115) residuals 
increased by 0.25 a foot but 19 well residual decreased (water levels decreased) when the 
hydraulic conductivities were lowered. The model sensitivity to changes in this 
parameter was closely linked to recharge. 


Vertical Conductance for Layer 1 


As a whole, the model was fairly insensitive to the changes in vertical 
conductance (VCONT) values between layers 1 and 2. The model showed much greater 
differences in simulated heads in urban areas than wetland areas. Decreases in VCONT 
generally produced higher water levels in layer 1 wells, indicating a restriction in 
downward movement of groundwater, primary recharge from rainfall and canal leakance, 
from layer 1 to layer 2. An interesting exception occurs at well G853 in coastal Broward 
County. Decreasing the VCONT produced lower water levels than the calibration did, 
indicating restriction of upward flow from layer 2. 


Vertical Conductance for Layer 2 


Changes made to the VCONT values between layer 2 to layer 3 had very little 
overall effect on calibrated water levels in the observation points. 
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Specific Yield for Layer 1 


The model is insensitive to changes in specific yield in the urban areas. Most of 
the wetland wells were not sensitive to changes in this parameter and mildly sensitive 
(changes of tenths of a foot) in the urban area. It should be noted that the specific yield 
for layer 1 is replaced by the specific yield specified in the wetlands package for the 
wetland areas. Therefore, in these sensitivity runs, the specific yield in the wetland areas 
is not modified. This increased sensitivity to specific yield during high recharge periods 
was also observed in most of the non-wetland wells. The model was far more sensitive to 
reducing the specific yield than increasing it. The minimum and maximum distribution 
dot diagram for decreasing the specific yield shows daily residuals for various wells 
indicating a change up to 4 ft at well PB-565. 


Surface-groundwater Interaction 


Kadlec coefficient 


Non-wetland wells were insensitive to changes in this parameter. Wetland wells 
were only slightly sensitive to Kadlec coefficient changes. Sensitivity fluctuations in 
individual wells could not be definitively explained, and could be the result of a complex 
interaction with the operational packages (diversions and redirected flow). In well 2B-Y, 
a high sensitivity to decrease Kadlec is being affected by both recharge and by 
operational features, such as outflows from the reinjected drainflow package. Well 2B-Y 
is located in Water Conservation Area 2B which is the smallest of the WCA's and is 
bordered on two sides with urban seepage collection systems resulting in high gradients. 


Specific yield in wetlands 


The model responded with increased sensitivity to decreased specific yield values. 
The entire model was insensitive to the 1.0 value for specific yield of the surface water 
body. 


Specific yield muck/aquifer underneath water body 


Non-wetland wells were insensitive to changes in the muck specific yield. 
Wetland wells were more sensitive to decreasing the value than increasing, with 
differences in water levels from calibrated levels up to 2 ft. Sensitivity fluctuations in 
individual wells could not be definitively explained, and could be the result of a complex 
interaction of operational and other factors. 


Horizontal hydraulic conductivity of the muck/aquifer underneath 
water body 


The model's response to changes in this parameter varied. In general, non- 
wetland areas were insensitive to that parameter, with a few exceptions (e.g PB-1642) 
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when the parameter is lowered. These wells showed increased sensitivity (lower water 
levels) that closely correlated with periods of high recharge. The non-wetland wells were 
largely insensitive to the tenfold increase in muck conductivity. 


The wetland wells were sensitive to changes in this parameter. The wetland wells 
located in the Everglades showed very similar sensitivity in both the decreased run and 
increase run, in both magnitude and direction of change from calibrated values. This is 
because the calibrated muck conductivity values are extremely high in those areas. 
Wetland wells did show a pattern of increased sensitivity during different times, but it 
was not as closely related to recharge as was seen in the non-wetland wells. Increase 
sensitivity to muck conductivity in modeled wetlands is due to water surface fluctuations 
(non-pounded areas) probably also influenced by water management operations. 


General head boundary conductance 


Sensitivity to this parameter was negligible except for areas very close to 
boundaries. The test runs showed more sensitivity when boundary conductance was 
decreased than when it was increased. 


River conductance 


The model is sensitive to changes in the river and drain conductance, with the 
Urban areas being more sensitive to increased in conductance. The mean residual range 
in the Urban area is -0.5 to 1.0 ft. The majority of mean residual values by increasing 
value were between -0.5 and 0.5 feet. The model was sensitive to reduction in river 
conductance but began to experience convergence issues in some areas when the 
conductance was increased. This is not unexpected, considering the abundance of canals 
and drains in the LEC and the high hydraulic conductivity of subsurface materials. 
When the canal conductance is decrease by a factor of 0.1, the frequency distribution 
table shows a slight bias to higher water levels which would be expected. The numerous 
canals that have been constructed in south Florida where built primarily to drain the land. 
So by limiting the canals ability to remove water, it would tend to increase water levels 
when conductance values where lowered. 


Many of the observation wells showed increased sensitivity to changes in river 
conductance during periods of high recharge. Some wells would respond with higher 
water levels during high recharge when the conductances were lowered, because the 
nearby canals were less able to get rid of the water. Conversely, those wells had lower 
water levels than the calibrated run when conductances were high, because the nearby 
canals lost more water through the open connection to the aquifer. In other wells such as 
JD12 or 2B-Y, the sensitivity to river conductance was more complex, being affected by 
both recharge and by operational features, such as outflows from the reinjected drainflow 
package. Well 2B-Y is also located in Water Conservation Area 2B which is the smallest 
of the WCA's and is bordered on two sides with urban seepage collection systems 
resulting in high gradients. It should be noted that instability was observed in the major 
wetland systems when the conductance was increased. Care should be taken when 
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increasing the conductance values for the perimeter canal systems that surround these 
systems. 
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CHAPTER 5 


Conclusions and Future Improvements 


MODEL LIMITATIONS 


A model is any device that represents an approximation of a field situation 
(Anderson and Woessner 1992). The model described in this report is a numerical 
groundwater flow model with coupled surface/groundwater flow in the wetlands area 
utilizing the well known computer code MODFLOW to approximate, on a subregional 
scale, the groundwater flow system within the Surficial aquifer system in southeast 
Florida. 


The Lower East Coast Subregional flow model developed in this report is based 
upon a simplified representation of the complex heterogeneous groundwater flow system. 
The governing equation utilized to solve the flow system is the continuity equation 
derived from the principle of mass balance incorporating Darcy’s Law. This equation 
assumes that the flow in laminar and does not reach turbulent conditions which is general 
true throughout the study area, for the groundwater system, with the possible exception of 
flow near some major production wells in the Biscayne aquifer. The horizontal and 
vertical discretization of the model assumes that the hydraulic properties are the same 
throughout the cell which a simplification of the aquifer properties considering the size 
and thickness of each grid cell. 


The spatial distribution of evapotranspiration and recharge is dependent upon 11 
and 26 stations respectively. Considering the sporadic nature of rainfall in south Florida 
some error may be introduced because the data sets do not necessary represent actual 
rainfall patterns. Recent advancements in meteorology, including the widespread use of 
Doppler radar should help in better refining rainfall patterns over a large area. However, 
this type of data is presently not available for simulations of conditions that occurred 20 
years ago. In addition the historical data from the monitoring sites has data gaps that 
must be filled from adjacent stations which introduce additional error to the model. 


The lack of water level data, especially in Martin and Palm Beach Counties 
introduces uncertainty into the model in certain areas. There is a complete lack of any 
available monitoring data throughout the Everglades Agricultural Area. There is only 
sporadic data available for northern Palm Beach and Martin Counties. Central and 
southeastern Palm Beach Counties are not as problematic due to the highly managed 
secondary canal systems. The actual water level data itself was in generally good 
condition with no significant problems detected. Some of the flow data from the 
structures appear to have some issues which may introduce some additional errors. 
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The simulation of flow through the C-18 canal utilizing the diversion package is a 
simplified method and can not be used for canal sizing or routing capabilities. This 
method was utilized to see 1f long term simulations could reasonably predict structure 
flows in a computationally efficient manner. Similarly, the redirected flow package and 
the diversion package utilized to move water through the Water Conservation Areas 
should include additional features to improve the routing of water through the system. 


The topography utilized in the model development was in a constant state of 
updating, modifying, adding, changing and general disarray. Errors associated with the 
topography primarily impact to the areas covered by the wetland systems in the model. 


The model is designed to be handle regional and basin level projects. Projects 
utilizing the model for local areas such as a wellfield should consider a partial 
recalibration of the model at the site specific location with a detailed collection of all 
historical data available, especially wellfield and individual well daily pumpage rates. 
The size of the grid should also be taken into consideration when evaluating local 
projects. 


In spite of the limitations, the model can provide an understanding of the 
movement of water in the study area. The model proved quite robust considering the 
long term daily calibration and verification periods in its overall performance during the 
simulation period. The calibration period included both a 1 in 100 year drought and a 1 
in 100 year wet events and the model did not react adversely to these extreme events. 
The LECsR model conceptualization and discretization was designed at a subregional or 
basin level scale. The special variability of the input parameters is also best described at 
a similar scale. Therefore, the model should be used for regional to subregional or basin 
level projects and interpretation of the results should also be at that scale. The hydraulic 
parameters utilized in the model can be used as initial estimates for more localized 
models. In addition, the model does provide a reasonable estimate of drawdowns 
associated with wellfield withdrawals and the ground/surface water interactions within 
wetland systems. 


CONCLUSIONS 


The Lower East Coast of Florida includes all or portions of southern Martin, Palm 
Beach, Broward, Miami-Dade, mainland Monroe and eastern Collier and Hendry 
Counties. The Surficial aquifer system 19 the principle aquifer in the region and includes 
the extremely prolific Biscayne aquifer with yields in excess of 7,000,000 gallons per day 
per foot of drawdown. This area also includes approximately 50 percent wetland systems 
including the world renowned Everglades system and the Everglades National Park. 
Also included in the region is the vast sugar cane fields of the Everglades Agricultural 
Area and the sprawling urban areas along the coast. A three dimensional numerical 
groundwater flow model with coupled surface-groundwater in the wetland areas was 
developed using the USGS code MODFLOW with additional packages to simulate 
conditions unique to south Florida. The model was developed by combining and 
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updating previous county level models which simulated coastal conditions and expanding 
the overall extent of these models to include the majority of the lower east coast of 
Florida. The overall active model domain is approximately 7,500 square miles. 


The sediments composing the Surficial aquifer system and the Biscayne aquifer 
are primary marine carbonate and clastic sediments of Pleistocene to Miocene age. The 
model was discretized into three layers utilizing Chronostratigraphic correlation of the 
sub-aerial exposure surfaces within these sediments resulting from the periodic 
submergence and emergence of the Floridan peninsula from fluctuations in sea level 
stands. A grid spacing of 704 foot by 704 foot was chosen to match the SFWMD 
regional models which can be used to provide internal boundary conditions for future 
simulations. The model was calibrated over an extend period of 14 years from January 
1986 through September of 1999 with daily time steps and stress periods. The model 
was primarily calibrated to observed heads with an overall mean error of 0.0 feet, a mean 
absolute error of 0.54 feet and root mean square error of 0.72 feet. The calibration results 
indicate a reasonable match between observed and measured water levels in most areas of 
the model domain. The model was verified from September, 1999 through December, 
2000 and produced similar results to the calibration period. The budget from the model 
indicates that recharge to the aquifer is the primary inflow to the system with 
evapotranspiration and canal drainage as the primary outflows. 


A sensitivity analysis was conducted in which 36 individual variables were 
analyzed where the values were increased and decreased within a predefined range of 
known values. The results were analyzed utilizing the mean head residuals, the 
maximum head residuals, the minimum head residuals and the standard deviation of the 
head residuals. The results were also evaluated separately for the developed and wetland 
areas. The results for the developed areas indicate that increases and decreases in the 
hydraulic conductivity of layers two and three and decreases in the river conductance 
were the most sensitive. In addition to these three parameters, noticeable individual wells 
showed biased sensitivity to changes to specific yield in layer one and the amount of 
recharge applied to the urban areas. The wetland systems were significantly more 
sensitive to changes in parameter values with noticeable changes to wetland stages when 
changes were made to all ET parameters, recharge, changes to the Kadlec number, the 
hydraulic conductivity of layer one underlying the wetland systems, changes to the 
hydraulic conductivity of layers two and three, the general head boundary and changes to 
the river conductance values. 


Some conclusions and recommendations from the model development effort are 
as follows: 


e Although the model covers a large area, the grid resolution is such that 
it can be used for both a basin level to regional analysis. 


e Use of the model at a local scale requires additional data collection and 
partial recalibration of the model prior to its use at a wellfield level. 
Depending upon the project, a finer grid size may also be necessary. 
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* The creation of a single model, as opposed to six county level models 
is easier to maintain, easier to simulate single or multiple projects, and 
ensures a consistent methodology for developing input data sets. 


* The model reasonably simulates the groundwater/surface water 
interactions in the wetland systems. Some potential improvements 
include the addition of a dual porosity component to increase overland 
flow in a preferential pathway and the modification of the wetlands 
code to include Manning’s roughness coefficients. 


* An unsaturated zone package similar to that developed by FAU (2001) 
and Geotrans (2004) should be included in future model refinement to 
eliminate the cumbersome pre-processing for the evapotranspiration 
and recharge packages. 


* The model requires extensive amounts of storage when ran for 15 — 40 
years on a daily stress period. Storage space needs to be made 
available prior to running any long term simulations. 


° The sensitivity of the model should include all active cells within the 
model domain and not just the selected observation nodes. This 
should significantly increase the sample size to provide a more 
statistically robust analysis of model sensitivity. 


* The model utilized saw grass as the reference crop in developing the 
potential evapotranspiration to be consistent with the SFWMD 
regional model. This approach is not consistent with the Water Use 
Permitting section of the District which utilizes grass as the reference 
crop. The potential evapotranspiration should be recalculated with 
grass as the reference crop. 


* The topography for the model comes from a myriad of different 
sources and methods of data collection. A uniform standard for 
collection and development of the topography needs to be 
implemented. 


* The model was developed utilizing MODFLOW 96. The various 
packages should be rewritten in Fortran 90 for inclusion in 
MODFLOW 2000/SEAWAT. The use of SEAWAT would allow for 
the simulation of the saline interface and improve calibration in the 
coastal observation wells. Several studies have been undertaken 
utilizing SEAWAT in Miami-Dade and Broward counties which could 
help in implementing the density dependent conditions. 


* The use of the diversions package for simulating a simplified structure 
flow routine should be implemented elsewhere in the model. This 
should help provide a better mass balance of the overall water budget 
of the study area. 


e Develop an artificial “average” and “1 in 10 year drought" rainfall and 
evapotranspiration data sets to assist Water Use Permitting in 
regulatory issues. 
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° Conditional stochastic simulation such as sequential Gaussian 
simulations should be investigated for the distribution of point data, 
such as the hydraulic conductivity, to help eliminate the “bulls eye” 
effect of standard inverse distance weighing technique utilized in the 
development of the input data sets. 


Implementing several of the above recommendations could improve the models 
performance in several areas. However, the model in its existing condition does 
provide a good understanding of the flow system and can be used to evaluate 
water management scenarios. 


References Cited 


Abtew, W., R.S. Huebner and S. Sunderland. 2002. Part 1: Hydrological Analysis of the 
2000 — 2001 Drought in South Florida. South Florida Water Management District, 
West Palm Beach, FL. 


Adams, K. 1992. A Three-Dimensional Finite Difference Flow Model of the Surficial 
Aquifer in Martin County, FL. South Florida Water Management District Technical 
Publication 92-02, DRE-310. 


Agassiz, A. 1895. Note on the Florida Reef. American Journal of Science, XLIX, 154- 
155. 


Agassiz, A. 1896. The elevated reef of Florida. Bulletin of the Museum of Comparative 
Zoology, XXVII, 26-62. 


Agassiz, L. 1852. Florida reefs, keys and coast. U.S. Coast Survey Report for 1851, 
appendix no. 10, pp. 145-160. 


Anderson, M.P. and W.M. Woessner. 1992. Applied Groundwater Modeling: Simulation 
of Flow and Advective Transport. Academic Press, San Diego, CA. 


Anderson, P.F., White, H.O., Mercer, J.W. and P.S. Huyakon. 1986. Numerical 
Modeling of Groundwater Flow and Saltwater Transport in the Vicinity of 
Hallandale, Floridt. Geotrans, Herndon, VA 


Abtew, W., et al. 2003. Evapotranspiration Estimation for South Florida. TechnicalPaper 
#407, SFWMD, West Palm Beach, FL. 


Boulton, N.S. 1963. Analysis of Data from Non-equilibrium Pumping Tests Allowing for 
Delayed Yield from Storage. Proc. Inst. Civil Engineers, 26:469-482. 


Draft 199 


DRAFT LEC subRegional MODFLOW Model Documentation 


Brion, L.M., E.R. Santee, C.J. Neidrauer, P.J. Trimble and N.C. Krishnan. 1999. A 
Primer to the South Florida Water Management Model (Version 3.5). Hydrologic 
Systems Modeling Division, SFWMD, West Palm Beach, FL, 233 pp. (get new 
version) 


Broecker, W.S. and D.L. Thurber. 1965. Uranium-series dating of corals and oolites 
from Bahaman and Florida Key limestones. Science, 149:58-60. 


Brooks, H.K. 1968. The Plio-Pleistocene of Florida, With Special Reference to the Strata 
Outcropping on the Caloosahatchee River. In Late Cenozoic Stratigraphy of Southern 
Florida — A Reappraisal, ed. R.D. Perkins. [publisher?] Miami, FL, pp. 3-42. 


Brooks 1974. Cited in Chapter 2 


Bush, J. and N. Gutierrez. 2004. Memorandum of Agreement Regarding Acceleration of 
the Comprehensive Everglades Restoration Plan. The Executive Office of the Governor 
and the South Florida Water Management District, Boynton Beach, FL. 


Butler, D. and D. Padgett. 1995. 4 Three-Dimensional Finite-Difference Ground Water 
Flow Model of the Surficial Aquifer System in St. Lucie County, Florida. South 
Florida Water Management District Technical Publication, DRE-326. SFWMD, West 
Palm Beach, FL. 


Causaras, C.R. 1985. Geology of the Surficial Aquifer System, Broward County, Florida. 
U.S. Geological Survey, Water-Resources Investigations Report 84-4068, Branch of 
Information Services, USGS, Federal Center, Denver, CO. 167 pp., 2 sheets. 


Causaras, C.R. 1987. Geology of the Surficial Aquifer System, Dade County, Florida. 
U.S. Geological Survey, Water-Resources Investigations Report 86-4126, Branch of 
Information Services, USGS, Federal Center, Denver, CO. 240 pp., 3 sheets. 


Chen, C. 1976. Flow resistance in broad shallow grassed channels. Journal of 
Hydrology, American Society of Civil Engineers, 102(3):307-322. 


Chescheir, G.M., R.W. Skaggs, J.W. Gilliam and R.G. Broadhead. 1987. The Hydrology 
of Wetland Buffer Areas for Pumped Agricultural Drainage Water. In The Ecology 
and Management of Wetlands, ed. D.D. Hook. Croom Helm, Beckenham, England, 
pp. 260—274. 


Chin, D.A. 1990. A Method to Estimate Canal Leakage to the Biscayne Aquifer, Dade 
County, Florida. U.S. Geological Survey Water-Resources Investigations Report 90- 
4135. Books and Open-File Reports Section, USGS, Federal Center, Denver, CO, 32 


Pp. 


Draft - 200 


LEC subRegional MODFLOW Model Documentation DRAFT 


Chin, D.A. 1996. Research Plan for the Study of Melaleuca Evapotranspiration. 
Technical Report No. CEN-96-1. [Publisher] University of Miami, Coral Gables, FL, 
133 pp. 


Cooke, C.W. and S. Mossum. 1929. Geology of Florida. In Florida Geological Survey 
20th Annual Report, ed. H. Gunter. Pp. 29-228. 


Cooke, C.W. 1945. Geology of Florida. Florida Geological Survey Bulletin No. 29. 


Cooke, C.W. 1952. Sedimentary deposits of Prince Georges County (Maryland) and the 
District of Columbia. Maryland Department Geology, Mines, and Water Resources 
Bulletin, 10:1-53. 


Cooper, H.H., Jr. and C.E. Jacob. A Generalized Graphical Method for Evaluating 
Formation Constants and Summarizing Well-Field History. American Geophysical 
Union Transactions, 27(4): 526-534. 


Cooper, R.M. 1987. An Atlas of Broward County Surface Water Management. South 
Florida Water Management District Technical Memorandum. SFWMD, West Palm 
Beach, FL. 


*Cooper, R. and J. Lane. 1987. An Atlas of Miami-Dade County Surface Water 
Management. South Florida Water Management District Technical Memorandum. 
SFWMD, West Palm Beach, FL. 


Cooper, R.M. and C.J. Neidrauer. 1989. A Two-Year Field Test of the Rainfall Plan: A 
Management Plan for Water Deliveries to Everglades National Park. Technical 
Publication 89-3, Resource Planning Department, South Florida Water Management 
District, West Palm Beach, FL, 98 pp. 


Cooper, R.M. and J. Roy. 1991. An Atlas of Surface Water Management Basins in the 
Everglades: The Water Conservation Areas and Everglades National Park. Water 
Resource Evaluation Department, South Florida Water Management District, West 
Palm Beach, FL, 84 pp. 


Council, G.W. 1998. A Lake Package for MODFLOW. Proceedings for MODFLOW 798, 
Colorado School of Mines, October 4-8, 1998, Volume II, pp. 675-682. 


Council, G.W. 1999. A Lake Package for MODFLOW (LAK2): Documentation and 
User's Manual, Version 2.2. HSI Geotrans, Roswell, GA, 132 pp. 


Dall, W.H. 1897. Notes on the Geology of Florida. American Journal of Science, 3rd 
Series, Vol. 24. 


Danish Hydraulic Institute. 1999. MIKE SHE Water Movement. [CITY] Denmark. 


Davis, S.M. and J.C. Ogden. 1994. Everglades: The Ecosystem and Its Restoration. St. 
Lucie Press, Delray Beach, FL. 


Draft 201 


DRAFT LEC subRegional MODFLOW Model Documentation 


EAS Engineering, Inc. 1996. Vegetation and Soils Mapping and Analysis, Lakebelt 
Ecological Studies, Dade County. Final Report to Dade County Department of 
Environmental Resource Management, [City], FL. 


Fish, J.E. 1988. Hydrogeology, Aquifer Characteristics, and Ground-Water Flow of the 
Surficial Aquifer System, Broward County, Florida. U.S. Geological Survey Water- 
Resources Investigations Report 87-4034, Books and Open-File Reports Section, 
USGS, Federal Center, Denver, CO, 92 pp., 10 sheets. 


Fish, J.F. and M.J. Stewart. 1991. Hydrogeology of the Surficial Aquifer System, Dade 
County, Florida. U.S. Geological Survey Water-Resources Investigations Report 90- 
4108. Books and Open-File Reports Section, USGS, Federal Center, Denver, CO, 
50 pp., 11 sheets. 


Freeze, R.A. and J.A. Cherry. 1979. Groundwater. Prentice-Hall, Englewoods Cliffs, NJ. 


Gavaskar, A., N. Gupta, B. Sass, R. Janosy and J. Hicks. 2000. Design Guidance for 
Application of Permeable Reactive Barriers for Groundwater Remediation. Prepared 
by Battelle, Columbus, Ohio, for the Strategic Environmental Research and 
Development Program, U.S. Department of Defense. 


GeoTrans. 2004. Surface/Vadose-Zone Package for Modflow User's Manual. GeoTrans, 
Inc., Roswell, GA. 


Giddings, J. B., R.L. Restepo. 1995. Surface/Vadose-Zone Package for Modflow User's 
Manual. SFWMD, West Palm Beach, FL. 


Giddings, J.B. and L.L. Kuebler. 2004. Seepage Mitigation Analysis for the Site 1 
Impoundment Utilizing the Lower East Coast subRegional MODFLOW Model. 
Model Application Section, Water Supply Planning and Development Division, 
South Florida Water Management District, West Palm Beach, FL. 


Giddings, J.B. 1999. Deposition and Chronostratigraphic Subdivision of the Quaternary 
Sediments Broward County, Florida. Florida Atlantic University, Boca Raton, FL. 


Griswold, L.S. 1896. Notes on the geology of southern Florida. Bulletin of the Museum 
of Comparative Zoology, XXVIII, 52-59. 


Hammer, D.E. and R. H. Kadlec. 1986. 4 Model for Wetland Surface Water Dynamics. 
Water Resources Research, 22(13):1951-1958. 


Hagen, S.C. 2005. Long-Term Tidal Analysis for Subregional Modeling. Winter Park, FL. 
38 pp. 


Hantush, M.S. 1964. Flow to Water-Table Wells. Advances in Hydroscience, 1:358-373. 


Harbaugh, A.W. and McDonald, M.G. 1996. Documentation for MODFLOW-96, and 
Update to the U.S. Geological Survey Modular Finite-Difference Ground Water Flow 


Draft - 202 


LEC subRegional MODFLOW Model Documentation DRAFT 


Model. Techniques of Water-Resources Investigations of the USGS. Open-File 
Report 96-485, USGS, Reston, VA. 


Harvey, J.W., S.L. Krupa, C. Gefvert, R.H. Mooney, J. Choi, S.A. King and J.B. 
Giddings. 2002. Interactions between Surface Water and Ground Water and Effects 
on Mercury Transport in the North-central Everglades. U.S. Geological Survey 
Water-Resources Investigations Report 02-4050. USGS, Tallahassee, FL, 82 pp. 


Hearty, P.J.; Vacher, H.L.; and Mitterer, R.M., 1992. Aminostratigraphy and ages of 
Pleistocene limestones of Bermuda. Geological Society of America Bulletin, 104, 
471-480. 


Heatherington, A.L. and P.A. Mueller. 1997. Geochemistry and Origin of Florida Crustal 
Basement Terrenes. In The Geology of Florida, eds. A.F. Randazzo and D.S. Jones. 
[PUBLISHER] Gainesville, Florida, pp. 27-38. 


Hill, M.C. 1992. A Computer Program (MODFLOWP) for Estimating Parameters of a 
Transient, Three-Dimensional, Ground Water Flow Model Using | Nonlinear 
Regression. U.S. Geological Survey Water-Resources Investigations Report 91-484. 
Books and Open-File Reports Section, USGS, Federal Center, Denver, CO, 358 pp. 


Hoffmeister, J.E., K.W. Stockman, and H.G. Multer. 1967. Miami Limestone of Florida 
and its recent [Bahamian] counterpart. Geologic Society of American Bulletin, 
78:17-190. 


Hoffmeister, J.E. 1974. Land from the Sea, University of Miami Press, Coral Gables, FL, 
143 pp. 


Hopkins, E. 1991. A Water Resource Analysis of the Jensen Beach Pensinsula, Martin 
County, Florida. South Florida Water Management District Technical Publication, 
DRE-301, SFWMD, West Palm Beach, FL. 


Hsieh, PA. and J.R. Freckleton. 1993. Documentation of a computer program to simulate 
horizontal-flow barriers using the U.S. Geological Survey modular three-dimensional 
finite-difference ground-water flow model. U.S. Geological Survey Open-File Report 
92-477, 32 pp. 


Irizarry-Ortiz, M.M. 2003. Memorandum. Selected Methodology for Long-Term (1965- 
2000) Solar Radiation and Potential Evapotranspiration Estimation for the SFWMM 
2000 Update. SFWMD, West Palm Beach, FL. 


Irizarry-Ortiz, M.M. 2003. Memorandum: Review of Methods for Long-Term (1965- 
2000) Solar Radiation and Potential Evapotranspiration Estimation for 
HydrologicModeling in South Florida. SEWMD, West Palm Beach, FL. 


Jensen, M.E., Burman, R.D., and Allen, R.G. (eds) 1990. Evaporation and Irrigation 
Water Requirements. ASCE Manuals and Reports No. 70, 350 pp. 


Draft 203 


DRAFT LEC subRegional MODFLOW Model Documentation 


Jones, M. 1999, 4 Reinjection Drainflow Package for MODFLOW. Albuquerque, NM, 
33 pp. 


Kadlec, R.H. 1990. Overland Flow in Wetlands: Vegetation Resistance. Journal of 
Hydraulic Engineering, 116:691-705. 


Kadlec, R.H., D.E. Hammer, 1S. Nam and J.O. Wilkes. 1981. The Hydrology of 
Overland Flow in Wetlands. Chem. Engr. Commun., 9:33 1-334. 


Keith and Schnars Survey, Inc. 1993. Survey of WCA-2 conducted by the Rock Miners 
c/o Paul Larsen. 


Keys, S.W. 1989. Borehole Geophysics Applied to Ground-Water Investigations. 
National Water Well Association, Dublin, OH, 313 pp. 


Klein, H. and C.B. Sherwood. 1961. Hydrologic Conditions in the Vicinity of Levee 30, 
Northern Dade County, Florida. Report of Investigations No. 24, Part I, Florida 
Geological Survey, Tallahassee, FL, 24 pp. 


Kohout, F.A. 1964. The Flow of Fresh Water and Salt Water in the Biscayne aquifer of 
the Miami Area, Florida. Seawater in Coastal Aquifers: United States Geological 
Survey Water Supply Paper 1613-C, Denver, CO. 


Krishnan, N. and C.A. Gove. 1993. South Florida Water Management Model (SFWMM) 
Coastal Boundary Conditions. Internal Memorandum, Lower East Coast Planning 
Division, South Florida Water Management District, West Palm Beach, FL. 


Krupa, A.J. 1997. Results of Literature Review and Assessment of Geologic Logs to 
Determine the Extent of the Dense Limestone Layer in the Upper Portion of the 
Biscayne Aquifer in the Pennsuco Area, Dade County. Draft Memorandum, Planning 
Department, South Florida Water Management District, West Palm Beach, FL, 27 pp. 


Land, L.F. 1972. Progress Report on Water-Resource Investigation in Palm Beach 
County, Florida. USGS Open file report 72-014, Tallahasee, FL. 


Langevin, C.D. 2001. Simulation of Ground-Water Discharge to Biscayne Bay, 
Southeastern Florida. Water-Resources Investigations Report 00-4251, USGS, 
Tallahassee, FL. 


Langevin, C.D., Shoemaker, W.B., and Guo, W. 2003. MODFLOW-00, the U.S. 
Geological Survey Modular Ground-Water Model — Documentation of the SEAWAT- 
2000 Version with the Variable-Density Flow Process (VDF) and the Integrated 
MT3DMS Transport Process (IMT). Open-File Report 03-426, USGS, Tallahassee, 
FL. 


Leach, S.D. and C.B. Sherwood. 1963. Hydrologic Studies in the Snake Creek Canal 
Area, Dade County, Florida. Report of Investigations No. 24, Part III, Florida 
Geological Survey, Tallahassee, FL, 33 pp. 


Draft - 204 


LEC subRegional MODFLOW Model Documentation DRAFT 


Lewandowski, J.A. 1993. Vegetation resistance and circulation modeling in a tidal 
wetland. Ph.D. diss., University of California, Berkeley, CA. 


Lichtler, W.F. 1960. Geology and ground-water resources of Martin County, Florida 
Report of Investigations No. 23. Florida Geological Survey, Tallahassee, FL, 149 pp. 


Lovejoy, D.W. 1983. The Anastasia Formation in Palm Beach and Martin Counties, 
Florida. Miami Geological Society Memoir, 3:58-72. 


Lukasiewicz, J. and Adams Smith, K. 1996. Hydrogeologic data and information 
collected from the surficial and Floridan aquifer systems, UEC planning area - Part 
1, Text. South Florida Water Management District Technical Publication, WRE-337, 
SEWMD, West Palm Beach, FL. 


MacVicar, T.K., T. Van Lent and A. Castro. 1984. South Florida Water Management 
Model Documentation Report. SFWMD Technical Publication 84-3. SFWMD, West 
Palm Beach, FL. 


Mansfield, W.C. 1939. Notes on the Upper Tertiary and Pleistocene Mollusks of 
Peninsular Florida. Florida Geological Survey, Geological Bulletin No. 18. 


Matson, G.C. and F.G. Clapp. 1909. A Preliminary Report on the Geology of Florida. 
Florida Geological Survey, 2nd Annual Report. 


McCollum, S.H., O.E. Cruz, L.T. Stem, W.H. Wittstruck, R.D. Ford and F.C. Watts. 
1978. Soil Survey of Palm Beach County Area, Florida. U.S. Department of 
Agriculture, Soil Conservation Service, [City, ST]. (get more current one) 


McDonald, M.G. and A.W. Harbaugh. 1988. A Modular, Three-Dimensional Finite- 
Difference Ground Water Flow Model. Techniques of  Water-Resources 
Investigations of the U.S. Geological Survey. Book 6, U.S. Government Printing 
Office, Washington, DC. 


Merritt, M.L. 1995. Simulation of the Water Table Altitude in the Biscayne Aquifer, 
Southern [Miami-[Dade County, Florida, Water Years 1945-89. U.S. Geological 
Survey Open File Report 95-337. U.S. Geological Survey Earth Science Information 
Center, Open-File Reports Section, Denver, CO, 88 pp., 9 sheets, 50 figures, 12 
tables. 


Merritt, M.L., 1996. Assessment of saltwater intrusion in southern coastal Broward 
County, Florida. U.S. Geological Survey Water-Resources Investigations Report, 96- 
4221, 133 pp. 


Miller, W.C. 1978. Effects of Bottom Sediments on Infiltration from the Miami and 
Tributary Canals to the Biscayne Aquifer, Dade County, Florida. Water-Resources 
Investigations Report 78-36, Books and Open-File Reports Section, United States 
Geological Survey, Federal Center, Denver, CO, 63 pp. 


Draft 205 


DRAFT LEC subRegional MODFLOW Model Documentation 


Miller, W.A. 1987. Lithology and Base of the Surficial Aquifer System, Palm Beach 
County, Florida. U.S. Geological Survey, Water-Resources Investigations Report 86- 
4067, USGS, Tallahassee, FL. 1 sheet. 


Mitterer, R.M. 1975. Ages and digenetic temperatures of Pleistocene deposits of Florida 
based upon isoleucine epimerization in Mercenaria. Earth and Planetary Science 
Letters, 28:275-282. 


Mock, Ross Associated Inc. 1992. Lake Worth Drainage Canal System. Ref. No. 
94730.98. West Palm Beach, FL. 


Mock, Ross Associated Inc. 1996. Public Facilities Report for Lake Worth Drainage 
System. West Palm Beach, FL. 


Model Application Group. 2003. L31N Seepage Management Pilot Project Study 
Utilizing the North Miami-Dade Subregional Ground Water Model. Water Supply 
Planning and Development Division, South Florida Water Management District, 
West Palm Beach, FL. 45 pp. 


Mortellaro, S. et al. [Other authors] 1995. Literature Review on the Effects of 
Groundwater Drawdowns on Isolated Wetlands. South Florida Water Management 
District Technical Publication, WRE-330. 


Nair, S.K. and M.M. Wilsnack. 1998. A Comparison of Two Approaches for Simulating 
Lake-Ground Water Interactions with MODFLOW. Proceedings for MODFLOW '98, 
Colorado School of Mines, October 4 — 8, Volume II, pp. 871 — 878. 


Nair, S.K., A.M. Montoya, D.E. Welter, M.M. Wilsnack, M.P. Switanek, J. Herr, D.M. 
Adam, L.M. Zamorano, J.I. Restrepo and J.T.B. Obeysekera (2000). South Palm 
Beach County Ground Water Flow Model. Hydrologic Systems Modeling 
Department, South Florida Water Management District, West Palm Beach, Florida. 


Neuman, S.P. 1975. Analysis of Pumping Test Data from Anisotropic Unconfined 
Aquifers Considering Delayed Gravity Response. Water Resources Research, 
II(2):329-342. 


NOAA. 1988. Web page: HYPERLINK "http://www.co- 
ops.nos.noaa.gov/data res.html" www.co-ops.nos.noaa.gov/data res.html 


NOAA. 1989. Web page: HYPERLINK "http://www.co- 
ops.nos.noaa.gov/data res.html" www.co-ops.nos.noaa.gov/data res.html 


NOAA. 1990. Web page: HYPERLINK "http://www.co- 
ops.nos.noaa.gov/data res.html" www.co-ops.nos.noaa.gov/data res bt! 


NOAA. 2000. Web page: HYPERLINK "http://www.co- 
ops.nos.noaa.gov/data res.html" 


Draft - 206 


LEC subRegional MODFLOW Model Documentation DRAFT 


HYPERLINK "http://www.co-ops.nos.noaa.gov/data res.html" 
WWW.co-ops.nos.noaa.gov/data res.html 


Osmond, J.K., J.R. Carpenter, and H.L. Windom. 1965. Th230/U234 age of Pleistocene 
corals and oolites of Florida. Journal of Geophysical Research, 70:1843-1847. 


Parker, G.G., and C.W. Cooke. 1944. Late Cenozoic geology of southern Florida, with a 
discussion of the ground water. Florida Geological Survey Geological Bulletin, 27, 
Tallahassee, FL, 119 pp. 


Parker, G.G., G.E. Ferguson and S.K. Love. 1944. Interim Report on the Investigations of 
Water Resources in Southeastern Florida with Special Reference to the Miami Area 
in Dade County. Report of Investigations No. 4, Florida Geological Survey, 
Tallahassee, FL, 39 pp. 


Parker, G.G., G.E. Ferguson and S.K. Love [et al. (eds.)] 1955. Water Resources of 
Southeastern Florida, with Special Reference to the Geology and Groundwater of the 
Miami Area. U.S. Geological Survey, Water-Supply Paper 1255, U.S.G.S., U.S. 
Government Printing Office, Washington, DC, 965 pp., 22 plates. 


Pendleton, R.F., H.D. Dollar and L. Law, Jr., 1976a. Soil Survey of Broward County 
Area, Florida. U.S.D.A. Soil Conservation Service, Tallahassee, FL. 


Pendleton, R.F., H.D. Dollar and L. Law, Jr., 1976b. Soil Survey of [Miami-] Dade 
County Area, Florida. U.S.D.A. Soil Conservation Service. [CITY, State] 


Perkins, R.D. 1968. Late Cenozoic Stratigraphy of Southern Florida — A Reappraisal. 
Second Annual Field Trip of the Miami Geological Society. Miami Geological 
Survey, Miami, Florida [pp. ??]. 


Perkins, R.D. 1977. Depositional framework of Pleistocene rocks in south Florida. In: 
Enos, P. and R.D. Perkins (eds.), Quaternary Sedimentation in South Florida, Part II. 
Geological Society of American Memoir, 147:131-198. 


Poeter, E.P. and M.C. Hill. 1998. Documentation of UCODE, A Computer Code for 
Universal Inverse Modeling. U.S. Geological Survey, Water-Resources Investigations 
Report 98-4080, Branch of Information Services, USGS, Federal Center, Denver, CO, 


116 pp. 


Randall, D. 1999. User's Documentation for the MODFLOW96 Trigger Package. Water 
Resources Management, Inc. 13 pp. 


Randazzo, A.F. 1997. The Sedimentary Platform of Florida: Mesozoic to Cenozoic. In 
The Geology of Florida, eds. A.F. Randazzo and D.S. Jones. University Press of 
Florida, Gainesville, FL, pp. 39-56. 


Draft 207 


DRAFT LEC subRegional MODFLOW Model Documentation 


Reese, R.S. and K.J. Cunningham. 2000. Hydrogeology of the Gray Limestone Aquifer in 
Southern Florida. U.S. Geological Survey, Water-Resources Investigations Report 
99-4213, Branch of Information Services, USGS, Federal Center, Denver, CO. 
244 pp. 


Renken, Robert A., Dixon, Joann, Koehmstedt, John, Ishman, Scott, Lietz, A.C., Marella, 
Richard L., Telis, Pamela, Rodgers, Jeff, and Memberg, Steven, 2005, Impact of 
Anthropogenic Development on Coastal Ground-Water Hydrology in Southeastern 
Florida, 1900-2000, U.S. Geological Survey Circular 1275, 77p. + CD. 


Restrepo, J.I. 1998. Personal communication with Restrepo, J.I. 1998. 


Restrepo, J.I. and J.B. Giddings. 1994. Physical[ly] Based Methods to Estimate ET 
[Evapotranspiration] and Recharge Rates Using GIS. In: Effects of Human-Induced 
Changes on Hydrologic Systems, American Water Resources Association, Bethesda, 
MD. [pp. 65-74] 


Restrepo, J.L. and A.M Montoya. 1997. MODFLOW Wetland Module: Final Report. 
Department of Geography and Geology, Florida Atlantic University, Boca Raton, FL, 


86 pp. 


Restrepo, J., Bevier, C. and Butler, D. 1992. A Three-Dimensional Finite Difference 
Ground Water Flow Model of the Surficial Aquifer System, Broward County, Florida. 
South Florida Water Management District Technical Publication, DRE-314. [Is this 
the same report as below? | 


Restrepo J.I., C. Bevier and D. Butler. 1992. A three dimensional finite difference ground 
water flow model of the surficial aquifer system, Broward County, Florida. Technical 
publication 92-05, West Palm Beach, Florida. South Florida Water Management 
District. 261 pp. [or 231 pp.?] 


Restrepo, J.I., A.M. Montoya and D. Garces. 2003. The Updated Utility Generation 
Package (USGEN) for the Groundwater Flow Model, MODFLOW. Developed by the 
Hydrological Research Center at Florida Atlantic University, Boca Raton, FL in 
cooperation with the South Florida Water Management District, West Palm Beach, 
FL. 


Restrepo, J.I., A.M. Montoya and J. Obeysekera. 1998. A Wetland Simulation Module for 
the MODFLOW Ground Water Model. Ground Water, 36(5):764-770. 


Restrepo, J.I., D. Garces, A. Montoya, J. Giddings and N. Restrepo. 2001. A three- 
dimensional finite difference groundwater flow model of the surficial aquifer system, 
south Miami-Dade County, Florida. Florida Atlantic University, Hydrological 
Modeling Center, in collaboration with the South Florida Water Management District, 


49 pp. 


Restrepo, J.1., J. Giddings, D. Garces and N. Restrepo. 2001. A three-dimensional finite 
difference ground water flow model of the surficial aquifer system, Broward County, 


Draft - 208 


LEC subRegional MODFLOW Model Documentation DRAFT 


Florida. Florida Atlantic University, Hydrological Modeling Center, in collaboration 
with the South Florida Water Management District, 51+ pp. 


Richter, W., E. Myers and K. Fanning. 1990. Bird Drive Everglades Basin Special Area 
Management Plan: Baseline Studies and Resource Evaluation. Dade County 
Department of Environmental Resources Management, Miami, FL, 11 pp. 


Russell, G.M. and E.J. Wexler. 1993. Hydrogeology and Simulation of Ground Water 
Flow near the Lantana Landfill, Palm Beach County, Florida. U.S. Geological 
Survey Water-Resources Investigations Report 92-4107, USGS, Palm Beach County, 
FL, 55 pp. 


Sanford, S. 1909. Geology of southern Florida. Florida Geological Survey 2nd Annual 
Report, Tallahassee, FL, pp. 177-231. 


Schroeder, M.C. 1954. Water Resource Studies: Water Resources of Palm Beach County, 
Florida. Report of Investigations No. 13. Florida Geological Survey, Tallahassee, FL, 


63 pp. 


Schroeder, M.C., H. Klein and N.D. Hoy. 1958. Biscayne Aquifer of Dade and Broward 
Counties, Florida. Report of Investigations No. 17, Florida Geological Survey, 
Tallahassee, FL, 56 pp. 


Scott, Thomas, 1992, A geological overview of Florida, Florida Geological Survey, Open 
File Report 50, 78p 


Scott, Thomas, 1990, The lithostratigraphy of the Hawthorn [i.e. Hawthorne] Group of 
peninsular Florida, Florida Geological Survey, Open File Report 36, 327-336 p 


Sellards, E.H. 1912. The Soils and Other Surface Residual Materials of Florida. In: 
Sellards, E.H. (ed.), Florida Geological Survey Annual Report 4; 1-79. 


Shaler, N.S. 1890. The topography of Florida. Bulletin of the Museum of Comparative 
Zoology, 16(7):143-144. 


Sherwood, C.B. and S.D. Leach. 1962. Hydrologic Studies in the Snapper Creek Canal 
Area, Dade County, Florida. Report of Investigations No. 24, Part II, Florida 
Geological Survey, Tallahassee, FL, 39 pp. 


Shine, M.J., D.G.J. Padgett, and W.M. Barfknecht. 1989. Ground Water Resource 
Assessment of Eastern Palm Beach County, Florida, Part I, Text. Technical 
Publication 89-4, DRE-278. South Florida Water Management District, West Palm 
Beach, FL. | 


*Shine, M.J., D.G.J. Padgett, and W.M. Barfknecht. 1989. Ground Water Resource 
Assessment of Eastern Palm Beach County, Florida, Part 11, Plates and Oversized 
Appendices. Technical Publication 89-4, DRE-279. South Florida Water 
Management District, West Palm Beach, FL. 


Draft 209 


DRAFT LEC subRegional MODFLOW Model Documentation 


Sklar, F., C. McVoy, R. VanZee, D.E. Gawlik, K. Tarboten, D. Rudnick, S. Miao and T. 
Armentano. 2002. Chapter 2: The Effects of Altered Hydrology on the Ecology of the 
Everglades. In The Everglades, Florida Bay, and Coral Reefs of the Florida Keys: an 
Ecosystem Sourcebook, eds. J.W. Porter and K.G. Porter. CRC Press, Boca Raton, 
FL, pp. xx-xx. 


Smajstrla, A.G. 1990. Agricultural Field Scale Irrigation Requirements Simulation 
(AFSIRS) Model, Version 5.5. Agricultural Engineering Department, University of 
Florida, Gainesville, FL. 


Sonenshein R. 1995. Delineation of extent saltwater intrusion in the Biscayne aquifer, 
eastern [Miami-]Dade County. U.S. Geological Survey Water-Resources 
Investigation. Report 96-4285, USGS [City, ST]. 


Sonenshein, R.S. and E.J. Koszalka. 1996. Trends in water-table altitude (1984-93) and 
saltwater intrusion (1974-93) in the Biscayne aquifer, [Miami-]Dade County, Florida. 
U.S. Geological Survey Open-File Report 95-705, 2 sheets. 


Sosnowski, R. 1994. Elevation Coverage Documentation. Draft Internal Report. Water 
Resources and Evaluation Department, South Florida Water Management District, 
West Palm Beach, FL, 14 pp. 


South Florida Water Management District. 1980. Interim Final Draft Report on 
Embankment Failure Florida Power and Light Company Martin Plan Cooling 
Reservoir. SFWMD Technical Publication, DRE-102, West Palm Beach, FL. 


South Florida Water Management District. 1996. Draft Data Analysis Task Report for 
ACME Improvement District. Contract No. C-3021 Amendment No 5, SFWMD, 
West Palm Beach, FL. 


South Florida Water Management District. 1996. Northwest Dade County Freshwater 
Lake Belt Plan: Making a Whole, Not Just Holes. Planning Department, SFWMD, 
West Palm Beach, FL, 29 pp. 


South Florida Water Management District. 2000a. Lower East Coast Regional Water 
Supply Plan. Water Supply Planning and Development Department, SFWMD, West 
Palm Beach, FL. 


South Florida Water Management District. 2000b. Minimum Flows and Levels for Lake 
Okeechobee, the Everglades, and the Biscayne Aquifer. Water Supply Department, 
SFWMD, West Palm Beach, FL. 


South Florida Water Management District. 2001. 4 Three-Dimensional Finite Difference 
Ground-Water Flow Model of the Surficial Aquifer System, Northern Palm Beach 
County, Florida, July 27 (Draft). SFWMD, West Palm Beach, FL. 


South Florida Water Management District. 2002. Technical Documentation to Support 
Development of Minimum Flows and Levels for the Northwest Fork of the 


Draft - 210 


LEC subRegional MODFLOW Model Documentation DRAFT 


Loxahatchee River and Estuary. Water Supply Department, SFWMD, West Palm 
Beach, FL. 


South Florida Water Management District. HYPERLINK http://www.sfwmd.gov/ 


South Florida Water Management District. Diversion Systems: Structures and Pump 
Stations. Prepared by Florida Atlantic University for the South Florida Water 
Management District. [Online] Available http://iweb/iwebB501/omod/hsm/modflow/ 
div system.htm, December 2005. 


Stephenson, L.W. 1912. Coastal Plain of North Carolina. North Carolina Geological 
Survey, 3: 266-290. 


Stodghill, A.M. 1983. Resistivity investigation of the Coastal Ridge Aquifer 
hydrostratigraphy, Martin County, Florida [electronic [resource]] / by Allan M. 
Stodghill. 


Swayze, L.J. and W.L. Miller. 1984. Hydrogeology of a Zone of Secondary Permeability 
in the Surficial Aquifer of Eastern Palm Beach County, Florida. Water-Resources 
Investigations Report 83-4249. U.S. Geological Survey, Palm Beach County, FL. 


Switanek, M.P. 2000. Data Acquisition, Review, and Analysis for the Lake Belt and 
Surrounding Areas, Miami-Dade County, Florida. Initial Draft Report. Water Supply 
Department, South Florida Water Management District, West Palm Beach, FL, 15 pp. 


Tuomey, M. 1851. Notice of the geology of the Florida Keys and of the southern coast of 
Florida. American Journal of Science, 11:390-394. 


University of Florida. 1992. Survey of WCA-1 conducted by Dr. Richardson. Arc/Info 
Grid coverage. 


U.S. Army Corps of Engineers. 1951. Central and Southern Florida Project - Part 4, 
Agricultural and Conservation Areas Supplement 7 - Design Memorandum 
Permeability Investigations by Well Pumping Tests. Office of the District Engineer, 
USACE, Jacksonville, FL. 


U.S. Army Corps of Engineers. 1960. Central and Southern Florida Project - Part 4, 
Agricultural and Conservation Areas Supplement 33 - General Design Memorandum, 
Conservation Area No. 3. June 22. USACE, Jacksonville, FL. 


U.S. Army Corps of Engineers. 1987. Wetlands Delineation Manual. Wetlands Research 
Program Technical Report Y-87-1 (on-line edition). USACE. 


U.S. Army Corps of Engineers and South Florida Water Management District. 1999, 
Central and Southern Florida [Flood Control] Project Comprehensive Review Study 
Final Integrated Feasibility Report and Programmatic Environmental Impact 
Statement. USACE, Jacksonville, District, Jacksonville, FL, and SFWMD, West Palm 
Beach, FL. 


Draft 211 


DRAFT LEC subRegional MODFLOW Model Documentation 


U.S. Department of Agriculture. 1978. Soil Survey of Palm Beach County, Florida. 
USDA, Soil Conservation Service [City, ST]. 


U.S. Geological Survey. 1988. Water Resources Data Florida Water Year 1988. South 
Florida Ground Water. Water-Data Report FL-88. Vol 2B, prepared in cooperation 
with the state of Florida and with other agencies. USGS, Miami, FL. 


U.S. Geological Survey. 1989. Water Resources Data Florida Water Year 1989. South 
Florida Ground Water. Water-Data Report FL-89. Vol 2B, prepared in cooperation 
with the state of Florida and with other agencies. USGS, Miami, FL. 


U.S. Geological Survey. 1990. Water Resources Data Florida Water Year 1990. South 
Florida Ground Water. Water-Data Report FL-90. Vol 2B, prepared in cooperation 
with the state of Florida and with other agencies. USGS, Miami, FL. 


U.S. Geological Survey. 1994. Water Resources Data Florida Water Year 1994. South 
Florida Ground Water. Water-Data Report FL-94-2B, prepared in cooperation with 
the state of Florida and with other agencies. USGS, Miami, FL. 


U.S. Geological Survey. 1995. Water Resources Data Florida Water Year 1995. South 
Florida Ground Water. Water-Data Report FL-95, prepared in cooperation with the 
state of Florida and with other agencies. USGS, Miami, FL. 


Walton, W.C. 1970. Groundwater Resource Evaluation. McGraw-Hill, New York, NY. 
Walton, W.C. 1987. Groundwater Pumping Tests. Lewis Publishers, Chelsea, MI. 


Wilsnack, M.M. 1993. A Compilation of Attributes for Selected Secondary Canals in 
Dade County. Lower East Coast Planning Division, South Florida Water 
Management District, West Palm Beach, FL, 17 pp. 


Wilsnack, M.M. 1995a. A Compilation of Attributes for District Canals in Dade County. 
Lower East Coast Planning Division, South Florida Water Management District, 
West Palm Beach, FL, 17 pp. 


Wilsnack, M.M. 1995b. The Use of ARC/INFO and MODFLOW to Evaluate the 
Feasibility of Using Limestone Mining Quarries in Southern Florida as Ground 
Water Recharge Basins. Proceedings of the Seventh Biennial Symposium on the 
Artificial Recharge of Ground Water: the Role of Recharge in Integrated Water 
Management, Tempe, AZ, pp. 199—213. 


Wilsnack, M.M. 1998. GIS Techniques for Incorporating Surface Water Features into 
MODFLOW Models. Hydrologic Systems Modeling Division, South Florida Water 
Management District, West Palm Beach, FL, 20 pp. 


Wilsnack, M.M. and S.K. Nair. 1998. Designing an ARC/INFO Database for a 
MODFLOW Model: A Case Study in the Limestone Mining Belt of Southern Florida. 


Draft - 212 


LEC subRegional MODFLOW Model Documentation DRAFT 


Proceedings for MODFLOW ’98, Colorado School of Mines, October 4-8, Volume 
II, pp. 649-656. 


Wilsnack, M.M., S.K. Nair and J. Obeysekera. 1997. Hydrologic Modeling of Initial Lake 
Belt Alternative Configurations: Ground Water Modeling. Hydrologic Systems 
Modeling Division, South Florida Water Management District, West Palm Beach, FL, 


66 pp. 


Wilsnack, M.M., D.E. Welter, A.M. Montoya, J.I. Restrepo and J. Obeysekera. 2001. 
Simulating Flow in Regional Wetlands with the MODFLOW Wetlands Package. 
Journal of the American Water Resources Association, 37(3):655-674. 


Wilsnack, M.M., D.E. Welter, S.K. Nair, A.M. Montoya, L.M. Zamorano, J.I. Restrepo 
and J.T.B. Obeysekera (2000). North Miami-Dade County Ground Water Flow 
Model. Hydrologic Systems Modeling Department, South Florida Water Management 
District, West Palm Beach, Florida, 40 pp. 


Yan, J.J., C. Bevier and K.A. Smith. 1993. Revision and Recalibration of the Three- 
Dimensional Ground Water Model of the Palm Beach County, Florida. South Florida 
Water Management District Technical Publication April 1993. SFWMD, West Palm 
Beach, FL. 


Yan, J.J., C. Bevier and K.A. Smith. 1995. Revision and Recalibration of the Three- 
Dimensional Groundwater Flow Model of Palm Beach County, Florida. South 
Florida Water Management District Technical Publication, DRE-322. SFWMD, West 
Palm Beach, FL. 


Draft 243 


LEC subRegional MODEL OW Model Documentation DRAFT 


APPENDIX A 
Hydrologic and Hydrologic Data 


Draft 1 


DRAFT LEC subRegional MODFLOW Model Documentation 
Table of Tables 
Table A.1. LECsR Biscayne or Surficial Aquifer Performance Test Results....................................... 3 
Table A.2. Aquifer Performance Test Results Based on Specific Capacity Well Measurements....6 
Table A.3. Summary Table of USGS Geologic Core Data Used in LECsR Model ............................. 9 
Table A.4. Conceptual Model layer Properties 23 
Table A.5. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 1986 
SOUS SFWMM E ETT 26 
Table A.6. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 1987 
(SOU E WMM Suet Edesa 27 
Table A.7. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 1988 
Ee SIE E WMM po ET 28 
Table A.8. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 1989 
(Source SEW MM ) cS 29 
Table A.9. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) bv ET Station for 1990 
Sonce SPWMM mE 30 
Table A.10. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
[991 (Source: SEWM M iaa ad 31 
Table A.11. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
[99 Z2 (SGU SE WMV oia 32 
Table A.12. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
OS eT c WME a aperta ee aia 33 
Table A.13. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
1994 (Source SEW MM EE 34 
Tabie A.14. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
1995 (Source -SEWNMUM.) aaa 35 
Table A.15. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
1996 (Source AT EE 36 
Table A.16. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
1997 (Source: SEW MM) O isis 37 
Table A.17. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
¡SOS Solier SF WMM) c E 38 
Table A.18. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) by ET Station for 
1999 (SOUT EEST EE 39 
Table A.19. LECsR Wet Marsh Potential Evapotranspiration Rates (in Inches) bv ET Station for 
ŻUQUSOviċe SEW MM aia 이이 40 
Table A.20. LECsR Precipitation (in Inches) by Rain Gauge Station for 1986............................... 41 
Table A.21. LECsR Precipitation (in Inches) by Rain Gauge Station for 1987............................... 42 
Table A.22. LECsR Precipitation (in Inches) by Rain Gauge Station for 1988............................... 43 
Table A.23. LECsR Precipitation (in Inches) by Rain Gauge Station for 1989............................... Ad 
Table A.24. LECsR Precipitation (in Inches) by Rain Gauge Station for 1990............................... 46 
Table A.25. LECsR Precipitation (in Inches) by Rain Gauge Station for 1991............................... 47 
Table A.26. LECSR Precipitation (in Inches) by Rain Gauge Station for 1992............................... 48 
Table A.27. LECsR Precipitation (in Inches) by Rain Gauge Station for 1993............................... 50 
Table A.28. LECsR Precipitation (in Inches) by Rain Gauge Station for 1994............................... 51 
Table A.29. LECSR Precipitation (in Inches) by Rain Gauge Station for 1995............................... 53 
Table A.30. LECsR Precipitation (in Inches) by Rain Gauge Station for 1996............................... 55 
Table A.31. LECsR Precipitation (in Inches) by Rain Gauge Station for 1997............................... 57 
Table A.32. LECsR Precipitation (in Inches) by Rain Gauge Station for 2000............................... 61 
Draft - 2 


LEC subRegional MODFLOW Model Documentation 


Draft 





DRAFT 


Table 27. LECsR Biscayne or Surficial Aquifer Performance Test 
Results 
Well Location? 
Well Name X Y HK” (ft/day) ` Thickness (ft) 
DADE-52 944824 566710 6357 40 
DADE-53 924256 569405 5000 50 
DADE-54 903062 570691 6750 40 
DADE-51 924256 569405 10000 40 
DADE-61 912358 589223 5500 50 
DADE-63 901469 593803 2600 50 
DADE-CC 802783 600824 6000 40 
HIGH-04 958794 756061 1250 50 
PB-1574 916970 759977 480 64 
Morikam 932971 759977 1300 50 
PB-1581 927970 739977 925 95 
NSID 885970 709977 100 50 
G1260 945912 722211 200 50 
G1213 922626 713008 300 50 
G1315 934080 710491 10 50 
G2147 950443 697779 20000 50 
G1316 922626 695513 250 50 
G1215 946815 708363 2500 50 
G820A 930721 678308 250 50 
PB491 952750 735464 1250 50 
PB1455 946922 735340 400 50 
PB145/ 957657 747485 1250 50 
PB1158 957657 754696 1250 50 
McCart 819085 1058463 33 40 
Allap 824356 1002412 58 60 
cauk 860786 997062 50 120 
jdsp 945537 973063 150 160 
jup 939904 943310 100 165 
pb15 906536 945312 35 76 
pbpc 892636 935362 36 55 
teq 952737 959763 250 220 
fpI2 808494 977251 53 45 
jdsp 928086 979688 100 90 
mobil 907136 1003863 80 53 
monrv 884286 939663 56 60 
mapwd 941987 939363 58 90 
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Well Name 
jonl 
mecca 
hobe4 
hobed 
acme 
palms 
boyn 
sys3 
lant 
sshor 
pb1598 
pb1603 
pb1578 
pb1571 
pb1544 
pb1576 
i95 
hood 
wpb 
jup 

riv 
pb1567 
pb1564 
pb1558 
pb1555 
pb1550 
pb1607 
lakel 
stlucie 
banyan 
MilesG 
Hydra 
Hobe 
Vista 
C23 
L65 
MDown 
LFarm 
Stuso 
Woods 
Piper 
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Well Location? 


X 
951537 
943837 
937986 
937136 
910909 
940419 
940801 
936822 
949266 
903754 
927408 
948839 
954014 
926587 
926983 
909055 
950360 
940246 
942435 
939904 
947545 
892975 
899423 
876043 
912643 
887449 
924443 
949266 
885536 
900636 
922236 
932086 
936086 
917636 
797236 
770636 
878086 
879523 
905486 
890086 
893586 


Y 
939663 
930362 
988363 
989063 
829724 
841728 
798815 
182632 
852089 
842403 
776815 
805434 
831220 
819624 
842649 
805822 
898648 
919985 
878397 
943310 
886611 
854962 
884686 
900610 
901730 
918848 
926544 
852089 
1004613 
1022013 
1022413 
1005463 
996363 
1015763 
1044063 
1008762 
1036513 
1018913 
1029513 
1041613 
1023663 


HK” (ft/day) 
91 
57 

126 
76 
130 
3800 
1600 
2000 
1300 
40 
380 
310 
560 
3200 
2800 
1300 
14 
1000 
1000 
70 
2000 
70 
60 
50 
1250 
40 
100 
600 
32 
42 
52 
49 
50 
33 
45 
53 
100 
78 
40 
33 
59 


Thickness (ft) 
110 
100 
180 
140 

20 
45 
70 
80 
50 
40 
84 
64 
14 
45 
75 
27 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
80 
84 
95 
115 
150 
40 
100 
80 
70 
120 
60 
100 
120 
90 
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Well Name 


GDC 
HR2 
SR76 
Stuart 
FPL1 
FPL3 


a. Location in Florida East State Plane coordinates, NAD 1983 Datum. 


Well Location? 


X 
8/2536 
882986 
901536 
902668 
808590 
799882 


Y 
1039463 
1045563 
1028163 
1038900 
991248 
995669 


b. Horizontal Hydraulic Conductivity. 


HK” (ft/day) 
62 
97 
33 
33 
53 
48 


Thickness (ft) 


65 
83 
100 
100 
35 
70 


DRAFT Appendix ALEC subRegional MODFLOW Model Documentation 


Draft 


Permit No. 
50-0036 7-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
90-00367-W 
50-00367-W 
50-00367-W 
50-0036 7-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
90-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 


Table 28. 


Well Location? 


X 
944308 
936714 
941056 
943104 
941076 
941101 
941114 
898364 
928988 
946842 
949338 
942220 
947000 
948388 
946280 
945533 
945287 
928988 
904937 
898647 
913388 
911824 
899306 
917423 
918650 
919745 
923041 
898096 
916644 
912936 
879958 
895432 
908903 
916142 
916529 
908266 


Capacity Weil Measurements 


Y 
691738 
749858 
749896 
749595 
749019 
756910 
745522 
713324 
721526 
721909 
719874 
717652 
713129 
712812 
700747 
698272 
695323 
702850 
706711 
695123 
694080 
685101 
679203 
6/8234 
679956 
677628 
677742 
669198 
668756 
661994 
648924 
650594 
653490 
647639 
642859 
633445 





Casing 


Depth (ft) 


110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 


Aquifer Performance Test Results Based on Specific 





Total Depth 

(ft) HK? (ft/day) 
162 2462 
162 329 
162 1469 
162 470 
162 529 
162 1361 
162 1069 
162 100 
162 250 
162 1000 
162 2111 
162 1000 
162 897 
162 700 
162 4400 
162 6364 
162 5500 
162 250 
162 100 
162 292 
162 3800 
162 183 
162 360 
162 1000 
162 1000 
162 1000 
162 1000 
162 1600 
162 4875 
162 4875 
162 2500 
162 4000 
162 4700 
162 1069 
162 1069 
162 7333 


Permit No. 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
06-00135-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-00367-W 
50-0036 7-W 
50-00367-W 

control 
control 
control 
control 
control 
control 
control 
control 
control 
control 
control 
control 
control 
control 
control 


Well Location? 


X 
868318 
884253 
895039 
934547 
934607 
903691 
912426 
924784 
924310 
925553 
931353 
933371 
858214 
846787 
844059 
902761 
848884 
849604 
847753 
888161 
889600 
891040 
855465 
866980 
848833 
826521 
936500 
821228 
844662 
821680 
824396 
818748 
821815 
947287 
946497 
946991 
945442 
945422 
945659 
946478 
947316 


Y 
622154 
627801 
626652 
624071 
623449 
611839 
588974 
597283 
609601 
612734 
604036 
603026 
609072 
617773 
653896 
581776 
553090 
548258 
545071 
543734 
538901 
544865 
494279 
496643 
421895 
416446 
653743 
401306 
394158 
385110 
384985 
385205 
381338 
707342 
707440 
708249 
707440 
708378 
706720 
706602 
706611 
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Casing 


Depth (ft) 


110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 


Total Depth 
(ft) 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
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HK? (ft/day) 
2000 
27500 
1069 
9500 
7368 
15000 
6000 
3421 
25000 
40000 
50000 
7368 
2000 
200 
1633 
6750 
50000 
50000 
50000 
10000 
5000 
10000 
45000 
15000 
5000 
8000 
12000 
75000 
75000 
25000 
25000 
25000 
25000 
125 
125 
125 
125 
125 
125 
125 
125 


Draft 


a. 
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Permit No. 
50-0036 7-W 


control 
control 
control 
control 
control 
control 
control 
calib 
calib 


calib 


Well Location" 


X 
826110 
821806 
814573 
948254 
948496 
949462 
949462 
948408 
903691 
866980 
946280 


Y 
407192 
391748 
388630 
708491 
707744 
707414 
707414 
707129 
611839 
496643 
700747 


Casing 


Depth (ft) 


110 
110 
110 
110 
110 
110 
110 
110 
110 
110 
110 


Total Depth 
(ft) 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 
162 


Location in Florida East State Plane coordinates, NAD 1983 Datum. 


b. Horizontal Hydraulic Conductivity. 


Draft 


HK” (ft/day) 
75000 
75000 
75000 

125 
125 
125 
125 
125 
15000 
5000 
4400 
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Table 30. Conceptual Model layer Properties 


NE | æ | 5 | 
| 24 | 99 | 53 | 
| 20 | 69 | 50 | 
[Car | 2% | 7 | 
| 1 | 5 | 6 | 
| 62 | 68 | 120 | 
| 40 | 15 | 60 | 
| 3 | 6 | 50 | 
| 35 | 105 | 80 
| 35 | 10 | 80 
| 20 | 3z | 9 | 
| 20 | 3z | 9 | 
| 20 | sz | 9 | 
| 20 | 3z | 9 | 
124 

130 


TOP2 


1 1 
1 


: 
O) 


I 
NO 
an 


' l ! 


= 
00 


I 


— 


— 


' I 1 I 1 ' [| 


d 
T 


N 


ES A | 
E NJN 


-29 
-23 
-12 
-16 
-15 
-10 
-10 
-10 
-10 
-24 


-15 
-2 


A-23 


A 


| e9 | -142 | 8281 | 6406 | 10000 | 
| 68. | -124 | 3318 |1431250| 244 ` 
30000.00 
| .68 | -136 | 11605 | 10000.00| 5.00 | 
1722.22 
| -99 | -169 |17717.14| 7600.00 | 38429 | 
21542.31 
| -a | -142 | 4474 |2900000| 101.06 | 
| -9 | -141 | 4730 | 24853 | 690 | 
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| id |THICKL1 BOT3 
| 63305 | 24 | 50 | 10 | -16 | oe | -186 | 7506.25 | 2000.00! 500 ` 
| G3324 | 20 | so | eo | -2 | -e | -142 
6-3315 | 38 | 48 | 90 | -30 
G-3304 45 
42 37 
G-3314 155 EV. 
| G-330 | 17 | 20 | 112 | -ə | aa | -141 | 2144.12 | 29500.00 | 137.59 | 
| G-3312 | 22 | 6 | so | -14 
6-3321 | 32 | 56 | 60 | -24 | -so | -140 | 8129.69 | 20000.00| 33500 | 
(633866 | 26 | 52 | 6 | 18 | -zo | -130 | 2544.23 |20000.00 | 35.00 | 
6-3313 | 36 | so | eo | -23 | - | -168 |2754.17 | 870000 | 3500 | 
G-3307 24 | -98 | -173 | 268594 | 8700.00 | 500.00 | 
| G3300 | 21 | 73 | 80 | -18 | -86 | -166 | 5.00 | 8700.00 | 7625 | 
| G3300 | 5 | 1 | 123 | 3 | -9 [| -13 | 4015.00 | 41292 | 370.73 | 
| LARGO | 32 | 92 | 80 | -24 
mo-177 | 6 | 3 | 79 | 2 [| - | am | 1000 | 1.00 
voam | 6 | 3 | 120 
HE-1110 WEM EN d | 443 | | 8 | 5 | | 138 | | 2500 | I wg. 7 | 21045 | 
HE-1116 5 | 6 | 141 
| P708 | 6 | 6 | 128 | o | -6 | -134 | 4167 | 2500 | 48.91 | 
PB-1704 
| G33055 | 3 | 3 | 128 | 8 | 5 | -123 | 500 | 500 | 20891 | 
| c2226 | 3 | 3 | 128 | 8 | s [| -123 | 450 | 5.00 | 29891. 
6-3308 | 7 | 8 [| 197 
| G-233 | 8 | 9 | 115 | 3 | e | -t | 4750 | 55.00 | 51174 | 
| G2329 | 8 | ə | 115 | 3 | -6 | -2 | 4750 | 5500 | 66826 | 
6-2340 | 8 | 9 | ns | 3 l-6 | -121 | 4750 | 5500 | 43696 | 
| G29144 | 8 | ə | 115 | 3 [| -6 [| -121 | 20313 | 5500 | 24130 | 
| FPLMA | 20 | 60 | 90 | 45 135 
| wi2zees | 20 | 70 | 90 | 15 | -55 | -145 | 4750 | 50.00 | 5000 | 
L-35TW1 -14 
| pp-674 | 32 | 113 | 90 | -15 | -128 | -218 | 1156 | 467.70 
40 13.00 
140 13.00 
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| THICKLA | | TOP2 | Tops | 8019 | HKL1 | HKL2 | ma 

| 4 | -4 | -124 | 1429 | 1000 | 8227 - 

| -8 | -20 | -130 | 1833 | 1000 | 6455 ` 

| 3 [| -æ | -180 | 1333 | 1000 | 8514 ` 
lo | as [| am | tooo | 1000 | 5744 ` 
18 
| ppisio | 20 | 51 | 18 | o | -51 [| -199 | 4500 | 1000 | 65304 - 
120 | 60 | o | -120 | -180 | 1300 | 2500 | 2500 - 
| JpsP | 30 | 90 | 10 | -18 | -108 | -210 | 13.00 | 2500 | 2500 | 
wa | 35 | 60 | 10 | -20 | -so | -180 | 500 | 2500 | 3750 . 
m9 | 22 | 99 | 30 | -12 
vam | 25 | 10 | 9 | z | -117 | -126 | 1300 | 6000 | 5600 | 
| PBi607 | 35 | 99 | 35 | 17 | -me | -151 | 1200 | 2500 | 2500 | 
mos | 24 | 40 | 82 | -2 | -2 | -124 | 3000 | 2500 | 2500 - 
| wt3 | 30 | 60 | 90 | 75 | 75 | -165 | 2500 | 2500 | 36141 | 
mios | 30 | 60 | 90 | -15 
wg | so | 60 | 90 | 35 | 96 | -185 | 2500 | 5000 | 5000 | 
| Mis | 30 | 60 | 90 | ~s | 75 | ~es | 2500 am | 2500 | 
m4 | 30 | 6 | 90 | 15 | -~s | -165 | 2500 | 3333 | 2500 ` 
| 01091 | 30 | 60 | 90 | 185 | 75 | ~es | 500 | 4167 | 5000 | 
Peace | 4o | so | eo | -24 [| -z4 | -134 | 0.00 | 2000 | 2000 ` 
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